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1 Introduction 



This review is devoted to the issue of hidden integrabiUty in quantum field theory, i.e. integrability in 
the systems that naively have not enough many conserved integrals of motion. This situation differs 
considerably from the picture familiar for the experts in the theory of integrable systems. In the typical 
integrable system, one deals with a well-defined (field theory) Lagrangian, with the fields given on a 
base space-time manifold. The fields at some fixed moment of time, i.e. given on a space slice of the 
base manifold, provide the dynamical degrees of freedom, while integrable flows are identified with 
the time direction in the base manifold. 

On the contrary, in the systems we are interested in, the phase space for the integrable dynamics 
is provided by a moduli space essential in the problem. The collective coordinates or moduli arise in 
the course of consideration of various issues and have essentially different origins. In particular, they 
can describe effective degrees of freedom important in one or another regime of the field theory not 
completely integrable by itself. The proper regime can be determined by the kinematics of the process 
or by the choice of specific topologically nontrivial data. In the cases when there are several collective 
degrees of freedom, a hidden integrability in the theory emerges as a rule. Its origin lies in symmetry 
properties of the moduli space. 

As a particular example, it appears that a hidden integrability governs the low-energy effective 
actions of A/" = 2 supersymmetric gauge theories, which in four dimensions were described by Seiberg 
and Witten |l|, |. 

The role of this integrability is twofold; it allows one to uncover the hidden symmetries and 
to define proper degrees of freedom. In what follows, we shall show how the low-energy sector of 
supersymmetric gauge theories is described via finite dimensional integrable systems. 

Restricted frames of the review does not allow us to provide the reader with all necessary infor- 
mation from the theory of integrable systems, supersymmetric gauge theories and nonperturbative 
string theory which might be relevant to the issue. This kind of problem arises since the integrability 
appears at the intersection of different branches of theoretical physics which use different tools and 
approaches. Therefore, we just recommend several (rather accidental) review papers, where the issues 
used in what follows as well as the proper list of references can be found. Other references can be 
found in the proper parts of the present review. 

The finite gap integration is perfectly presented in , Whitham hierarchies are discussed in |Q, || , 
while background of the Seiberg- Witten approach is described in |^ . Nonperturbative string theory is 
reviewed in ^ . Necessary background from algebraic geometry used in the brane approach is described 
in while a generic point of view on integrability in (nonperturbative) field/string theory can be 
found in |^]. An intensive discussion of integrability in Seiberg- Witten theory and related structures 
(Whitham hierarchies and WDVV equations) can be found in the book that contains recent 
comprehensive reviews. Note also that the very detailed review may serve a good complimentary 
text to the present paper. 

1.1 Moduli and nonperturbative degrees of freedom 

Description of the strong coupling regime remains one of the most important problems of quantum 
field theory. Since the perturbation theory does not work in this region, one can mostly reply on 
constructing some effective degrees of freedom which would allow one to develop a new perturbation 
theory with a small parameter distinct from the coupling constant of the initial theory. A nice example 



of the existence of different descriptions at different energy scales is quantum chromodynamics (QCD). 
At high energies, quarks and gluons giving the proper degrees of freedom are described by the Yang- 
Mills Lagrangian coupled to fermions in the fundamental representation, while at low energies it is a 
chiral (cr-model) Lagrangian that describes the low-energy effective action, with the proper degrees of 
freedom given by the colorless mesons and baryons. Note that the parameters of the chiral Lagrangian 
can not be calculated from the first principles (at the moment) and are fixed phenomenologically. 

The general structure of low-energy actions is fixed by the symmetries. It is expected that the 
effective actions enjoy a complete set of the Ward identities of the initial field theory. Some of 
the Ward identities may be anomalous. However, since the anomalies are related to the level crossing 
phenomena, not renormalized and, therefore, can be equally well described in the UV and IR terms [ p^ , 
the anomalous Ward identities could still impose restrictions on the structure of effective theories. For 
instance, the chiral anomalous Ward identities fix the Lagrangian of mesons, while the conformal Ward 
identities yield the effective Lagrangian for the dilaton in ordinary as well as in A/" = 1 supersymmetric 
QCD 

Another phenomenon that has much to do with the issues discussed is universality. Effective 
actions typically enjoy this feature which implies that several theories different in the UV regime can 
fiow to the same theory in the infrared region. The reason for such an universality to occur is the 
rigidity of effective actions due to the symmetries. Moreover, the symmetry in the IR regime may be 
even higher then that in the initial UV theory. 

The symmetry nature of effective actions results in nontrivial consequences. Presence of an un- 
derlying integrable systems seems to be the most important of all. As a rule, the partition function 
calculated in the effective theory appears to be the r-function of some integrable system. On the 
other hand, the r-function is a generating function for the conservation laws in integrable system and 



simultaneously the generating function of correlators in some quantum field theory |9|, 14 1. 

Identification of the variables in the integrable system describing the effective action is a difficult 
problem. An invariant way to introduce these variables in non-topological theories is still missing. 
However, there are several examples of two-dimensional theories when these variables can be identified 
with the transition amplitudes between vacuum states induced by nonperturbative configurations. 
The couplings and sources serve as the "space-time variables" in integrable systems. Examples of 
the systems related to the 2d Toda lattice and similar hierarchies (KP, KdV etc) can be borrowed 
from 2d tr-models 15 1 and 2d gravity ||l6|, 17, |l^, |l^, ||]. Selecting the concrete solution to the 
nonlinear integrable equations is a separate problem. It is typically done by an additional equation, 
string equation [p^P| which is one of the Ward identities imposed on the theory ^ . 

Although one can fix the general structure of effective action from the symmetry arguments, 
it is also interesting to evaluate it by the direct integration over heavy modes in the initial field 
theory. It is implied that there is some energy scale and all modes above this scale are integrated 
out. This integration should be done over the perturbative configurations and the scale serves a cutoff 
in the corresponding Feynman loop diagrams. It is also necessary to sum over the nonperturbative 
configurations, their typical sizes being restricted by this scale. The nonperturbative configurations 
important in different dimensions are: instantons, monopoles, vortexes as well as solitons of different 
codimension in string theory. 

The important feature of nonperturbative effects is that they give rise to the moduli spaces of the 
corresponding solutions. The very reason for the moduli space to emerge is the high symmetry of the 

^ The string equation does not fix the possible discrete ambiguity of solutions . 



theory. Therefore, the moduli space itself encodes symmetry properties. For instance, the dimension 
of the instanton moduli space is fixed by the space-time symmetries - translations and rotations 
as well as by the gauge symmetries |23| There are many other important moduli spaces, for 

instance, the moduli spaces: of the complex structures of Riemann surfaces, of the flat connections on 
a fixed surface, of the holomorphic vector bundles etc. Each of these appears when considering some 
nonperturbative configurations and can be derived via the reduction procedure from the instanton 
moduli space. Thus, the universal instanton moduli space comes as a most general object arising in 
field theories in four or less dimensions. 

As we explained, the integration over moduli spaces involved into the derivation of effective actions 
should ultimately give rise to integrable theories. The same kind of moduli spaces have much to do 
with integrable systems also from another point of view. Namely, the phase spaces of integrable 
systems typically coincide with moduli spaces or cotangent bundles to the moduli spaces of the same 
type. Simple examples are given by solutions of the KdV hierarchy related to the moduli space of the 
complex structures or by solutions of the 2d Toda lattice hierarchy related to the moduli of the flat 
connections. Another important example is given by the Chern-Simons theory defined on a surface, 
whose phase space is the moduli space of the flat connections on this surface [ p5[ , while any gauge 
invariant observable can serve as a Hamiltonian. In terms of integrable systems, the problem of 
calculation of the nonperturbative contributions into effective actions is nothing but the calculations 
of some vacuum expectation values of observables on moduli spaces (see as an example). 

A large class of theories related with the moduli spaces is known as topological theories. They were 
introduced into quantum field theory by Witten j2^] to obtain invariants of the vacuum manifolds in 
field theories. The first examples of topological theories have been elaborated in four dimensions for 
a-models and gauge theories. In the cr-model case, the topological theory provides invariants of the 
holomorphic maps, while in the gauge theory it gives invariants of the instanton moduli space |2^. At 
the moment, there are a lot of examples of topological theories in different dimensions, e.g., tr-models 
in d = 2, Chern-Simons theory in d = 3, topological Lagrangians in d = 4 etc. 

In all the cases, the topological theories can be obtained from the nontopological Af — 2 SUSY 
theories via the twisting procedure |p9|| . The partition function of the topological theories appears to 
be some invariant of the corresponding moduli space. 

There are also some examples of field theories with the topological subsectors related solely with 
the vacuum states. An example is given by the J\f ^ 1 supersymmetric gauge theory, where the 
perturbative /3 function js^ as well as some correlators can be calculated exactly js^. At last, let 
us mention that the associativity (WDVV) equations which typically emerge in topological theories 
also appear in Af = 2 SUSY theories. This implies some hidden topological structures (subsectors) in 
these theories as well 

Since the moduli spaces arising in different applications are typically finite dimensional, the topo- 
logical theories that provide finite dimensional moduli spaces are of the prime interest. It is these 
moduli spaces where the dynamics of integrable many-body systems with nontrivial conservation laws 
is developing. Gauge invariant observables with appropriate couplings serve the Hamiltonians. 

Note that the partition function for the perturbed topological theory, i.e. the generating function 
for the correlators is the most interesting and informative object. It is simultaneously the generating 
function for topological invariants of the moduli space manifold [ p7[ and the r-function of an integrable 
hierarchy |§. 



1.2 Seiberg-Witten solution 



New possibilities to deal with nonperturbative physics at strong coupling come from the exact deriva- 
tion of the low energy effective action for the Af = 2 SYM theory ||. In principle, to get the exact 
answer one needs to perform the summation of the infinite instanton series. Therefore, the result |^ 
due to N.Seiberg and E.Witten was the first example of the exact derivation of the total instanton 
sum. 

However, the authors did not calculate instantonic series, using instead indirect and quite sophis- 
ticated methods. The calculation of the exact mass spectrum of stable particles was their second 
essentially new result. What is important, this theory is non-trivial in one loop and belongs to the 
class of asymptotically free theories. In fact, the situation is a little bit unusual, since, despite the fact 
that the direct instanton summation has not been performed yet, there are no doubts in the results 
derived which survives under different tests. 



In order to get the effective action unambiguously, the authors of [gj, |^ used the three ideas: holo- 
morphy | |3^ , duality and their compatibility with the renormalization group flows. Holomorphy 
implies that the low energy effective action can be only of the form J hnJ^i^) depending on a single 
holomorphic function J- called prepotential. The duality principle which provides the relation between 
the strong and week coupling regimes is natural in the UV finite^ say A/" = 4 SUSY , theory, where 
the evident modular parameter built from the coupling constant and the 0-term 

Ani ,^ ^. 
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is not renormalized and coincides with its asymptotic value. To construct duality transformation 
in asymptotically free theory explicitly, one needs an artificial trick, to introduce into the game an 
additional object, a higher genus Riemann surface whose period matrix coincides with the matrix of 
couplings depending on the values of condensate in a given vacuum. 

In the Af = 2 SYM theory this scheme is roughly realized as follows First of all, in the 

Lagrangian there is a potential term for the scalar fields of the form 



\^(0)-Tr[(/),0+]2, (1.2) 

where the trace is taken over adjoint representation of the gauge group (say, SU{Nc)). This potential 
gives rise to the valleys in the theory, when [(/>, 0+] = 0. The vacuum energy vanishes along the valleys, 
hence, the supersymmetry remains unbroken. One may always choose the v.e.v.'s of the scalar field 
to lie in the Cartan subalgebra </> ~ diag{ai, Un). These parameters can not serve, however, 
good order parameters, since there is still a residual Weyl symmetry which changes but leave 
the same vacuum state. Hence, one should consider the set of the gauge invariant order parameters 
Uk =< Trcf)'' > that fix the vacuum state unambiguously. Thus, we obtain a moduli space parametrized 
by the vacuum expectation values of the scalars, which is known as the Coulomb branch of the whole 
moduli space of the theory. The choice of the point on the Coulomb branch is equivalent to the choice 
of the vacuum state and simultaneously yields the scale which the coupling is frozen on. At the generic 
point of the moduli space, the theory becomes effectively Abelian after the condensation of the scalar. 

As soon as the scalar field acquires the vacuum expectation value, the standard Higgs mechanism 
works and there emerge heavy gauge bosons at large values of the vacuum condensate. To derive 
the effective low energy action one has to sum up the loop corrections as well as the multi-instanton 

■^In fact, one should require more: not renormalized and, therefore, non-running coupling constants. 



contributions. The loop corrections are trivial, since they are vanishing beyond the one loop due to the 
supersymmetry. The explicit summation over instantons has not been done yet, however, additional 
arguments allow one to define the effective action indirectly. 

The initial action oi M — 2 theory written in A/" = 2 superfields has the simple structure S{'^) = 
Imr J Tr^*^ with r = ^ + ^ . Now the N — 2 supersymmetry implies that the low energy effective 
action gets renormalized only by holomorphic contributions so that it is ultimately given by a single 
function known as prepotential S'e//(5') = Im J The prepotential is a holomorphic function of 
moduli Uk (or v.e.v.'s) except for possible singular points at the values of moduli where additional 
massless states can appear disturbing the low energy behaviour. 

Thus, the problem effectively reduces to the determination of one holomorphic function. If one 
manage to fix its behaviour nearby singularities, the function can be unambiguously restored. One of 
the singularities, corresponding to large values of v.e.v.'s, i.e. to the perturbative limit is under control 
(since the theory is asymptotically free). All other singularities are treated with the use of duality 
and of the non-renormalization theorems for the central charges of the SUSY algebra. A combination 
of these two ideas allows one to predict the spectrum of the stable BPS states which become massless 
in the deep nonpcrturbative region and are in charge of all other singularities. 

The duality transformation can be easily defined in the finite M = A SUSY theory just as the 
modular transformations generated by r — > and t — > r + 1. This makes a strong hint that the 
duality can be related with a modular space of some Riemann surfaces, where the modular group 
acts. In fact, the naive application of duality meets some serious difficulties. The reason is that, 
in the asymptotically free theory, one has to match the duality with the renormalization group^. 
This is non-trivial, since now r depends on the scale which is supposedly involved into the duality 
transformation. The solution of this problem is that one is still able to connect the duality and 
modular transformations if considering the theory at different vacua connected by the duality. Then, 
the duality acts on the moduli space of vacua and this moduli space is associated with the moduli 
space of the auxiliary Riemann surface, where the modular transformations act. 

At the next step, one has to find out proper variables whose modular properties fit the field theory 
interpretation. These variables are the integrals^ of a mcromorphic 1-form dS over the cycles on the 
Riemann surface, a, and a? 



dS, 

(1-3) 

aB = <f dS, 

(where i,j — 1, iVc — 1 for the gauge group SU{Nc)). 

^In fact, there are more problems with duahty even in A/" = 2 SUSY theories. In particular, it unifies into one duality 
multiplet the monopole and gauge boson supcrmultiplets, while they have different spins, see, e.g., ||^, |6j. 
■* We define the symbols ^ and res with additional factors {27rj)~^ so that 

reso— = -resoo— ^ J ~^ 

This explains the appearance of 2-iTi factors in the Riemann identities and in perturbative prepotentials in sections 5 
and 6. Accordingly, the theta-functions are periodic with period 2-Ki, and 

— — — = i-Kdije(i,\r) 

aiij 

since periods of the Jacobi transformation = dia^i belong to 27ri (2 + ^Z)- 
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These integrals play the two-fold role in the Seiberg-Witten approach. First of all, one may 
calculate the prepotential T and, therefore, the low energy effective action through the identification 
of and dT /da with a defined as a function of moduli (values of condensate) by formula (1.3). 
Then, using the property of the differential dS that its variations w.r.t. moduli are holomorphic one 
may also calculate the matrix of coupling constants 



daidan 



(1.4) 



The second role of formula (1.3) is that, as was shown these integrals define the spectrum of 
the stable states in the theory which saturate the Bogomolny-Prasad-Sommerfeld (BPS) limit. For 
instance, the formula for the BPS spectrum in the SU (2) theory reads as 



A^n,m = \na{u) + maD{u)\ , (1.5) 

where the quantum numbers n, m correspond to the "electric" and "magnetic" states. The reader 
should not be confused with the spectrum derived from the low energy behaviour which fixes arbitrarily 
heavy BPS states. The point is that the BPS spectrum is related to the central charge of the extended 
SUSY algebra ^'j^j {„} = a^n' +afm'^ and, therefore, has an anomaly origin. On the other hand, 
anomalies are not renormalized by the quantum corrections and can be evaluated in either of the UV 
and IR regions. 

Note that the column (ai,af) transforms under the action of the modular group SL{2,'^) as a 
section of the linear bundle. Its global behaviour, in particular, the structure of the singularities is 
uniquely determined by the monodromy data. As we discussed earlier, the duality transformation 
connects different singular points. In particular, it interchanges "electric" , and "magnetic" , 
variables which describe the perturbative degrees of freedom at the strong and weak coupling regimes 
(of the theory at different vacua) respectively. Manifest calculations with the Riemann surface allow 
one to analyze the monodromy properties of dual variables when moving in the space of the order 
parameters. For instance, in the simplest SU{2) case, on the u-plane of the single order parameter 
there are three singular points, and the magnetic and electric variables mix when encircling these 
points. Physically, in the theory with non- vanishing 0-term (not pure imaginary r), the monopole 
acquires the electric charge, while the polarization of the instanton medium yields the induced dyons. 



1.3 Seiberg-Witten theory and theory of prepotential 

Although the both auxiliary objects, the Riemann surface and the differential dS have come artificially, 
an attempt to recognize them in SUSY gauge theories results into discovery of the integrable structures 
responsible for the Seiberg-Witten solutions. 

Now let us briefly formulate the structures underlying Seiberg-Witten theory. Later on, we often 
refer to Seiberg-Witten theory as to the following set of data: 

• Riemann surface C 

• moduli space Ai (of the curves C), the moduli space of vacua of the gauge theory 

• meromorphic 1-form dS on C 

How it was pointed out in p^ , this input can be naturally described in the framework of some 
underlying integrable system. 



To this end, first, we introduce bare spectral curve E that is torus — + g2X^ + 93 for the 
UV-finiteQ gauge theories with the associated holomorphic 1-form duj — dx/y. This bare spectral curve 
degenerates into the double-punctured sphere (annulus) for the asymptotically free theories (where 
dimensional transmutation occurs): a; — > w + y — s- w — l/w, duj — dw/w. On this bare curve, 
there are given a matrix-valued Lax operator L{x,y). The corresponding dressed spectral curve C is 
defined from the formula det(i — A) = 0. 

This spectral curve is a ramified covering of E given by the equation 



P(A;x,y)-0 (1-6) 
In the case of the gauge group G = SU{Nc), the function T-" is a polynomial of degree Nc in A. 

Thus, we have the spectral curve C, the moduli space A4 of the spectral curve being given just 
by coefficients of V. The third important ingredient of the construction is the generating 1-form 
dS = Xdco meromorphic on C ("=" denotes the equality modulo total derivatives). From the point of 
view of the integrable system, it is just the shortened action "pdg" along the non-contractablc contours 
on the Hamiltonian tori. This means that the variables at in (1.3) are nothing but the action variables 
in the integrable system. The defining property of dS is that its derivatives with respect to the moduli 
(ramification points) are holomorphic differentials on the spectral curve. This, in particular, means 
that 



ddS , ,^ 
— = dw^ (1.7) 

where du^i are the canonical holomorphic differentials]^ Integrating this formula over B-cycles and 
using that au ~ dT jda^ one immediately obtains ( |l.4|) . 

So far we reckoned without matter hypermultiplets. In order to include them, one just needs to 
consider the surface C with punctures. Then, the masses are proportional to the residues of ^5* at the 
punctures, and the moduli space has to be extended to include these mass moduli. All other formulas 
remain in essence the same (see |^ for more details). 

Note that it is still unknown whether the set of data above is identically equivalent to an integrable 
system. Inversely, it is associated only with specific integrable systems, which can be called analytic 
integrable systems. Although being wide enough, this class of integrable systems far docs not cover 
the systems integrable in the Liouvillc sense. 

The prepotential T and other "physical" quantities are defined in terms of the cohomology class 
of dS*, formula (l.S). Note that formula ( |l.4| ) allows one to identify the prepotential with logarithm 
of the T- function of the Whitham hierarchy |Q: T = logTwhitham- The other important property of 
the prepotential related to the feature of the Seiberg-Witten framework which merits the adjective 

^The situation is still unclear in application to the case of fundamental matter with the number of matter hyper- 
multiplets Nf equal to 27Vc. In existing formulation for spin chains the bare coupling constant appears rather as a 
twist in gluing the ends of the chain together (this parameter occurs only when Nj = 2Nc) and is not immediately 
identified as a modulus of a bare elliptic curve. This problem is a fragment of a more general puzzle: spin chains have 
not been described as Hitchin systems; only the "2 X 2" Lax representation is known for them, while its "dual" x 
one is not yet available. 

^I.e. satisfying the conditions 
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"topological" has much to do with the associative algebras Q. Namely, it turns out that the pre- 
potential of Seiberg-Witten theory satisfies a set of the Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) 
equations. These equations have been originally presented in (in a different form, see below) 



F,F-^Fk = FkF-^F, (1.8) 

where i^i's are matrices with the matrix elements that are the third derivatives of a unique function 
F of many variables Oj's (prepotential in the Seiberg-Witten theory) parameterizing a moduli space: 



Although generally there are a lot of solutions to the matrix equations (1.8), it is extremely non- 
trivial task to express all the matrix elements through the only function F. In fact, there have been 
only known the two different classes of non-trivial solutions to the WDVV equations, both being 
intimately related to the two-dimensional topological theories of type A (quantum cohomologies [p2|) 
and of type B {J\f = 2 SUSY Ginzburg-Landau theories that were investigated in a variety of papers, 
see, for example, and references therein). Thus, the existence of a new class of solutions connected 
with the four-dimensional (and higher dimensional) theories looks quite striking. 



1.4 Duality, branes and gauge theories 

The progress in the description of the low energy SUSY gauge theories initiated the investigation 
of duality in string theory. It appeared that, in order to formulate the duality transformation in 
string theory, one has to introduce new objects called D-branes Their main property is that the 
open strings can end on the D-brane. Therefore, the Dirichlet boundary condition can be imposed. 
Moreover, the charge in the RR sector of the theory is attributed to the D-brane. Using D-branes, a 
stringy duality can be self-consistently described, and actually three different types of dualities can 
be introduced Q. 

Since the superstring theory is anomaly free only in ten dimensions, the question about the role 
of remaining dimensions naturally arises. The results of the last years definitely indicate that these 
dimensions should be identified with the zero modes of the scalar fields presented in 4d field 
theories. In the most symmetric situation in d = 4, which is the J\f — A SUSY gauge theory there 
are 3 complex scalars whose vacuum expectation values should be considered as the six compactified 
dimensions. Generically the structure of the field space is very complicated and its metric is the 
subject of numerous investigations. 

From the field theory point of view. A/" = 4 is maximal admissible supersymmetry, since it is the 
case when the gravitational degrees of freedom are still decoupled |Q from the gauge sector of the 
theory, the case of our prime interest. 

Now, in order to analyze a generic Ad gauge theory, one has to introduce an arbitrary gauge group, 
to ask for the proper SUSY and to define all the parameters in the theory such as masses and couplings. 
It appears that the account of the brane degrees of freedom allows one to fulfill this program at full. 
The important properties of branes which allow one to obtain these results are the presence of an 
f7(l) gauge field on their worldvolumes as well the existence of brane bound states ]4^ . 

Let us briefly describe how the 4d field theory can be formulated as a theory on the brane world 
volume of some dimension. At the first step, we start with the superstring theory in ten dimensions 
of type II, and branes are the BPS states breaking half of the supersymmetries (this gives us A/" = 

1 n 



4 SUSY in Ad). In one approach, the initial theory is considered in the flat space and SUSY is 
broken by the background brane configuration ]4^ . In the other approach, known as "geometrical 
engineering" [^8| the supersymnietry is broken by a special choice of metric of the six dimensional 
compactification manifold. For instance, the metric of the K3 manifold gives rise to breaking of the 
half of supersymmetries. 

As a corollary of definition, the strings ending on the D-brane induce an U{1) gauge field living on 
its worldvolume. At the next step, one should produce the general gauge group, say SU{Nc), which 
is absent in the original string theory. Once again two approaches are possible; in the first approach, 
the gauge group arises from the bound state of branes |Q , while in the second approach it originates 
from choice of the compactification manifold with the singularity structure related to the gauge group 
root system p9| . The higher rank gauge group obtained from the brane bound state can be explained 
as follows. Each of the branes carries the C/(l) index while the open strings can stretch between 
them. Therefore, having the stack of parallel branes, each string can be associated with an U{Nc) 
matrix realizing the U{Nc) gauge theory^. If one now splits the branes, there are massive modes of 
the gauge fields in the spectrum of open strings which correspond to the strings stretched between 
different branes and the corresponding masses are proportional to the distance between branes in 
the embedding space. This is a realization of the Higgs mechanism. When the distance between 
branes vanishes, rank of the corresponding "unbroken" gauge group increases. In the "geometrical 
engineering" approach the gauge group comes from the wrapping of branes over the manifold with 
the corresponding singularity structure [^0| . 

Note that the scalars in field theory play the role of coordinates of the brane or the bound state 
of branes in the 6d space. 

The next step of this rough (rather pictorial) construction of field theory is to introduce massive 
hypermultiplets and couplings. To this end, one should further break the supersymmetry]^ It can be 
done by adding more, background branes, of different codimensions and non-parallel to the stack of 
Nc branes. 

Say, one can introduce hypermultiplets by adding branes of other codimensions so that the coor- 
dinates of these branes in 6d compactification manifold correspond to the values of their masses . 
Couplings are introduced similarly and their values are identified with the distance between the back- 
ground branes in one picture or by the size of the compact manifold in the "geometrical engineering" 
approach |48 . 



1.5 Seiberg-Witten solution and branes 

Realization of the Seiberg-Witten solution of the Af = 2 SUSY gauge theory in brane terms seems 
to be the most impressive realization of the approach described. In the type IIA framework, the 
theory is solved along the following line |5^. The SU{Nc) gauge theory arises on the worldvolume 
of Nc parallel D4 branes embedded into the flat lOd manifold. The worldvolume of the branes is five 
dimensional and stretches along, say, xo,xi, a;2, xa, xe-directions, the branes being of finite size in the 
a;4-direction, with the first four coordinates identified with the coordinates in our (infinite) 4d world. 

Now one has to describe the compact six dimensions in field theory terms. To this end, one first 

^It can be shown that the total U{1) factor decouples. 

®Note that, due to the conformal invariance, the = 4 SUSY theory depends on the coupling constant trivially and 
can not contain massive hypermultiplets. 
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assumes that all the branes are placed at the origin with respect to three coordinates^, say, xy, ss, s^g- 
Now, the two remaining directions not stretching along the brane, Xi and X5 form the holomorphic 
variable, X4 + 1x5 which is associated with the vacuum expectation value of the scalar in the J\f — 2 
SUSY gauge theory, that is to say, the coordinate on the Coulomb branch of the moduli space. At last, 
the compact coordinate of brane, xq is just the coordinate on (real section of) the auxiliary Riemann 
surface. 

Since the D4 branes are restricted along the xg-direction, one needs some other, background 
branes where the D4 branes should end on. It is to be two additional branes of the NS type. It 
is just this configuration which breaks the Af ~ 8 SUSY in ten dimensions (32 generators) down to 
the Af — 2 SUSY in four dimensions (8 generators). The distance between the background branes 
along the a;6-direction is identified with the coupling constant in the gauge theory Sxq = p- , while the 
back reaction to the background brane can be attributed to the perturbative renormalization of the 
coupling. In order to incorporate the fundamental matter hypermultiplets into the theory, one has to 
introduce additional Nf D6 branes whose coordinated^ along X4 + 1x5 are related with masses of the 
hypermultiplets. On the other hand, in order to add the adjoint matter, the (complexified) coordinate 
xq should be made periodic. 

The configuration considered is singular, since it is impossible to describe the intersection of the 
D4 and background NS5 branes smoothly. To resolve the singularities, Witten has suggested to lift the 
configuration into d = 11 manifold which plays the same role for membrane theory as the d — 10 space 
for superstrings. The theory defined in d = 11 is known as M-theory and it has the c? = 11 supergravity 
as its low energy limit |5^. It is important that there are solitons in M-theory, that is, M5 branes 
and membranes (M2 branes). In the non-singular M-theory picture of A/" = 2 SUSY gauge theories, 
the eleventh dimension is assumed to be compact and there is a single M5 brane wrapped around 
it. This M5 brane is assumed to be topologically non-trivial and to have 2 compact and 4 infinite 
dimensions, these latter being associated with our 4d world. With three coordinates {xj, xs,xg) being 
set to zero, the two compact dimensions of the M5 brane are embedded into the four-dimensional 
space (of X4,X5,xe,xio coordinates) in a tricky way. This embedding is, however, holomorphic and 
describes the auxiliary Riemann surface. Being projected onto ten dimensions, different pieces of this 
two-dimensional surface are associated with the D4 and NS5 branes. 

Note now that a similar picture has been also developed in the M-theory framework for J\f — 1 
gauge theories [ p3| which manifest confinement due to the monopole condensation. This gives some 
hope for relevance of the brane picture also for the ordinary, non-supersymmetric QCD. 

In the other IIB approach [54[ | related to the previous consideration via the T duality transformation 
, the embedding manifold is topologically non-trivial and branes are assumed to be wrapped around 
some cycles on the surfaces, for instance, on tori (elliptic curves) embedded into K3 manifold @. 
This approach is less transparent, although the topological theories on the curves which are associated 
with the intcgrablc many-body systems arise within this framework automatically. Similarly to the 
way of resolving the singular IIA configuration by lifting it into higher dimensional M-theory, the 
IIB configuration has to be lifted to the d = 12 theory known as F-theory |5^. The additional two 
dimensions are assumed to be compactified onto the torus whose modular parameter is identified with 
the (imaginary and real parts of) UV values of the coupling constant and the 0-term in the gauge 
theory. 



''This is necessary to preserve, at least, JV = 2 SUSY. 
''More precisely, the distance to the D4 branes. 
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2 Many-body systems and gauge theories in two and three 
dimensions 



In this section we demonstrate how the many-body systems naturally arise in two and three dimen- 
sional gauge theories and the phase space of the system of particles is interpreted in terms of the 
topological gauge degrees of freedom. From what follows it becomes clear that the 2d YM theory 



is closely related to the Calogero-Moser system |5q, 59, 60 , while the Chcrn-Simons theory to its 



relativistic counterpart the Ruijsenaars-Schneider system. 

2.1 Calogero many-body system and Yang-Mills theory in two dimensions 

Let us show how the Lagrangian formulation of the Calogero systems can be related to the YM theory 
on a cylinder lo. Remind that the Calogero system is defined as an N particle system with the 



pairwise interaction 58, p9| 



^^^^ (2.1) 

v"^ is the coupling constant. As we already mentioned, the phase spaces of integrable systems are 
typically related with group manifolds. Therefore, it should not be surprising that the Lagrangians 
under consideration have a group origin. 

Let us start with the central extension g of the loop algebra £g, where g is a semi-simple Lie 
algebra. It is defined as the space of pairs (0, c), where (f) is the g- valued scalar field given on the circle 
and c is the central element, just a real number. The Lie algebra structure is introduced by formula 

[((/.i(^),ci),(02(^),C2)] = (^[M^),cb2{v)]J^^<^i,d^h >) , (2.2) 

where < ., . > denotes the Killing form on the algebra g. 

Now we set the cotangent bundle to this algebra, T*'g to be a "bare" phase space. T*g is defined 
to be a direct sum of g and g*, i.e. the space of sets {A, k; <j>, c), where A is the gauge field on the 
circle and k is dual to c. The pairing between g and g* is given by 

< {A,K);{(f>,c) >^ <(f),A>+cK, (2.3) 
On this cotangent bundle, one can define the natural symplectic structure 



fi= / Tr{S(l)ASA)+5cASK (2.4) 
The adjoint and coadjoint actions of the loop group CG on g and g* are given by formulas 

-1 f rr.,./ J,„-1; 



(0(^),c) -> i^9iip)H^)gi^r\J^^ Tr(-0.g-^a^.g) + cj (2.5) 

{A, k) {gAg-^ + ngd^g-^n) (2.6) 
This action clearly preserves the symplectic structure and thus defines a moment map 

^l■.T*g^g\ (2.7) 
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which sends {<f), c; A, n) to (nd^cj) + [^7 0]; 0). 

Note that the choice of level of the moment map has a simple physical meaning: in the theory of 
the angular momentum it just fixes the sector of the Hilbert space. In the gauge theory, one selects 
the gauge invariant states. 

Let us now discuss how to choose the proper level of the moment map in order to get further the 
Hamiltonian reduction w.r.t. this level. Let us remind that in the case of finite-dimensional group G 
one takes an element J from g*, which has the maximal stabilizer different from the whole G. It is 
easy to show that the representative of the coadjoint orbit of this element has the following form 

Ju = v ^ iec + e-a), (2.8) 

where v is some real number, e±a are the elements of the nilpotent subalgebras, n± C g corresponding 
to the root a, and A+ is the set of positive roots. Let us denote the coadjoint orbit of J^, as O^, and 
the representation of G arising upon the quantization of O^, as R„. For a generic J S g*, let Gj be 
the stabilizer of J in the coadjoint orbit Oj, thus Oj = G/Gj. 

Now in order to deal with the affine case, one just needs to substitutes the roots by their afRne 
counterparts. There is also a more geometrical standpoint: because of the vanishing level k =p of 
the moment map, the coadjoint orbit of its generic value is huge, while the stabilizer of the element 
{J{(p);0) S g* is "a continuous product" Ili^gsi ^Jiv) which is very small as compared with the whole 
loop group CG unless the group Gj(;p) coincides with G for almost all 4> G S^. 

All this implies that one can reasonably choose the value of the moment map to be 

fi = {J[fi],0) : J[fi]{^) ^ S{^)J, (2.9) 
The coadjoint orbit of fi is nothing but the finite dimensional orbit O^. 

Now in order to perform the Hamiltonian reduction, one has to resolve the moment mapP| 

nd^cj) + [A, 0] = J^((^) = J(^) J,. (2.10) 

It can be done in the following way. First, we use the generic gauge transformation g{ip) to make A 
to be a Cartan subalgebra valued constant one-form D. However, we are still left with the freedom 
to use the constant gauge transformations which do not touch D. Actually the choice of D is not 
unique and is parametrized by the conjugacy classes of the monodromy exp{^D) e T C G. Let 
us fix a conjugacy class and denote ixi the elements of the matrix D = iQ. Decompose g-valued 
function (j) on the on the Cartan- valued Pif) G t and nilpotent-valued 4'±{ip) G n± components. 



Let (j)a =< 4',ea >■ Then ( p.lCj ) reads: 

Kd^P ^ Siip)[J^]^ (2.11) 

Kd^^^+<D,q}^>=S{^)[JS]^, (2.12) 
where is AdX^^^{J^), [J].y denotes the Cartan part of J and [J]a —< J,ea >■ 



From (2.11) one gets D = const and [J^]^ = 0. Therefore, one can twist back to J^. Moreover, 



( 2.12 ) implies that {if ip ^ ) 4>a{^) can be presented as 



^^This is nothing but the Gauss law. 



(f>aiip) = exp (^-^ < D,a >j X Ma, 



(2.13) 



where Ma is a locally constant element in g. From ( 2.12 ) one sees that Ma jumps, when (p passes 
through 0. The jump is equal to 



, 2tt 

exp < D,a > 



Ma = [J^]a. 



(2.14) 



The final conclusion is that the physical degrees of freedom are contained in exp{—^Q) and P 
with the symplectic structure 



and one has 



to = —Ti(SP A SQ) 
2tti 



(2.15) 



exp(-^ <Q,a>) 



exp(-2^ < Q,a >) - 1 

With this manifest expression, one easily restores the whole (pif) adding the constant Cartan part. 
Note that, while the diagonal elements qi of Q determine particle coordinates, the constant Cartan 
part which can not be fixed from the moment map is parameterized by particle momenta pi. In these 
variables, (f>{(p) is nothing but the Calogcro Lax matrix and the particles are the Calogero ones. 

Now if one takes a simple Hamiltonian system on T*g with a Hamiltonian invariant under the 



action (2. 5), (2. 6), this induces a relatively complicated system on the reduced symplectic manifold, 
T*T. As a Hamiltonian on T*g, one can choose 



n= I dipP{(j)) 



(2.16) 



where P((/)) is any invariant polynomial on the Lie algebra g. This is, indeed, the standard way to 
construct Hamiltonians from the Lax matrix in integrable system. 

Let us choose the Hamiltonian to be the quadratic Casimir 



H2 = I dip < 41,(1) >, 
'51 



(2.17) 



and look at the action of this system. Taking into account the constraint, i.e. adding to the action 
the moment map with the Lagrangian multiplier, one obtains 



S 



dipdt{TvAdt(l)- <<f),<j)> --Trylo(Ka^0 + [A, cj)] - 5{ip)J^)). 



Denoting Ai = A and gvAi = l/'s, one can rewrite this action in the form 



S 



d^x 



-Tr {4>Fai 



(2.18) 



(2.19) 



where Fqi is the component of the gauge field strength tensor and C is the contour passing through 
ip = along the time axis. This action of induced gauge theory, after integrating over the scalar field 
(/) can be rewritten as the action of the 2d YM theory with an additional Wilson line inserted 



1 c 



(2.20) 



On the reduced manifold one obtains from the Hamiltonian ( ^.17 ) the trigonometric Calogero- 
Sutherland system [p9| with the Hamiltonian 
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E 



sm < y, a > 



For instance, for G = SU{N) one obtains the Hamiltonian with the pairwise interaction 



while for G = SO{2N) one has 



sin (q, - qj) 



(2.21) 



(2.22) 



1 



1 



sin {q^ - qj) 



1 



1 



sin {qi) sin [q^) 



(2.23) 



where VitV2 are the coupling constants. In the general case, one has as many different coupling 
constants as many orbits the Weyl group has in the root system. 

Generalization of the consideration to the generic Wilson line which amounts to the spin trigono- 
metric Calogero system was considered in [|3 . 



2.2 Gauged G/G a-model and Ruijsenaars systems 

In this subsection we show that the Calogero system have a natural relativization which also admits an 
interpretation in field theory terms. The corresponding field theory (given on the phase space before 
the reduction) is nothing but the gauged G/G sigma model which is equivalent to the Chern-Simons 
Lagrangian the central charge of the Kac- Moody algebra being just the deformation parameter, 
from the loop group point of view. 

Note that this relativization arises in different physical problems, for instance, in two dimensional 
systems p6t and in gauge theories in higher dimensions Hence, the unified Lagrangian 

description for these different physical problems is of evident importance. 

As in the non-relativistic system, we perform the Hamiltonian reduction procedure, however now 
the unreduced phase space is the cotangent bundle to the central extended loop group G. As a 
symplectic manifold, it is presented by the following data: (g : S*^ ^ G, c G U{1)\ A E ft^ {S^) ^ g* , k G 
R). 

The group acts as 



g -> hgh ^,A~>hAh ^ + nhd^h ^ 
K — *■ K, c c X S{g, h), 
where S is constructed from the L/(l)-valued 2-cocycle T{g, h) on the group CG {T{g,g~^), 



(2.24) 



Sig,h)^T{h,g)Tihg,h-'). (2.25) 
Then, the group invariant symplectic structure looks as follows 



n= Tv [A{g-'Sgf +SAA g-'Sg] + Kd^gg-\5gg-^ f - K5{d^g)g-Hgg-^] + c'^Sc A Sk 

(2.26) 

The generalization of the Gauss law can be found from the group action 

ti{g, c; A, k) = {gAg-^ + ngd^g-^ - A, 0) (2.27) 

At the next step, we fix the level of the moment map. It is natural to assume that its color 
structure does not change keeping the integrability intact 

fi{g, c; A, n) = iv ( ^Id - e «> e+ ) 5{lp). (2.28) 



As before, we make A a diagonal matrix D hy & generic gauge transformation h{ip) modulo the action 
of the afRne Weyl group. 

The moment map equation now looks 

gDg-^ + ngd^g-^ - D = iv (^Ifd - f ® 5{^) (2.29) 

where / = h{(p)e is the vector from with the unit norm < /, / 1. Using the residual gauge 
symmetry, one can make / G R^, and using (2.2£) one immediately gets 



5 = exp(^I?)G(^)exp(-^Z?) ; d^G ^ --GS{if), (2.30) 

where J = w(-^Id — / <E) /^)- It is useful to introduce the notation for the monodromy of connection 
D: Z = exp(— ^Z?) — diag{zi, . . . , z^), Yli^i — 1' Zi — exp( ^"'^' ). It satisfies the boundary 
condition 

G-^ZG = exp (^^^ Z (2.31) 

where G = G(+0). 

Solving this equation is a separate interesting issue, see |6^. The result looks as follows. First, 
introduce the characteristic polynomial P{z) of the matrix Z 

Piz)^'[[iz-z,). (2.32) 

i 

Define also 

^ F(A±M^F(z^ .2 33) 

^ ^'>- -i)zP'{zy 

where A = e^^. At A ^ 1, the rational functions Q^(z) tend to In these notations, the matrix 
G can be put in the form 

G., = -A-^-^^ e^P'{Q+{z,)Q-{z,)y^' = 

A Z^i — Zj 

^in(^)^ sin(^) sin(^) ^^''^^ 



e 



J--'- sin(^) sin(^) 
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where pi are the quantities that can not be determined from the moment map equation and, similarly 
to the Calogero case, are identified with particle momenta, while qi are particle coordinates, the 
symplectic structure on the reduced manifold being 5pi A Sqi. 

In order to construct Hamiltonians, one may fix any ^c?-invariant function xi9) ^^^d define 



= / dipxig). (2.35) 

Such a Hamiltonian of the dynamical system on T*G is evidently invariant under the group action. 

For instance, one may take ^±(5) = Tr{g±g^^) (where the trace, Tr is taken over the fundamental 
representation of SU{N)) obtaining 

H±=J2{e^''^±e-^^^)l[f{q,,) (2.36) 



where the function /(g) is 

siin^ {nh' / kN) 
sin^(7rq/K) 

One may recognize in this Hamiltonian the Hamiltonian of the trigonometric Ruijsenaars-Schneider 
model |6^, while the solution for g{ip) constructed above gives the Ruijsenaars Lax matrix. Re- 
defining qi nqi and -di — > di/ n, which does not change the symplectic form, and taking the limit 
K — s- 00, we return back to the Calogero-Sutherland system. 

Now we come to the question what is the field theory our system corresponds to. This time we 
choose the Hamiltonian to be zero so that the term / pdq serves as a free action and we also add to 
the action the constraint, the moment map with the Lagrangian multiplier 

S{A, g)^ I dipdt Tr [A^g-^dtg + ^dtg ■ g-' ■ d^g ■ g'' 

^ n (2-38) 

-nd-^d^g ■ g-^dg ■ g-^f) - -A^ {ng-^d^g + g-^A^g - Ai ~ J6{ip)) 

where we have omitted the term Kc^^dtC. This action is the action of the G/G gauged cr-model with 
the Wilson line inserted. 

Now we present the equivalent representation of this theory in terms of the Chern-Simons theory 
with some Hamiltonian. The Chern-Simons theory is the gauge theory given on a three manifold 
which we choose to be the product of an interval and a torus X = I x T^, with the action 



IK 



[ d^xTT{A/\dA+-A/\A/\A). (2.39) 
Jx 3 

We are going to consider the theory with an inserted Wilson line. The phase space of this theory is 
the moduli space of the flat connections on the torus with a marked point with prescribed conjugacy 
class of the monodromy of connection around the marked point. The path integral is of the form 

< wi|Tr(Fexp J A)\v2 > exp{-Scs{A)). (2.40) 

The monodromy class U around the marked point is fixed by the highest weight t) of the represen- 
tation Ri,, U — exp(-j^p^diag(f)j)). The monodromy condition on the torus S = that relates the 
monodromies along the cycles, gA,gB and that around the marked point gc is 



9AgB9A 9b =9c (2.41) 

With the condition on the conjugacy class of gc and taking into account that for the representation 
the signature is hke v-diag^N, 0, 0), one can get equation ( 2.31| ) (see |^ for detailed explanations). 



Let us also point out the equivalence of the sector of the Chern-Simons theory on S x S*^ with 
observables depending on Aq, dt being tangent to S*^, and of the gauged G/G theory. This can be 
proved, say, by gauge fixing constant and diagonal. 

The equivalence provides one with a tool for computation of the spectrum of the Ruijsenaars- 
Schneider model which is contained in the spectrum of the relativistic free particles on the circle. 
In the pure Chern-Simons theory, the particles obey the Fermi statistics, since the phase space is 
^rpAT-i ^ T^^^)/5jv, the moduli space of the flat connections on the torus. Here T^^^ is the Cartan 
torus in SU{N), and Sn is the symmetric group. Since the wave fmiction vanishes at the diagonal in 
the N — 1-dimensional torus, these particles really behave like fermions. 

The Shrodinger equation for the relativistic particle is the difference equation, since cos( ^^]^ dq) 
is a finite-difference operator. Its eigenfunctions on the circle are e^'^*"* with the eigenvalues 

En = cos[^^) (2.42) 
Then, the complete spectrum of the model reads 

H -^"^ (2.43) 

i 

with additional conditions that follow from the invariance En = En+n+N and the symmetry group 
action 

(/t -h iV) > Tijv > . . . > Tij > . . . > Til > (2.44) 
^n,<{K + N) (2.45) 

i 

2.3 Elliptic Calogero model and gauge theories 

There is also another generalization of the trigonometric Calogero model, the elliptic model. In next 
sections we shall see that just the elliptic systems are responsible for nonperturbative dynamics in 
the vacuum sector of the softly broken J\f — ^ SUSY gauge theories, i.e. of the N = 2 SUSY gauge 
theories with the adjoint matter hypermultiplet. 

In this subsection we show how the elliptic models are described via Hamiltonian reduction of the 
holomorphic theories in three dimensions Let us consider a Hamiltonian system whose phase 

space is the cotangent bundle to the algebra of SL{N, (C)-valued functions on the elliptic curve with 
modular parameter r, T*g^^: 

i(l),c;Kd + A), (jj-.Er ^ si NiC), AG r2(°'i)(S^) sZjv(C), c, k G C 

Let us also denote through du; the holomorphic differential with integrals over the A- and B- cycles 
equal to 1 and t, respectively. 

On this phase space acts naturally the current group SLn{C)^^ 



Hz,z) ^ g{z,z)(l){z,z)g{z,z) ^ 
A(z,z)dz — > g{z,z)A{z,z)g{z,z)^^ + Kg{z,z)dg{z,'z)^^ 

K ^ c ^ c + duj /\ Ti{(l)gdg) 



(2.46) 



preserving the holomorphic symplectic form on r*g^^ 

n = 6cASK+ duj Tr(5(^ A SA) (2.47) 
Then, the moment map has the form 

= Kd(j) + \A, (j)] (2.48) 

In order to get the system of interacting particles, one is to enlarge the phase space by adding a finite 

dimensional coadjoint orbit of S'L(Af,C)- A generic complex orbit S'L(iV, (C)/ (^S'L(^ (D) x (C*)') 
is given by a matrix of rank / which can be presented as 

; 

/^.=E<«i (2.49) 

a=l 

Then, the moment map M on the enlarged manifold X^, has the zero level and reads as 

+ \A, 0] - J,j5{z,z)dz A dzuj. (2.50) 

where = {5ij — i)fij- Using gauge transformations, one can obtain for the matrix elements of </> 

Kd(f>tj + qij4>tj = Jij6{z,z), (2-51) 

where = Ui — aj , and Oi € C are diagonal elements of A. From this equation one obtains three 
consequences: 

• Jii = 

• 0ii = Pi = const, i = 1, . . . , N; J2iPi ^ 

• (pij = -Jij exp ( — '- \ 5^ — f^hjX^ foi' * 7^ jj where T2 =Imr, 6'* is the 6'-function with odd 

characteristics]^. 

The first two statements merely follow from the diagonal components of equation ( ^.5l| ), with pi 
being just arbitrary constants of integration (restricted by X^iP* ~ since we consider the SL{N, (C) 
group). To prove the third statement, let us make the substitution 

(t>ij{z,z) = exp ^!^!zlf_^^ (z, z), (2.52) 
so that ipij has the nontrivial monodromy 

^•^We also use the notation Oi for this 0-function. 



on 



4'ijiz + 1,Z + 1) = Ipzjiz, Z) 

_ _ ^ _ (2-53) 

ipijiz + T,z + t) ^ e « ""-'tpijlz^z) 



According to ( 2.51 ), we have d^^ij = Jij5{z, z), i.e. ij^ij is a meromorphic section of the holomorphic 
linear bundle with the pole at the origin. 

The matrix cj) obtained is nothing but the Lax operator for the elliptic spin Calogero system jTl], Q . 
In order to get spinless Calogero case, one needs to choose the simplest orbit I = 1 (which describes 
the fundamental representation). In this case, the factor matrix Jij in (j) can be trivially gauged away 
to an arbitrary scalar —nv. One should also introduce the new rescaled variables qi = — 7rai/KT2 
identified with coordinates in the many-body system, pi being the corresponding momenta, and one 
obtains finally the Lax operator for the elliptic Calogero system which being calculated at z ^ 0, 
coincides with the Lax matrix found by Krichever earlier fz^ . 

Now invariants of the algebra yield Hamiltonians of the integrable system. The simplest one 
coincides, up to a constant, with the second Hamiltonian of the elliptic Calogero-Moser system 



H= d'z-TTcbiz,zy ^ d'zl^-pi + u 



J<3 

2 j 



(2.54) 



Let us remark that there is an essential difference between elliptic model and its degeneration at the 
quantum level. It appears that it is impossible to get the quantum spectrum of the elliptic model 
analytically. The only tool available is the Bethe anzatz approach which yields the equation for the 
eigenvalues Jz^ . Some quantum aspects of these systems within the representation theory framework 
were discussed in [ [zs] . 

One may want to obtain an elliptic deformation of the two-dimensional YM theory. The action 
proposed in [TO is of the form 



Sr= / ujATT{cj)Ft-~eq^') (2.55) 

where 5*^ is a time-like circle. One should certainly also insert a (vertical) Wilson line Pexp (i^ / At) 
in representation R^. 

Let us emphasize that the elliptic Calogero-Moser model is an example of the Hitchin systems given 
on the moduli space of the holomorphic vector bundles |7^, ^ ^ for the surfaces with marked 
points. The elliptic Calogero-Moser system has a natural relativization - the elliptic Ruijsenaars- 
Schneider model which is generalization of the Hitchin systems. A description of the elliptic Ruijsenaars- 
Schneider system via Hamiltonian and Poissonian reductions can be found in . 

Recently the Hamiltonian reduction procedure for the elliptic Calogero models considered above 
for An case has been extended to other classical groups in |5l| . 



3 Particle systems and supersymmetric gauge theories 

After having presented a manifest construction of the integrable systems, we are ready to show how 
these integrable many-body systems emerge in the consideration of effective low-energy M = 2 su- 



persymmetric gauge field theories. As it will be clear later, the role of the many-body systems is to 
provide the proper degrees of freedom in the vacuum sector of the theory. The fact that the complete 
set of data fixing the low-energy sector of the theory is governed by the many-body system definitely 
implies that the number of the effective degrees of freedom in the vacuum sector of the supersymmetric 
gauge theory is finite. In this section, we demonstrate how the corresponding particle systems can be 
recovered, while later we comment on the role of these degrees of freedom within the string theory 
context. 

3.1 Seiberg-Witten theory and integrable systems 

As explained in the Introduction, the set of data a Riemann surface, the corresponding moduli space 
and the differential dS that gives the physical quantities, in particular, prepotential is naturally 
associated with an integrable system. In the case of the SU{2) pure gauge theory, which has been 
worked out in detail in [|l| this correspondence can be explicitly checked [^7| . In particular, one may 
manifestly check formulas (^^), since all necessary ingredients are obtained in In this way, one 
can actually prove that this Seiberg-Witten theory corresponds to the integrable system (periodic 
Toda chain with two particles, see below). However, to repeat the whole Seiberg-Witten procedure 
for theories with more vacuum moduli is technically quite tedious if possible at all. Meanwhile, only 
in this way one could unquestionably prove the correspondence in other cases. 

In real situation, the identification between the Seiberg-Witten theory and the integrable system 
comes via comparing the three characteristics 

• number of the vacuum moduli and external parameters 

• perturbative prepotentials 

• deformations of the two theories 

The number of vacuum moduli (i.e. the number of scalar fields that may have non-zero v.e.v.'s) on 
the physical side should be compared with dimension of the moduli space of the spectral curves in 
integrable systems, while the external parameters in gauge theories (bare coupling constants, hyper- 
multiplet masses) should be also some external parameters in integrable models (coupling constants, 
values of Casimir functions for spin chains etc). 

As for second item, we already pointed out the distinguished role of the prepotential in Seiberg- 
Witten theory (which celebrates a lot of essential properties which we shall discuss later). In the 
prepotential, the contributions of particles and solitons/monopoles (dyons) sharing the same mass 
scale, are still distinguishable, because of different dependencies on the bare coupling constant, i.e. on 
the modulus t of the bare coordinate elliptic curve (in the UV-finite case) or on the Aqcd parameter 
(emerging after dimensional transmutation in UV-infinite cases). In the limit r —f ioo {Aqcd 0), 
the solitons/monopoles do not contribute and the prepotential reduces to the "perturbative" one, 
describing the spectrum of non-interacting particles. It is immediately given by the SUSY Coleman- 
Weinberg formula : 

^pert(a)- J2 (-)^TrK(a + Af^)'log(a + MO ^3_^) 

reps R,i 

Seiberg-Witten theory (actually, the identification of appropriate integrable system) can be used to 
construct the non-perturbative prepotential, describing the mass spectrum of all the "light" (non- 



stringy) excitations (including solitons/monopoles). Switching on Whitham times presumably 
allows one to extract some correlation functions in the "light" sector (see below). 

The perturbative prepotential (3J) has to be, certainly, added with the classical part. 



^ciass = -ra^ (3.2) 

since the bare action is S{'i!) = Inir / Tr^''^. r here is value of the coupling constant in the deeply UV 
region. It plays the role of modular parameter of the bare torus in the UV finite theories. On the other 
hand, in asymptotically free theories, where the perturbative correction behaves like a^log{a/ Aqcd), 
T is just contained in Aqcd, since can not be invariantly defined neither in intcgrable nor in gauge 
theory. 

As we already mentioned the problem of calculation of the prepotential in physical theory is simple 
only at the perturbative level, where it is given just by the leading contribution, since the /3-function 
in A/" = 2 theories is non-trivial only in one loop. However, the calculation of all higher (instantonic) 
corrections in the gauge theory can be hardly done at the momentp|. Therefore, the standard way of 
doing is to make an identification of the Seiberg-Wittcn theory and an integrable system and then to 
rely on integrable calculations. This is why establishing the "gauge theories ^ integrable theories" 
correspondence is of clear practical (apart from theoretical) importance. 

Now let us come to the third item in the above correspondence, namely, how one can extend 
the original Seiberg-Witten theory. There are basically three different ways to extend the original 
Seiberg-Witten theory . 

First of all, one may consider other gauge groups, from other simple classical groups to those 
being a product of several simple factors. The other possibility is to add some matter hypermultiplets 
in different representations. The two main cases here are the matter in fundamental or adjoint 
representations. At last, the third possible direction to deform Seiberg-Witten theory is to consider 5- 
or 6-dimensional theories, compactified respectively onto the circle of radius i?5 or torus with modulus 
Rb/Re (if the number of dimensions exceeds 6, the gravity becomes obligatory coupled to the gauge 
theory). 

Now we present a table of known relations between gauge theories and integrable systems, Fig.l. 
Before coming into details, let us briefly describe different cells of the table. 

The original Seiberg-Witten model, which is the 4d pure gauge SU{N) theory (in fact, in their 
papers ||, |^, the authors considered the SU{2) case only, but the generahzation made in [|6) is quite 
immediate), is the upper left square of the table. The remaining part of the table contains possible 
deformations. Here only two of the three possible ways to deform the original Seiberg-Witten model 
are shown. Otherwise, the table would be three-dimensional. In fact, the third direction related to 
changing the gauge group, although being of an interest is slightly out of the main line. Therefore, 
we only make several comments on it. 

One direction in the table corresponds to matter hypermultiplets added. The most interesting is 
to add matter in adjoint or fundamental representations, although attempts to add antisymmetric 
and symmetric matter hypermultiplets were also done (see [ p5[ for the construction of the curves and 
for the corresponding integrable systems). Adding matter in other representations in the basic 
SU (N) case leads to non- asymptotically free theories. 

^^In order to check the very ideology that the integrable systems lead to the correct answers, there were calculated first 
several corrections The results proved to exactly coincide with the predictions obtained within the integrable 

approach. 

oo 
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Figure 1: SUSY gauge theories ■(=> integrable systems correspondence. The perturbative limit is 
marked by the italic font (in parenthesis). 

Columns: Matter in adjoint vs. fundamental representations of the gauge group. Matter 
in adjoint representation can be described in terms of a larger pure SYM model, either with higher 
SUSY or in higher dimensional space-time. Thus models with adjoint matter form a hierarchy, nat- 
urally associated with the hierarchy of integrable models Toda chain ^ Calogero ^ Ruijsenaars ^ 
Dell (mm H H, |0[ IIHI, H, H, ||. similarly, the models with fundamental matter [|, || are 
related to the hierarchy of spin chains originated from the Toda chain: Toda chain ^ XXX ^ XXZ 
XYZ gl^gll, ||. 

Note that, while coordinates in integrable systems describing pure gauge theories and those with 
fundamental matter, live on the cylinder (i.e. the dependence on coordinates is trigonometric), the 
coordinates in the Calogero system (adjoint matter added) live on a torus|^. However, when one takes 
the perturbative limit, the coordinate dependence becomes trigonometric. 

Lines: Gauge theories in different dimensions. Integrable systems relevant for the descrip- 
tion of vacua of d = 4 and d — 5 models are respectively the Calogero and Ruijsenaars ones (which 
possess the ordinary Toda chain and "relativistic Toda chain" as Inozemtsev's limits |lOC|1 ), while 
d = 6 theories are described by the double elliptic (Dell) systems. When we go from 4d (Toda, 
XXX , Calogero) theories to (compactified onto circle) 5c? (relativistic Toda, XXZ , Ruijsenaars) 
theories the momentum-dependence of the Hamiltonians becomes trigonometric (the momenta live on 
the compactification circle) instead of rational. Similarly, the (compactified onto torus) 6d theories 

^■^Since these theories are UV finite, they depend on an additional (UV-regularizing) parameter, which is exactly the 
modulus T of the torus. 



give rise to an elliptic momentum-dependence of the Hamiltonians, with momenta living on the com- 
pactification torus. Since adding the adjoint hypermultiplet elliptizes the coordinate dependence, the 
integrable system corresponding to 6c? theory with adjoint^^ matter celebrates both the coordinate- 



and momentum-dependencies elliptic. A candidate for "the elliptic Toda chain" was proposed in [101|. 

Further in this section we show how the Riemann surfaces and the meromorphic differentials arise 
when considering a concrete integrable many-body system and explain in more details different squares 
in the table of correspondence between Seiberg-Witten theories and integrable systems. However, of 
all 6d theories we describe here only the XYZ chain case, postponing the discussion of the most 
general adjoint matter system, Dell system (the right bottom cell of the table) till section 7 where 
this integrable model naturally emerges within the framework of duality. 



3.2 Ad pure gauge theory: Toda chain 

We start with the simplest case of the 4c? pure gauge theory, which was studied in details in the 
original paper for the SU{2) gauge group and straightforwardly generalized to the SU{Nc) case in 

m. 



This system is described by the periodical Toda chain with period Nc ||37 , 102 , 102 ] whose equations 
of motion read 



(3.3) 



dt dt 

with the periodic boundary conditions Pi+N^ — Pi^ Qi+N^ = Qi imposed. In physical system, there 
are Nc moduli and there are no external parameters {Aqqd can be easily removed by rescaling). In 
the Toda system there are exactly conservation laws. These conservation laws can be constructed 
from the Lax operator defined for any integrable system. The eigenvalues of the Lax operator do not 
evolve, thus, any function of the eigenvalues is an integral of motion. 

There are two different Lax representation describing the periodic Toda chain. In the first one, 
the Lax operator is represented by the Nc x Nc matrix depending on dynamical variables 



I 



e2 



Pi 

1(92-91) 



i_p3(9i-9«e) 



W 



,^(92-91) 
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1(93-92) 



The characteristic equation for the Lax matrix 
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-P3 






PNa 



(3.4) 



V{X, w) = detAT^xJV, {^'^^{w) - A) = 
generates the conservation laws and determines the spectral curve 



(3.5) 



1 



+ - = 2P^v,(A) 
w 



(3.6) 



^^Let us point out once more that by adding the adjoint matter we always moan soft breaking of higher 
supersymmetries . 



where Pn^ W is a polynomial of degree Nc whose coefficients are integrals of motion. If one restores 
the dependence on Aqcd in this spectral curve, it takes the form 



^+^^2P^(A) (3-7) 
w 

However, we typically re-scale it away, for exception of s.6.3, where we discuss the manifest RG 
dynamics w.r.t. Aqcd- 

The spectral curve ( |3.6| ) is exactly the Riemann surface introduced in the context of SU{Nc) gauge 
theory. It can be also put into the hyperelliptic form 

2Y = w--, y2^p2(A)-l (3.8) 
w 

The integrals of motion parameterize the moduli space of the complex structures of the hyperelliptic 
surfaces of genus Nc — 1, which is the moduli space of vacua in physical theory. 

If one considers the case of two particles {SU{2) gauge theory), P(A) = — u, u = p'^ — coshg, 
p = Pi — —p2, q ~ <li — 02- Thus we see that the order parameter of the SUSY theory plays the 
role of Hamiltonian in integrable system. In the perturbative regime of the gauge theory, one of the 
exponentials in coshg vanishes, and one obtains the non-periodic Toda chain. In the equation ( [3.6| ) 
the perturbative limit implies vanishing the second term in the l.h.s., i.e. the spectral curve becomes 
the sphere w = 2P]y^{X). 

The alternative Lax representation appeals to introducing the local operators at each site of the 
chain 

( \ -p^ e« 

Lf^(A)= |, z=l,...,iVe (3.9) 

y -e-«' 

This Lax operator "makes a shift" to the neighbour site so that the auxiliary linear problem is 

L,(A)vI/,(A)-vi/^+i(A) (3.10) 

where ^'i(A) is a two-component Baker- Akhiezer function. The periodic boundary conditions are 
easily formulated in terms of this Baker- Akhiezer function and read as 

*,+w^(A) =?«*,(A) (3.11) 

where w is a free parameter (diagonal matrix) . One can also introduce the transfer matrix that shifts 

i to i + Nc 

r(A) = LAr^(A...Li(A) (3.12) 
Then, the boundary conditions imply that r(A)5'i(A) = w'^i{X), i.e. 

det2x2(Tiv,(A)-u.) =0 (3.13) 
providing us with the same spectral curve 



= w^ - wTtTn^{X) + detTN^{X) ^ - wTtTn^{X) + 1, 

1 1 (3.14) 

i.e.: ViX.w) = Trrjv,(A) - w - - ^ 2Pjv,(A) -w--=0 

w w 

The identification of tlie two representations can be also done at tfie level of their linear problems. 



Indeed, let us consider the Lax operator ( |3.9| ) with the linear problem ( 3.1C ) and the boundary 



conditions (3.11). If we parameterize = ( * ) , then the linear problem ( |3.10| ) can be rewritten 

' X^ J 



as 



and, along with the periodic boundary condition ( ^.ll| ) reduces to the linear problem C{w)^ = A<I> 



for the Nc x Nc Lax operator (3^) with the iVc-component Baker-Akhiezer function $ = {e '^^^'^tpi}- 

After having constructed the Riemann surface and the moduli space describing the Ad pure gauge 
theory, we turn to the third crucial ingredient of Seiberg-Witten theory that comes from integrable 
systems, the generating differential dS. The general construction was explained in the Introduction, 
this differential is in essence the "shorten" action pdq. Indeed, in order to construct action variables, 
tti one needs to integrate the differential dS = '^^Pidqi over Nc — 1 non-contractable cycles in the 
Liouville torus which is nothing but the level submanifold of the phase space, i.e. the submanifold 
defined by values of all Nc~ 1 integrals of motion fixed. On this submanifold, the momenta pi are 
functions of the coordinates, qi. The Liouville torus in Seiberg-Witten theory is just the Jacobian 
corresponding to the spectral curve (|l.6| ) (or its factor over a finite subgroup) . 

In the general case of a ^-parameter family of complex curves (Riemann surfaces) of genus g, the 
Seiberg-Witten differential ^5* is characterized by the property 5dS = X]i=i Suidvi, where dvi{z) are 
g holomorphic 1-differentials on the curves (on the fibers), while 5ui are variations of g moduli (along 
the base). In the associated integrable system, Ui are integrals of motion and tt^, some g points on 
the curve are momenta. The symplectic structure is 

g g 
^ da^ A dpi"" = ^ dw, A dv^{■Kk) (3.16) 

i—l i,k—l 

The vector of the angle variables, 

pf- = ^ / du;, (3.17) 

k=l •' 

is a point of the Jacobian, and the Jacobi map identifies this with the g-th power of the curve, 
Jac = C®9 , Here doji are canonical holomorphic differentials, dvi = X]f=i ^"-^i §a ^'^'^'^ details 
on the symplectic form on the finite-gap solutions can be found in |l04| . 

Technical calculation is, however, quite tedious. It is simple only in the 2-particle {SU{2)) case, 
when the Jacobian coincides with the curve itself. In this case, the spectral curve is 

w + -^2(X'^-u) (3.18) 
w 

while u — p^ — coshq. Therefore, one can write for the action variable 



Of? 



dw ,„ , dw 



a ~ f pdq = f \Ju — cosh qdq = f A — , dS — X (3.19) 



where we made the change of variable w = e'^ and used equation ( 3.1§| ) 



Now one can naturally assume that this expression for the differential dS is suitable for generic 
Nc- A long calculation (which can be borrowed, say, from the book by M.Toda |105| ) shows that this 
is really the case. One can easily check that the derivatives of dS w.r.t. to moduli arc holomorphic, 
up to total derivatives. Say, if one parameterizes^ Pn^{X) = — A^^ + SAr^_2A^'=^^ + ... and note that 
dS — Xdw/w = XdPM^{X) /y = PN^{X)dX/Y + total derivatives , then 

ddS^X^ (3.20) 
dsk Y 

and these differentials are holomorphic if fc < Nc — 2. Thus, iVc — 1 moduli gives rise to A^c — 1 
holomorphic differentials which perfectly fits the genus of the curve (we use here that there is no 
modulus sat^-i). 

It turns out that the form of dS = Xdw/w is quite general and does not change even in more 
complicated cases of spin chains, which we discuss in the next subsection. 

To conclude this subsection, let us briefly consider the theories with other gauge groups. First of 



all, we present an invariant algebraic formulation of the Toda Lax operator (3.4). It can be obtained 



much along the line discussed in the previous section, i.e. via the Hamiltonian reduction |6(|, or even 



obtained from the elliptic Calogero model by a degeneration [100| (we discuss it later). Here we just 
write down the answer which looks like [|lO|, |§ 



Na-l 



Ciw) = }_J [e"'-i (e„^ + e„J + p^h,] + u;e„„e"«-i + le_,„e"«"i (3.21) 



w 



where the sum runs over all simple roots ai of the affine algebra A^''^-^, Gq,, is the highest (long) root, 
hi span the Cartan subalgebra and the vector q has as its components the particle coordinates. Note 



that the generators in (3.4) are taken in the simplest fundamental representation of the algebra. 
However, this requirement is not essential for the whole construction and can be removed. 

This expression is straightforwardly extended to any simple affine Lie algebra^ Nc — 1 being 
substituted by the rank of the group. Now, in order to construct Seiberg-Witten theory with a simple 



gauge group G one needs to consider the Toda chain for the corresponding dual affine algebra [102|. 
Say, the gauge theory with the group Sp{2n) is described by the affine Toda for (al'^) etc. 

Having the Lax operator, one may repeat the procedure of this subsection in order to get the 
spectral curve and the corresponding moduli space. As for the generating differential dS, it preserves 
the form dS = Xdw/w for all groups. 

Note that for any, at least, classical simple group there is also a 2 x 2 Lax representation generalizing 
(p^). In this case, the Lax operator is given on the Dynkin diagram as a substitute of the closed 
chain above and the definition of transfer matrix is to be extended to include the reflection matrices 
[ |l07| . Therefore, one can immediately construct some Seiberg-Witten theory in this way. We refer 



for the details to [108| 



^^The absence of the term A^"^"^ is due to the SU{Nc) group and corresponds to the total momentum equal to zero, 
i.e. to the center mass frame. This is important for further counting of holomorphic differentials. 

Sometimes one needs to add non-unit coefficients, depending on the root, in front of exponentials. 



oo 



3.3 Adjoint hypermultiplet added: Calogero system 



As far as there exist two different Lax representations, there are, at least, two possible generalizations 
of the periodic Toda chain (in fact, more as we shall see later). The first generalization is related 
to the Nc X Nc representation (3.4). On physical side, it corresponds to the simplest deformation of 



the Af — 2 SUSY gauge theory by adding the adjoint hypermultiplet |87|, |90|, |89|. If the matter 
is massless, the theory is finite and the coupling constant can be defined at arbitrary scale. We 
mainly discuss the massive matter and find the corresponding dynamical system with the additional 
parameter, ft appears that the deformation by the massive hypermultiplet can be described in terms 
of the elliptic Calogero model. The hypermultiplet mass plays the role of the coupling constant in the 
Calogero model, while the UV value of the field theory coupling constant defines the modulus of the 
torus where the Lax operator is defined. When one sends the mass of the hypermultiplet to infinity, 
in field theory the dimensional transmutation procedure is applied, which has its counterpart in the 
many-body system and is equivalent to the transition back from Calogero to the Toda system. 

The matrix Nc x Nc Lax operator for the Calogero system was found in section 2 and read as 



\ i aeA J 



( Pi F{qi-q2\0 

F{q2-qi\£.) P2 

V FiqN^-q,\0 FiqN^~q2\0 



F{qi-qNM) \ 
F{q2 ~ qNM) 



PN^ J 
(3.22) 

where the first sum runs over the Cartan subalgebra of SU{Nc), while the second one runs over all 
roots (not only over the simple ones, in variance with the Toda case). The special functions that 
appear above are defined [O 



Fm 



<^{q - 



(3.23) 



where (t(^) denotes the (Weierstrass) cr- function |l09[] . One identity used throughout the calculations 
with the Calogero (and Ruijsenaars below) system is 



a{u — v)a{u + v) 
where the Weierstrass p-function is defined as^ 



= p{v) ~ p{u). 



1 



sinh (^ -|- tot) 



(3.24) 



(3.25) 



This Lax operator (3.22) can be just read off from the Calogero Lax operator obtained in s.2.3 (after 
a trivial gauge thransformation). The only external parameter of the integrable system in this case, 
the coupling constant v is proportional to the only external parameter of the physical theory, the 
hypermultiplet mass M^. 



^*In accordance with the standard definition, 109 1 p{^) is the doubly periodic Weierstrass function with periods 
2uj and 2lj' , and t = We remark, however, that although the Weierstrass function depends on two periods, the 
homogeneity relation 

p{tz\tLLl,tul') = t~'^ p{z\uj,0j') 

enable us to arbitrary scale one of these. All our results here are ultimately independent of such scaling and this allows 
us to choose the real period to be tt (that is a; = 7r/2). We make this choice throughout. 



The spectral curve S*^"' 

deW^xiv. (/:'^"'(0~A) =0 (3.26) 
covers the "bare" spectral curve, torus E{t) 

= {x - ei){x ~ e2){x - 63) (3.27) 
manifestly parametrized by the variable x — p{£,), y = — |p'(^). 

Y^Cai ^^j^ |-|g presented as a polynomial V{x,y;X) = of degree Nc in A, the manifest (quite 
involved) expression for it can be found in ^ |llC[] . However, this manifest expression can be 
always obtained by a degeneration from the Ruijsenaars spectral curve discussed later. Amazingly, 
this latter curve can be written in a quite compact form. 

In the simplest case of two particles, the curve S^?"' is 

V{x,y;\) = -u + iy^p{£,) ^0 (3.28) 
The generating differential dS on the Calogero curves has the form 

dS = Xdcu, duj = d^ = -2— (3.29) 

y 

Here dio is the holomorphic 1-differential on the bare spectral torus and one can check that the 
derivatives w.r.t. moduli of this differential are holomorphic on Y^^'^K 

Now let us make the dimensional transmutation, i.e. reproduce the periodic Toda chain from the 
Calogero system [100, To this end, first, we degenerate the bare spectral torus r ioo. 



Then, to provide the nearest neighbor interaction, we introduce the homogeneous coordinate "lattice" 
with the large distance A between sites 

qj - J A + 0, (3.30) 

To see explicitly what kind of interaction emerges in this limit, it suffices to look at the Weierstrass 
function giving the Calogero potential, although the procedure can be easily repeated for the Lax 



operator. The Calogero potential has form (see (2.54)) 



Vix.j) ^ ly^Y^Pilij) (3.31) 

In the limit under consideration we introduce the renormalized coupling constant 1^ = 1^0 exp(A) so 
that nuo will be ultimately nothing but Aqcd- Now, choosing A ~ t and taking the limit r — > ioo, 
we see that only m = term survives in the sum ( 3.25 ) in the potential ( 3.3l| ) so that the resulting 



potential reads as 

Na-l 

Voicl^^) = ^ e*'+^-^- (3.32) 

and describes the open (non-periodic) Toda chain. 

In order to get the periodic Toda chain, one needs to fix A = — . In this case, the m = — 1 term in 



2q) also contributes into the sum and one finally obtains the potential 



VTcm^Vo + vle't''-'t'- (3.33) 

describing the periodic Toda chain. 



Now, the bare torus (3.27) under the procedure described degenerates into a sphere with two 



punctures, at w = and w — oo: y ^ w — x w -\- l/w, while the generating differential dS 



turns into Xdw/w [duj — > dw/w) coinciding with (3.19) 



Returning to physical theory, making the limit to the pure gauge theory one needs to send the 
hypermultiplet mass m to infinity keeping finite A^'' ~ M^^e^"^" . This is exactly the limit we just 
discussed in integrable system. 

To conclude our consideration, let us note that the Calogcro Lax operator considered above can 



be generalized to other groups [111|, and the corresponding Seiberg-Witten theory can be constructed 



[ |112| . However, since the exact construction requires entering many technicalities we skip it here. 
3.4 Fundamental matter added (SQCD): XXX spin chain 

The second possible generalization is related to the 2x2 representation. This generalization on physical 
side corresponds considering theories with massive hypermultiplets in the fundamental representation. 
To preserve the asymptotic freedom one has to demand for the number of matter hypermultiplets, 
Nf not exceed 2Nc. In particular, the theory with zero /^-function, i.e. UV finite one corresponds 
to Nf = 2Nc. On integrable side, this theory is described by the inhomogeneous twisted XXX spin 
chain, while theories with Nf < 2Nc are described by its degenerations jo^, |38| . 

The periodic inhomogeneous sl{2) XXX chain of length Nc is given by the 2x2 Lax matrices 



L,{X) = i\ + X^)■l + Y,Sa.^■<J'' (3.34) 

cr" being the standard Pauli matrices and Xi being the chain inhomogeneities, and periodic boundary 
conditions are imposed. The linear problem in the spin chain has the same form ( pO| ) as in the Toda 
case. 



One can also introduce the transfer matrix as in ( ^.12 ) which provides the spectral curve equation 



{ \i.l2 ) and generates a complete set of integrals of motion. 



Integrability of the spin chain follows from the quadratic r-matrix relations (see, e.g. ]113| ) 

{L,(A) f Lj(A')} = <5y [r(A - A'), L,{X) (g, L,(A')] (3.35) 
with the rational r-matrix 

1 ^ 

r(A) = -^cT^^fT" (3.36) 

a— 1 



The crucial property of this relation is that it is multiplicative and any product like ( 3.12| ) satisfies 
the same relation 

{T(A) ? r(A')} = [r{X - A'), T(A) (g, T(A')] (3.37) 
The Poisson brackets of the dynamical variables Sa, a — 1,2,3 (taking values in the algebra of 



functions) are fixed by (3.35) and arc just 



{Sa; Sb} = -if-abcSc (3.38) 

i.e. the vector {Sa} plays the role of angular momentum ("classical spin") giving the name "spin- 
chains" to the whole class of systems. Algebra ( 3.38| ) has an obvious Casimir function (an invariant, 
which Poisson commutes with all the spins Sa)^ 

3 

i^' = S2 = ^ SaSa (3.39) 



The spectral curve ( 3.13 ) is explicitly now 



^_^ Qw,(A) ^2Pa,^(A), 2Pjv.(A) =Trr(A), Qjv, (A) = detr(A) (3.40) 
w 

or in the hyperelliptic parameterization 



^P],^{\)~Qn,{\) (3.41) 
Zeroes of QNf{^) define the masses of the hypermultiplets. Since 

det2x2i.(A) = (A + A,)' - (3.42) 

one gets 



QNfW = det2x2T(A) = I]det2x2i»(A) = H ((^ + ^»)' " = 



i=l 



(3.43) 



= [](A-m+)(A-m-) 

1=1 

where we assumed that the values of spin K can be different at different sites of the chain, and 



m 



f^~X,±K,. (3.44) 



Thus we obtain that Nf is generally equal to 2Nc. Note that the hypermultiplet masses being external 
parameters in the gauge theory are also such in the spin chain. Indeed, there are still Nc — 1 integrals 
of motion that parameterize the spectral curve and the moduli space (i.e. parameterize the Coloumb 
branch of gauge theory), but the system also depends on 2Nc additional external parameters, Nc 
Casimir functions and Nc inhomogeneities. 



While the determinant of the monodromy matrix (3.42) depends on dynamical variables only 
through the Casimirs Ki of the Poisson algebra, the dependence of the trace Tr2x27^(A) is less trivial. 
Still, it depends on Sa^ only through Hamiltonians of the spin chain (which are not Casimirs but 
Poisson-commute with each other) - see further details in . 

Let us note that we have some additional freedom in the definition of the spin chain and the spectral 



curve. Namely, note that r-matrix (3.36) is proportional to the permutation operator P(X Y) = 
Y ^ X. Therefore, it commutes with any matrix of the form U ® U . Thus, one can multiply Lax 
operator of the spin chain by arbitrary constant matrix without changing the commutation relations 
and conservation laws. Moreover, one can also attach a constant (external magnetic field) matrix V 



to the end of the chain (to the A'c-th site). This is the same as to consider more general boundary 
conditions - those with arbitrary matrix V~^. This is why such a model is called twisted. 

The described freedom allows one to fit easily the form of the spectral curve proposed in [|[ |ll4| 



w 

P{X)^l[{\-cl>,), Q{X) = h{h + l)l[ix 

3 = 1 \ 



i=l 



h{r) 




where r is the bare curve modular parameter and 9i are the theta-constants. 
It can be done, e.g., by choosing the matrices U and V to bc[^ 



(3.45) 



1 



a,; 



V 



I.e. 



1 



detV \\oii 




(3.46) 



Now let us briefly consider degenerating our UV finite system to an Nf < 2Nc case. This can 
be done in the standard way j^] by sending I masses mi,..., mi to infinity while keeping Ag^^j = 



e^'^'^mi...mi finite. After this procedure the modular forms disappear from (3.45) so that Aqcd 
emerges instead. 

Degenerations of the system can be studied at a single site (for the sake of brevity, we omit the 
index of the site). Let us consider the Lax operator L = UL (see (3.34), ( 3.46| )) with spins satisfying 



the Poisson brackets (3.38). We are going to send a to zero. It still reserves two possibilities to get 



nontrivial Lax operator. The first possibility, when the both masses ( 3.44 ) disappear and one reaches 
the pure gauge theory, is described by the periodic Toda chain. In order to get it, one needs to redefine 
5*+ — > ^5+, then to send a to zero and to remove, after this, the inhomogeneity by the shift of Sq. 
This brings us to the Lax operator of the form (we introduce new notations Sq — S3, S± — Si ± 182) 




so that the Poisson brackets are 



(3.47) 



{S±,Sq}^±S±, {5+,5_} = 
This algebra is realized in new (Heisenberg) variables p and q 



(3.48) 



(3.49) 



This leads us finally to the Toda chain Lax operator ( |3.9| ) and the Toda spectral curve (3.14) 
i^To fit ( ^.45| ), we also need to shift Ai ^ - ^ ^ . m;. 



Now let us return to the second possibility of the asymmetric degeneration, when one of the masses 
( ^.44| ) remains in the spectrum while the second one goes to infinity. One can understand from ( 3.44 ) 
that, in contrast to the Toda case, this degeneration requires a special fine tuning of the Casimir 
function and inhomogeneity, so that both of them go to infinity but their sum (difference) remains 
finite. In the Lax operator, it can be done in the following way. Let us redefine 5+ — s- ^S-^- and 
So —Sq. This means that the Poisson brackets take the form 



{5+,5_} = -25o, {^±,5o}-0 



(3.50) 



Now in order to preserve the finite Lax operator ( 3.34 ), one needs to take care of its element iii(A). 
This can be done by rescaling Ai ^Xi and fixing Ai + 5o to be a • sq. This brings us to the Lax 
operator 



X + so \ / A + so 

S+ Xi — So J y 5*+ — 25*0 
The determinant of this Lax operator is equal to (A — m) where m is the finite mass 



(3.51) 



TO — 2soSo + S+S- (3.52) 
in perfect agreement with ( t3.44| ). Let us note that 2so<5'o + S+S- is also the Casimir function of the 



algebra (3.5C), since {S'±,so} = iS±. 



Let us note that another (equivalent) way to count all possible degenerations |115| is to consider 
2x2 Lax operator of the most general form linear in A, which satisfies the Poisson brackets ( 3.35| ) 



with the rational r-matrix (3.36) and to determine all Casimir functions with respect to this Poisson 
brackets. Then, all possible degenerations are determined by vanishing the Casimirs functions. 

Let us describe another possible interpretation of formula ( 3.44 ). Consider a free particle in d = 4 
with the dispersion law 



and note that in the complex momentum space there is the level crossing at p = in the massless case 
and on the surface = —rr? in the massive situation. According to general logic of the Berry phase 
phenomenon, if one considers the momentum space as a parameter space each level crossing point is 
associated with monopoles and each level crossing surface with "defects" of higher dimensions. 

Therefore, the masses define the locations of singularities in the momentum space, indeed. In 
section 8, we shall see that at the points of the space of scalar field zero modes with coordinates equal 
to masses, there are branes of different dimensions located. This fits the above arguments, since, 
due to the gauge invariance, momenta and gauge fields enter the equations of motion combined in 
covariant derivatives. 

The monopole singularity at p = in the massless case results in an additional term in the 
rotation operator in the momentum space, the phenomenon well-known for the rotation operator 
in the monopole background. Moreover, the monopole contribution to the rotation operator allows 
one to define a chirality transformation for the spin 0, 1 ,2 particles in free theory as the rotation 
when carrying over closed contour in the momentum space. In the interacting theory, the chirality 
transformation is generically anomalous and, therefore, such a simple reasoning fails. However, the 
anomaly itself can be described in terms of the Berry phase framework in the space of fields both for 

O A 



the external and internal anomalies. Then, one has again a satisfactory qualitative agreement, since 
the chiral invariance is broken either by mass or by anomaly: in both cases the violation is related with 
the appearance of singularities in the "generalized momentum space" . An analogy with the Peierls 



model discussed in this context in [116 provides an additional support for this point of view. 

To conclude this long subsection, let us mention two other deformations of the original gauge 
theory. First of all, as in all other cases, one may consider the theory with a different gauge group 
G. This case is treated similarly to the Toda case by introducing reflection matrices and considering 



the spin chain on the Dynkin diagram of g V |108|. Then, the spectral curves can be described by the 
general formula 

P{X,w) = 2P{X)-w-^^ (3.53) 

w 

Here P{X) is the characteristic polynomial of the group G for all G ^ Cn, i.e. 

P(A) = det(.g - XI) = 1[{X - A,) (3 54) 

i 

where determinant is taken in the first fundamental representation and 's are the eigenvalues of the 
algebraic element g. For the pure gauge theories with the classical groups ]102( |, Q{X) ~ A^'* andp^ 



A„„i: P(A) = [](A - A,), s = 0; 

i=l 
n 

i?„: P{X) = Xl[{X^ - X^,), s = 2; 



„ (3.55) 
C„: P(A) = n(A2-A,2)- A, s = ~2; 

i—1 

n 

Z?„: P{X) = l[{X'-X^), s^2 

1=1 

For exceptional groups, the curves arising from the characteristic polynomials of the dual affine alge- 
bras do not acquire the hyperelliptic form. 

In order to include Np massive hypermultiplets in the first fundamental representation one can 
just change A^^ for Q(A) = \^'Y\^J:^{\ - rnj if G = A„ and for Q(A) = A^'* H^IK-^^ - if 



G = B„,C„,i?„ [114, 117, 108[ 



Another possibility of deformation, letting the gauge groups be products of simple factors and 
bi-fundamental matter be added, was proposed in [5l| ]. This system is described by a higher, sl{p) 
XXX magnet [ ^ , with the Lax operator given by the p x p matrix at each site 

L,{X)=K''''Sa,^Xa + iX + X,)■l (3.56) 

where Sa.i are dynamical variables, Xa are generators of the sl{p) algebra and K"-'' is its Killing form. 

Thus, we described the families of spectral curves for theories with fundamental matter. Note 
that, at any XXX spin chain, the generating differential has the same form as in the Toda theory 



(p.l£), dS = Xdw/w. 



-*In the symplectic case, the curve can be easily recast in the form with polynomial P(A) and s = 0. 



3.5 5d SQCD: twisted XXZ chain 



Now we describe the integrable system behind the 5d theory with fundamental matter hypermultiplets. 
The theory is considered with one dimension compactified onto the circle of radius R5, i.e. given on 
the space x S^. We start from the UV finite theory, i.e. that with Nf ~ 2Nc hypermultiplets. 
The corresponding integrable theory is the inhomogeneous XXZ spin chain |98[ p9| . 

The Lax matrix for the XXZ (sl(2)) spin magnet has the form 

fie^" - fi-'^e-^° 2S'_ 

2S+ /ie~^° - ^l-^e^° 

Thus, this Lax operator is given on a cylinder. In fact, from the point of view of integrable systems, 



some more natural spectral parameter is n — e^, and the Lax operator (3.57) becomes clearly 
trigonometric. This Lax operator is intertwined by the trigonometric r-matrix 

^(C) = ^-j — r (o"! (8) CTi + 0-2 (8) cr2 + cosh7rCo-3 ® 0-3) (3.58) 
smh 7r(, 



so that the Poisson bracket of the Lax operators ( 3.35) gives rise to the Poisson brackets of the algebra 
of Si's: 

{S+, So} = ±5±; {S+, S.} = sinh25o (3.59) 
The second Casimir function of this algebra is 

C2 = cosh 25o + 2S+S- (3.60) 



The non-linear commutation relations (3.59) are ones from the quantum deformed algebra Uq{sl{2)). 
This generalizes the fact that the XXX magnet is described by the Poisson brackets that reproduce 
the classical sl{2) algebra. In these Poisson brackets, the Plank constant h turns out to be inessential 
and can be put equal to unity. In fact, it is proportional to the radius i?5 of the space-time circle in 
the corresponding 5d SUSY theory and can be easily restored with the replace of any generator Si by 
R^Si. Hereafter, we omit i?5 from all the formulas. 

Following the standard procedure, now we consider the chain with sites with the Lax operators 



( p. 57 ) associated with each site and commuting with each other, introduce the inhomogeneities Ci 



which depend on the site of chain by the replace Q Q d and impose the periodic boundary 
conditions. On this chain, the transfer matrix acts and the periodic boundary conditions generate the 
spectral curve and the conservation laws. 

The manifest form of the curve in the XXZ case can be derived using explicit expression for the 
Lax matrix ( 3.57 ) and readsQ 



Q2N. (e2«) 



= 2e-''^'^-T..^^PN^{e'^)^2e-''^^~^S' (^e^^cC + . . . + C.) (3.6I) 



Changing the variables e^^ = /i^ = A, w ^ e ^'w, this curve can be recast in the hyperelliptic 

form in A variables: 



In fact, the coefficient in front of Pnci instead of 2, is equal to 2cosh||, e It is, however, the integral of 
motion that can be put equal to any number. We always fix 2_^^ 5o,z to be zero. 



w - 



Q2N^ (A) 



(3.62) 



while in terms of the "true" spectral parameter C, this curve looks considerably more tricky. However, 
one can work with the variable A taking into account that it lives on a cylinder or sphere with two 



marked points. In equation (3.62), the polynomial 

Q2iv.(A) = n - 2K.,e'^^\ + e^^') ^\{(\ 



•im>: 



defines the masses rnf^' of matter hypermultiplets (cf. with (|3.44| )) 



if = ± ^log {K^ 



K? 



1 



(3.63) 



(3.64) 



Note that, being written in terms of the variable A, the curve (3.62) is very similar to the curves arising 
for Ad theories. However, the difference coming from different generating 1-differentials dS turns out 
to be very crucial (see, for example, the discussion of residue formula and perturbative prepotentials 
in section 5). 



It is easy to compare the curves (3.62) with other curves appeared in literature |llq| , rewriting 
them in the form 



n.sinh(C-m.)^^-|-r^.^^ 



ai) 



(3.65) 



2ai 



where we rescaled w — > e^'''''^^i''^w, denoted the roots of the polynomial Pn^{X) ( 3.61 ) A^ 
and made use of formula ( |3.64 ) and the manifest form of the leading and the constant terms in this 
polynomial. Comparing ( 3.65| ) with (3.61), one finds, in particular, that 



(3.66) 



The condition (3.66) is not, however, absolutely necessary in the context of XXZ chains and 5d 
theories. It emerges only in the standard XXZ chain. In the context of Seiberg-Witten theory 
one needs rather to consider, similarly to the Ad case, the twisted XXZ model [119, |9^, |9^. It is 
characterized by the Lax operator 



fie""-' - - a^^l-'^e~'^'>^'+'^' 2S-^i 
with generally non-unit constants a^. Pi. These constants provide an arbitrary coefficient in front 



(3.67) 



of the product in the r.h.s. of ( 3.65 ) (which can be also fixed by the integral of motion '^f" So,i) 
and, thus, break the condition ( [3.66[ ). They are also important for a careful treatment of dimensional 
transmutation (i.e. decreasing the number of massive hypermultiplets), see below. 

The maximally degenerated case is related to the pure gauge theory when all the masses become 
infinite. In this case, one gets the relativistic Toda chain, s.3.6. 

There certainly exist all intermediate degenerations when all but Nf masses are brought to infinity. 
This system, corresponding to the gauge theory with Nf < 2Nc massive hypermultiplets, is described 



by the chain with some of the sites degenerated and some of them not. The corresponding spectral 
curve has the form 



a i 

or 

^ ^ n«^sinh(C-m„) ^ 2e(^-iv,/2)c+E-'-i/2E"'> fTsinh(C- a.) (3.69) 

i 

(where we again rescaled w — > e^'-^^^^S '""'w; this curve also coincides with that in |118|). As we 



aheady noted, the particular coefficient can be, in principle, replaced by any other 



number, say put equal to unity, which literally corresponds to It does not influence the result 

and we will discuss the reason below. 



Note that the form of the spectral curve ( 3.65 ) is perfectly designed for taking the 4d limit. Indeed, 



one can restore i?5-dependence in this formula multiplying each d and mass parameter by R^. In 



terms of algebra (3.59) it means that each generator has to be multiplied by and so does the 
spectral parameter. Then, one immediately reproduces the results of the previous subsection. 

As in the Ad case, the description of gauge theory in terms of the spin chain allows an extension 
to the group product case. It is described by higher sl{p) XXZ magnets given by the Lax operator 

^(a^) = ® (3.70) 

and 

Lu = Me^°- - ^l''e-^"■\ = 2Sj„ {i ^ j). (3.71) 

Here S^^i are associated with the vectors e^, realizing the simple roots Ui — ei — Ci+i, and has the 
only non- vanishing matrix element The generators Si satisfy non- linear Poisson brackets that 

can be read off from the commutation relations of the quantum algebra Uq{sl{p)) like it has been done 
in ( |3.59 ). These Poisson brackets can be certainly obtained from the quadratic relation (3.35) with 



the trigonometric r-matrix for the sl{p) case |113|. Note that, specially degenerating this system, one 
can easily reproduce the (p, q)-Weh constructed in |l20| that corresponds to different toric varieties. 

Now let us say some words on the generating differential dS in 5c? theories with fundamental 
matter. 

The general scheme presented for M theory can be almost literally transferred to the hd M = 2 
SUSY gauge models with one compactified dimension. It can be described through involving trigono- 
metric r-matrices and L-operators instead of rational ones, i.e. coming from Yangian to affine algebras. 
It means that now it is natural to consider the both parameters Q = \\og\ and logw as coordinates on 
the cylinder. 

Similar arguments imply that, instead of differential dS — \d\ogw, one now has to consider the 
differential |6^, H 

d5 ^ C- - logA^ (3.72) 

w w 



oo 



so that, despite the similarity of the 5d spectral curves with Ad ones, periods of dS'^^'> are different 
from those of dS'^'*-'. Note that the derivatives of this differential w.r.t. moduli again give holomorphic 
differentials. 

In fact, one can interpret the 5d theory as the Ad theory with infinite number of (Kaluza-Klein) 
vector multiplets with masses Af„ = n/Rs and infinite number of analogous (Kaluza-Klein) funda- 
mental hypermultiplets. Then, one can equally consider either the just described 5d picture, or the 
4c? picture that involves Riemann surface presented as an infinite order covering (see the spectral 
curve equation in terms of the variable C ( 3.6l| )) with infinitely many punctures. This latter picture 



can be effectively encapsulated in the usual hyperelliptic Riemann surface ( 3.62 ) with finite number 
of punctures, but, as a memory of infinitely many multiplets, the spectral parameter A now lives 
on the cylinder. Meanwhile, the differential dS^^'' now evidently should be of the form ( |3.72 ) which 
"remembers" of its 4c? origin. 

3.6 5d pure gauge theory: relativistic Toda chain 

We are going to get the pure gauge theory in 5d from the theory with fundamental matter via the 
dimensional transmutation, sending masses of all the hypermultiplets to infinity. Therefore, we study 
the maximal degeneration of the XXZ spin chain This degeneration is obtained from the 

twisted XXZ chain introduced above and is described by the relativistic Toda chain |^ (see also 
[^). The twisting means just introducing some new parameters into the Lax operator, which are 
central elements of the Poisson bracket algebra. A special fine tuning of these parameters allows one 
to match smoothly different regimes and limiting cases of the XXZ spin chain. 

As it has been explained in ||3^ , there are two equivalent ways to twist the integrable system. One 
of them, which was applied in the XXX case, is to multiply the Lax operator by an arbitrary constant 
matrix U. This is possible in the XXX case, since the corresponding (rational) r-matrix commutes 
with the tensor product U ®U . The situation is much more restrictive in the trigonometric case, when 
the r-matrix commutes only with the matrices U of the very special form (say, it can be proportional 
to any one of the Pauli matrices). Therefore, in the trigonometric case we apply the second way of 



doing |11£, 115 1, that is, we consider the Lax matrix of the general form but with some prescribed 



dependence of matrix elements on the spectral parameter. Then, the Poisson bracket ( 3.35| ) dictates 



Poisson brackets of matrix elements and, in particular, fixes some coefficients to be centers of the 
Poisson bracket algebra. 

More explicitly, we fix the Lax matrix to be of the form 

m =1 1 r J (3-^3) 



and to be intertwined by the same r-matrix ( 3.58 ). Then, up to inessential total normalization, one 
can rewrite the Lax operator in the form 

where p, v and 7 are constants (centers of algebra). 

In this case, the commutation relations are slightly changed to 



{S+, So} = ±5±; S^} = - {lye'^" - pje-'''°) 



(3.75) 



and the second Casimir function is equal to 



C2 = ^ [pie 



-2So 



(3.76) 



thus, the determinant of the Lax operator (3.74) is stiU the quadratic polynomial with coefhcient 
being the Casimir function. 

Now we are ready to demonstrate how the limit to the relativistic Toda chain can be done. To 



make the Lax operator (3.74) looking more similar to the relativistic Toda Lax operator, we multiply 
it by the function e"^". In the whole chain it results into the factor which is the integral of 

motion and can be put zero. In particular, in the relativistic Toda case this integral is equal to the 
full momentum of the system. Although being so inessential, this redefining of the Lax operator still 



requires to modify some expressions [115 . For instance, the new Lax operator 



m = 



So 



2S^e 



ip. - vp^^e^^° 



(3.77) 



is intertwined by the new r-matrix (here r is the r-matrix from (3.58)) 



r'*"') = r + i (I as - CT3 ® I) (3-78) 
which is called twisted. This name comes from the quantum generalization of this matrix which can 



be obtained from the standard quantum trigonometric i?-matrix by twisting [115 



^ Fi2{l)R{u)F^^\i)' F12 = F^,^ = exp{i(a3®I-I®(73)} 

An important property of the Lax operator ( |3.77 ) is that its determinant is no longer dependent on 
only the Casimir function but also on e'^". The product of this quantity over all the sites is still 
certainly an integral of motion. 



(3.79) 



Now one can consider the particular reduction [115[ 



I.e. 



7 = i/ = 0, p=l 



(3.80) 



L{p) 



pe~"" — p ^ 2S-e 



o2So _ 



So 



(3.81) 



Taking into account the commutation relations (3.75), that is, {iS'+,S'_} = 0, one can realize the 
algebra as S± — e"*"*. So ~ p, {p, 9} = 1 and reproduce the standard Lax operator for the relativistic 



Toda chain |l2l[[ . 

The spectral curve corresponding to this degenerated case can be easily written 



1 1 

w + - = TTr^Pjv.(A) 



w 



(3.82) 



An 



Now one can also obtain from the theory with Nf — 2Nc the theory with less number of massive 
multiplets degenerating the general XXZ chain at several sites. The curves for such theories we 
already discussed in the previous subsection. 



At last, let us note that the differential dS constructed for the XXZ spin chain (3.72) does not 
change when degenerating to the relativistic Toda chain. 



3.7 5d adjoint matter: Ruijsenaars-Schneider system 

Now we are going to study the system emerging upon adding in 5d theory the adjoint matter hy- 
permultiplet. This UV finite system can be described differently as follows. By starting with a five 
dimensional model one may obtain four dimensional Af — 2 SUSY models (with fields only in the 
adjoint representation of the gauge group) by imposing non-trivial boundary conditions on half of the 
fields: 



(f>{x5+R5) = e^''(l){x5). (3.83) 

If e = one obtains M — 4: SUSY in four dimensions, but when e 7^ 2Tm this is explicitly broken 
to A/" — 2. The low-energy mass spectrum of the four dimensional theory contains two towers of 
Kaluza-Klein modes: 



TTTl 6 'T' TTTL 

M„ = — and Af„ = — — , neZ. (3.84) 

According to the proposal of ||6^], this five dimensional theory may be associated with the elliptic 
Ruijsenaars-Schneider integrable model ^ . In various double-scaling limits it reduces to systems 
we already studied: 



(a) If i?5 ^ (with finite e) the (finite) mass spectrum (3.84) reduces to a single point M„ = 0. This 
is the standard four dimensional M = 2 SUSY pure gauge model associated with the periodic Toda 
chain. As we already know, in this situation M ^ 2 SUSY in four dimensions is insufficient to ensure 
UV-finiteness, thus r — > ioo, but the phenomenon of dimensional transmutation occurs whereupon 
one substitutes the dimcnsionless r by the new dimensional parameter Aq^^ — e^'^" (e/ R^)^" . 

(b) If i?5 ^ and e ^ MRc, for finite M, then UV finiteness is preserved. The mass spectrum 
( ^.84 ) reduces to the two points M„ = and M„ — M. This is the four dimensional YM model with 



J\f = 4 SUSY softly broken to Af = 2. The associated finite-dimensional integrable system is, as we 
know, the elliptic Calogero model. Case (a) is then obtained from (b) by the double scaling limit. 



(c) If i? 7^ but € ioo the mass spectrum (5.84) reduces to a single Kaluza-Klcin tower, 
Mn — irn/R^, n E Z. This compactification of the five dimensional model has JV = 1 SUSY and is 
not UV-finite. Here r ioo and e ~> ioo, such that 2ttt — Nee remains finite. This system as we 
discussed in the previous subsection is described by the relativistic Toda chain. 

(d) Finally, when i?5 7^ and e and t arc both finite one distinguished case still remains: e = tt/2^ 
Here only periodic and anti-periodic boundary conditions occur in the compact dimension. This is 



the case analyzed in [122 1. It is clearly special from both the point of view of Yang-Mills theory and 
integrable systems. 

■^■^The case e = of fully unbroken five dimensional J\f = 2 supersymmetry is of course also distinguished, but trivial: 
there is no evolution of effective couplings (renormalisation group flows) and the integrable system is just that of Nc 
non-interacting (free) particles. 



Thus, the eUiptic Ruijsenaars-Schneider model is a remarkable completely integrable system whose 
various limits include the (finite dimensional) Toda and Calogero-Moser models. Here we shall review 
only a few of its salient features. More comprehensive accounts of its structure and applications may 



be found in |123, 124|. The SU{Nc) model has an explicit Lax representation with Lax operator 
[pT[ p2[ |125[ p6| 



= c(C|e) 



(3.85) 



Here 



(3.86) 



where the pi and qi are canonically conjugate momenta and coordinates, {pi, qj} — Sij. The commuting 
Hamiltonians may be variously written 



E n 



JC{1.....^} 



F(^qij\e) 



E 



■+Ph- 



l<ii<...<ik<N^ 



n 

a<b 



1 

p{qiaib) - P(e)' 



(3.87) 

(The final product is taken to be unity in the case k = 1.) These Hamiltonians arise in the description 
of the spectral curve. Note that, after trigonometric degeneration, these Hamiltonians are proportional 
to those obtained in s.2.2. We discuss this degeneration in more details in section 5. 



The integrability of the model depends on the functional equations satisfied by F fclj. 



the 



connection with functional equations and integrable systems is part of a larger story [127| 



Some comment on the parameters ^ and e in these formulae is called for. Here ^ is precisely the 
spectral parameter associated with the bare spectral curve, where the Lax operator is given. Further, 
the additional parameter e in (3.85), ( p. 86 ) above is the same parameter we introduced in (3.82) 



characterizing the boundary conditions. Actually the integrability of (3.87) does not require e to be 
real, but such a choice guarantees the reality of the the Hamiltonians. The identification of these 
two parameters has been simplified by our choice of the real period of p being tt so that both ( 3.83| ) 
and ( [3.85 ), (3.86) are manifestly tt periodic. The "non-relativistic limit" e ^ which leads to the 
Calogero-Moser system means we can identify the mass M = e/i?5 of the gauge multiplet with the 
Calogero-Moser coupling constant. The special point e = tt/2 oj) singled out earlier is now a 
half-period of the p-function, and at this point p'(§) = 0. 

Now inserting the Lax operator ( 3.85 ) into the determinant (^ 



detjv.xjv. (C^'iO - A) = 



(3.88) 



one gets the spectral curve equation. Its coefficients are expressed through the Hamiltonians ( 3.87 ) 
and parameterize the corresponding moduli space. 

To conclude constructing Seiberg-Witten theory we should determine the differential dS. It is 
evident from our previous discussions that it is to be of the form 



dS = —logXd^ 

^5 



(3.89) 



,10 



Indeed, on one hand, it should be proportional to since describes the adjoint matter (compare with 
the Calogero case). On the other hand, the proportionality coefficient is to be i?5 ^logA, since the 
Ruijsenaars-Schneider system describes a 5d theory. 

Now let us discuss more the spectral curve arising. From the point of view of classical integrability 
the overall normalization factor c(^|e) of the Lax operator does not change the integrals of motion 
apart from scaling. This normalization factor does however lead to the rescaling of A and this effects 
the explicit form of the 1-form dS*.]^ We shall see below that significant simplifications occur with 
the choice 



cm = com ^ , , , - / /^^!w 4-.^ = ^c^mc-m 



(3.90) 



where 



^ cr(£)cr(e) 1 



Similar issues of normalization enter into discussion of separation of variables 128 |. 

Before generalizing to the higher rank situation it is instructive to first consider the case of gauge 
group SU(2). In the formulas of the previous section this corresponds to setting J\f = 2 and working 
in the center-of-mass frame pi + P2 = 0. We define p = pi and q = qi — q2. 



Using the explicit form of the Lax operator the spectral curve (3.88) is found to be 



A^ - cuX + c^ipiO - P(e)) = A^ - cuA + ^ = 0, (3.92) 

Co 



where 



u = Hi (p, q)^2 coshpv/p(9) ^ p(e) (3-93) 



is the Hamiltonian (3.87). For the choice c(^|/i) = co(^|/i) this simplifies to yield 

A^ - MCoA +1 = 0, (3.94) 

or simply 

Observe that with our choice ( 3.90| ) this equation, which describes the Seiberg-Witten spectral curve, 
may be expressed in the form 

u = H{\og\i). (3.96) 
Comparison with ( 3.89| ) shows that the generating differential Rc,dS — logAd^ takes the iovm pdq here. 



We now consider the general SU{Nc) model. The first step is to evaluate the spectral curve (3.8S 



for the Lax matrix ( 3.85 ). Expanding the determinant about the diagonal yields 



Note that although the choice of normalization c(^|e) (or, equivalently, rescaling A) can change the manifest 

d\ 
A 



expression for dS, it affects neither the symplectic 2-form ^ A — , nor the period matrix, which is the second derivative 



of the prepotential. 



AO 




e^-+-+^-det(,,)lA:i^^=0, (3.97) 

l<ii<...<ifc<-/Vc 



and Ruijsenaars |6]J, p2] expressed the determinants appearing here in terms of the Hamiltonians 



( p. 87] ) by means of a generahzed Cauchy formula. This expansion was re-obtained in ||69| using a 
simple fact that generalized Cauchy formulae can be derived in terms of free-fermion correlators via 
Wick's theorem for fermions. The Ruijsenaars-Schneider model (and its spin generalizations) may 



be understood |12q , 129 1 in terms of a reduction of the Toda lattice hierarchy. The Hirota bilinear 
identities of that hierarchy may be expressed in terms of free-fermion correlators [130, |2^, and are 
the origin of those here. We believe that the free field expansions will ultimately lead to a better field 
theoretic understanding of the appearance of these integrable systems. The machinery of free-fermion 
correlators has already found use, for example, in calculating within the context of the Whitham 
hierarchy for these integrable systems |]8^ . 

The formulas for the determinants are 



Dr 



detf 



(-) 



M^-i) a''-^^ - e)cr(e + (fc - l)e) 



n 



<7^(9i„ijCT^(e) 



(3.98) 

where Ik — {ii, . . . , ik\ denotes the set of k indices and a, 5 — 1, fc. The right hand side of this 
expression is in fact a doubly periodic function of all the arguments ((/i , ^ and e) and so may rewritten 



in terms of Weierstrass functions. In particular, upon using (3.24), we find that 

'^'(9.„.J'^'(e) T-r 1 



n 



n 



(3.99) 



thus recovering the Hamiltonians (3.87) noted earlier. 

Now we are ready to obtain the final formula for the spectral curve in the Ruijsenaars-Schneider 



system. Upon substitution of (3.98) into (3.97) the curve takes the form 



N 



(3.100) 



fe=0 



fe=0 



Here we have collected the g-independcnt factors in (3.9S) into where 



Dkm = (-) 



k(k~^ gfc-i(^-e)g(^ + (fc- 1)£) 
cr''^(^)aMfc-i)(e) 



l_ a(C + (fc-l)e) i- 

r,k 



(T(e-e) ■afc('=^-2)(e) 



(3.101) 



Again this is doubly periodic function in both ^ and e and so is expressible in terms of the Weierstrass 
function and its derivative, though the explicit formulae are rather complicated. Although it may 



appear from the last expression in (3.101) that (3. IOC) simplifies when c = c_, such a choice would 
break the double periodicity of our spectral curve and so is inappropriate for the fully elliptic modelf^. 

Part of the difficulty in dealing with the elliptic Ruijsenaars-Schneider and Calogero-Moser models 
is the complicated nature of these spectral curves. For example, in the case of SU{i) the spectral 



curve (|3.88| ) is 



^■'However, in the perturbative limit when one of the periods becomes infinite this choice is then available. 



- cu\^ + c^v\{p{^) - p(e)) + 
Qp'(Op'(e) - 3p(0p'(e) + P'(e) + 2ap(0p(e) + 2/3(p(0 + p(e)) + 27 j = (3.102) 
u = Hi = H+,v = H2 = H^ 



where 



H± = e^P^ V P(9i2) - P(e) Vp(9i3) - p(e) + e=^^' V P(9i2) " p(e) Vpfea) - p(e)+ 

(3.103) 

+e=^P3 \/ p(gi3) - p(e)\/ p(g23) - p(e), Pi+P2+P3 = 0. 

Again the issue is whether some choice of the normahzation c might simphfy matters. With our choice 
of ( |j.90 ) for example, the third term in ( 3.102| ) turns into just vX, but no drastic simplification occurs 



in the constant term unless at the special point e — (where p'(e = tt/2) = 0) when - modulo c - 
the whole equation becomes linear in x = p(^). 

To conclude the discussion of the Ruijsenaars-Schneider system, we separately consider the special 
point where e equals the real half-period lu — Tr/2. Using the fact that cr(^ + 2uj) = — cr(^)e^''^^^"'"'^-' 



with 771 = ^{lo) [10£], we find that 

EE c+(ek) - acoi^lu;) = a'c^i^W) (3.104) 

with a = e^'"''^ . Then 

r c^'^H^'^"-'''/' k even, 

DMe) = <^ (3.105) 
I c^'H^'^'-'^^'coH;'/' k odd. 

Here Tii = 2ef + 6263 is independent of ^ and we note that in the trigonometric limit Tii — 1. The 
choice c ~ Co now essentially removes all ^-dependence in the spectral curve apart from a term linear 
in Co that multiplies the odd Hamiltonians. Similar to the SU{3) case, we choose c — —icq. By 
absorbing the (inessential) constant factors into the Hamiltonians, hk = i^Tii ^'"^ '^^^'^Hk for even k 
and hk = i''-^n^^^^n^''''''~^^^^Hk for odd k, one obtains the non-perturbative spectral curve ( pocj ) 
at the special point e — lu: 

EitM-A)"-" _P(A)-(-)"P(-A) 



where 

P(A) = 5:/^.(-A)^-'= = A^-t-... + l (3.107) 

k 

These formulae may also be easily expressed in terms of the Jacobi elliptic functions. For example, 
p(6 - P(w) = (ei - e3) ™:i[";fc] and, therefore, co(C|'^) = ^/=t^^!^^i here sn{u, k) and cn{u, k) are 
the Jacobi functions [109|, u = (ei — 63)^/^^, k is elliptic modulus. 



To conclude this subsection, note that one can immediately generalize the Ruijsenaars-Schneider 



model to other groups, see |131 



A C 



3.8 Fundamental matter in Qd: XYZ chain 



At the last subsection of this long and detailed consideration of different Seiberg-Witten theories, we 
make the first step towards 6d theories, that is 6d theory with the fundamental matter. This is the 
theory with two extra compactified dimensions, of radii R5 and Rq. The same argument as we used 
for the transition from four to five dimensions, i.e. taking account of the Kaluza-Klein modes, allows 
one to predict the structure of results in the 6d case as well. Indeed, one expects that coming from Ad 
(5d) to 6d means one should replace the rational (trigonometric) expressions by the elliptic ones, the 
(imaginary part of) modular parameter being identified with the ratio of the compactification radii 

This would give rise to " elliptization" of an integrable systems. In terms of 2 x 2 Lax representation. 



the "elliptization" corresponds to replacing by elliptic r-matrices and L-operators the rational (3.34) 



and trigonometric ( |3.57 ) ones. In particular, it means that now it is natural to consider the parameter 



C = ^logA as co-ordinate on a torus, i.e. one of the spectral parameters lives on a torus, while the 
other one remains on a cylinder. 

In fact, we know the only system providing proper elliptization of the XXX spin chain, that is, 



the XYZ Sklyanin chain [132 . It was propose as a candidate for the integrable system behind the 6d 



theory in |97, BSl 



Now we briefly describe the XYZ chain. All necessary facts on the elliptic functions can be found 



in [109|. The Lax matrix is now defined on the elliptic curve E{t) and is explicitly given by (see 113 



and references therein): 

3 

iSki(^) ^sH + tY, WaiOS'^aa (3.108) 



0=1 

where 



To keep the similarities with |97|, |93, l99|] , we redefine in this subsection the spectral parameter C i— 



2K' 

where K = , = ?£'?(0), fc^ = so that fsT ^ f as r ^ ioo. This factor resuhs into 

the additional multiplier tt in trigonometric functions in the limiting cases. 



The Lax operator (3.108) satisfies the Poisson relation (3.35) with the numerical elliptic r-matrix 



r(C) = — I X]a=i ^a(C)o'a ® cTq , which implies that 5*°, 5" form the (classical) Sklyanin algebra [132|: 

=i(eb-ec)5^^= 

(3.110) 

[S'^^S^] = -iS^S" 

with the obvious notation: ahc is the triple 123 or its cyclic permutations. 

One can distinguish three interesting limits of the Sklyanin algebra: the rational, trigonometric 
and double-scaling limits |M. We are interested in only trigonometric limit here, since it describes 



degeneration to the 5c? case. In this limit, t +i(yo and the Sklyanin algebra (3.110) transforms to 

{S^^S^"} = 0, {S^ ,S°} = iS^S^ , {S^ ,S^} ^ ~2iS'^S^ 
^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ (3.111) 

{S^,S^} = -iS'^S^, {S^,S''} = iS''S^, {S^,S^} = -iS''S^ 



The corresponding Lax matrix is 



Lxxz = S°l - (§Vi + S''a2 + coshTrCS^aa 

sinnTTC V 



(3.112) 



and the r-matrix 



sinh TiC, 



(cFi (g) CTi + (72 0-2 + cosh 7rCo-3 (Ki 0-3) 



(3.113) 



Now one can note that the algebra ( p.lllD admits the identification So = e + e~ , 5*3 



Qtrig Qitrig 

e~ — e . With this identification and up to normahzation of the Lax operator (3.112), we 



finally get the Lax operator ( 3.5? ) of the XXZ chain with the Poisson brackets (|3.5 
The determinant det2x2-^^'*'(C) is equal to 



det2x2i''^HC) ^S'o+J2 - Y.Sl^K^ M^piO ^K~M' 



(3.114) 



where 



a=l 



(3.115) 



are the Casimir functions of the Sklyanin algebra. 

Now, in. order to construct the spectral curve, note that the determinant of the monodromy 
matrix T^^ {(,) is 



AT, 



Q(C) = ciet2x2rA^.(C) = ndet2x2^"'^'(C - G) = n (^^ - ^^^^(^ - ^^)) 



i=l 



(3.116) 



Qo = Y[A, 



9.(c-m<+')e.(c-m<-') 



where we used the Appendix of [g9| in order to rewrite Ki — Mfp{( — Q) as Ai — ' g'i^^^_^:^ 

Here 



and Ki is a solution to the equation 



Similarly to the 5d case. 



(3.117) 



(3.118) 



Nf 



(3.119) 



As a particular example, let us consider the case of the homogeneous chain |97| (all Q — in (3.116)). 
Trr;v^(C) is the combination of polynomials: 



P(C) = Pol\'-lJx)+yPol\%_,.{x), 



(3.120) 



where [^] is integral part of and coefficients of Pol^^^ and Pol'^^^ are Hamiltonians of the XYZ 
model. The spectral equation ( |3.13 ) for the XYZ model is: 



(3.121) 



where for the homogeneous chain P and Q are polynomials in a; = p(C) ^-^d V — \p'{C)- Eq. ( p. 121 ) 
describes the double covering of the elliptic curve E{t): with generic point C, G E{t) one associates 



the two points of Y.-^^^ , labeled by the two roots w± of equation ( 3.121 ). The ramification points 

±VQ, or r i 



correspond to w+ = w- = ±y/Q, or Y ~ ^ [j^ ^ w) ^ V-P^ — Q = 0. Note that, for the curve 
( |3.121 ), X = oo is not a branch point, therefore, the number of cuts on the both copies of E{t) is Nc 
and the genus of the spectral curve is iVe + Ipl 

In the case of a generic inhomogeneous chain the trace of monodromy matrix looks more compli- 
cated, being expression of the form 



Trr^,(C) = Tr [] (C " CO = Tr ^ [] ^-.^^"W„.(C - G) 



Na>i>l 



{a,} N^>i>l 



(3.124) 



a, = 0,a,; a, = 1,2,3; Wo(C) = 1 

It is clear from ( ^.109 ) that each term in the sum ( 3.124 ) (and thus the whole sum) has simple poles 
in all inhomogeneities C = Ci- Thus, the trace of the monodromy matrix can be represented as 



P(C) = ^TrT..(C) = Pon|lML^ 
E a* = E = I E "^a 



(3.125) 



(3.126) 

i i a 

where should be identified with (shifted) gauge moduli. The equation ( p. 121 ) can be finally rewritten 



as 



^ + vo — =r;v; — —. — r — ^o: 



wn:=i^*(c-c. 



nf=i^*(c-o) 



(3.127) 



We see that in the 6c? case it is necessary always to require Nf — 2Nc in order to have a single-valued 



meromorphic functions on torus. We also meet again the constraint (3.126). In the 5d case, it was 



'Rewriting analytically E'^ 



the set of holomorphic 1-difTerentials on S 

x°'dx 



as a system of equations 
3 

y^ = Y[{x-ea), Y^ = P^ 

a = l 

can be chosen as 



■Q 



dx 

V = , Va 

y 



yY 



= 0, 



yi3 



r^dx 
Y 



13 = 0,. 



Nc 



(3.122) 



(3.123) 



with the total number of holomorphic 1-differentials 1 -I- ( ["^] ^" -'^) + ( [ ^"2 ^ ] ~^ ~ ~^ ^ being equal to the genus 
of S^yz. 



removed, along with the condition Nf — 2Nc twisting the spin chain. Since no twisted XYZ chain is 
known we are forced to preserve both constraints^ 



The curve ( 3.127 ) has no longer that simple expression in terms of polynomials of x and y similar 



to ( 3.1201 ). Thus, in the inhomogeneous case it is not convenient to work in terms of "elliptic" variables 



X and y (analogous to A in the 5d case) instead we deal with the variable C G C/r. 

Given the spectral curve and integrable system, we complete the description with writing down 
the generating 1-differential dS. It can be evidently chosen in the following way 

dS^YZ c^^^^ -logwdC (3.128) 
w 

since, as any spin chain, it is to be proportional to dw/w and to C due to six dimensions (in 5 
dimensions, C was just logA and in 6 dimensions A is replaced by the pair (x, y) related with ( via 
p(C) and p'(C))- Return to our oversimplified example of the homogeneous chain. Now, under the 
variation of moduli (which are all contained in P, while Q is moduli independent), 

,(d5--) ^ ^dC = ^=:PL==dC = ^SP (3.129) 
' w ^ ^P(C)2 - 4Q(C) yY ' 



and, according to (3.123), the r.h.s. is a holomorphic 1-differential on the spectral curve (3.121) 



4 Integrability and J\f = 1 supersymmetric YM theories 

4.1 Superpotentials in J\f = 1 theories from J\f = 2 spectral curves 

In this section we shall discuss gauge theories with Af = 1 SUSY and their relation with the M = 2 



SUSY theories. For the detailed review of A/" = 1 theories we recommend [133|. Here we explain some 
features of A/" = 1 theories important for treating them within the integrable framework. 

The key difference between theories with M — 1 and TV = 2 SUSY is the absence of vacuum valleys 
in the M — 1 case. Therefore, the Coulomb branch of the moduli space is absent and the number 
of the vacuum states is finite. The order parameter which distinguishes between the vacua can be 
identified with the gluino condensate. Due to the discreteness of vacua the spectrum of BPS states 
now is different. Instead of the particles saturating the BPS bound one has domain walls interpolating 
between vacua as well as stringy objects which can be identified with the strings themselves or with 
the domain wall intersections. These objects turn out to be related with the corresponding central 



charges in the algebra oi Af = 1 supersymmetry |134, 135] 



The main target is to be an object that presumably could be determined exactly. Instead of 
the prepotential in A/" = 2 theory, the function which enjoys holomorphy is the superpotential W; 
therefore, one could attempt to evaluate it. There are generically perturbative and nonperturbative 
contributions to the superpotential. As for the perturbative contribution, the /3-function unlike the 
Af = 2 case gets contributions from all loops. However, the higher loops IR contributions can be 



calculated exactly [ 136 |. There is also a crucial difference between nonperturbative contributions in 



Af — I and Af ~ 2 theories. Indeed, configurations of all possible instanton numbers contribute 



26 T 



Formula (3.126) implies that, while the sum of all the would-be gauge moduli still remains a constant, it is no 



longer zero. Therefore, one would rather associate with the gauge moduli the quantities shifted by the constant 
2^ ^ma- We do not know the proper field theory interpretation of this condition. 



An 



into the prepotential in the M — 2 case, while only one-instanton terms contribute into the A/" = 1 
superpotential due to the non- anomalous R-symmetry. 

Despite all this difference one may expect links between the exact results in A/" = 2 theories and 
the superpotential. These expectations are true indeed and below we sketch the relation following 



[ [137| and especially |138|. The starting point is the Seiberg-Witten curve P{x,y) = in A/" = 2 
theory with some matter content. Then, one should fix singular points on the Coulomb branch of the 
moduli space which correspond to the vacua in A/" = 1 theory. These singular points are in one-to- 
one correspondence with vanishing of the discriminant of the curve. Then, comparing the equations 
in A/" = 2 theory and the equation dW — 0, where the derivative is taken w.r.t. the appropriated 
composite field, one observes that these equations coincide if the coordinate on the spectral curve is 
identified with some composite operator. 

Consider J\f = 1 supersymmetric SU{2) gauge theories in four dimensions, with any possible 
content of matter superfields, such that the theory is either one-loop asymptotic free or conformal. 
This allows one to introduce 2Nf matter supermultiplets in the fundamental representation, Qf, 
i = l,...,2Nf and Na supermultiplets in the adjoint representation, <i>°'', a — 1,...,Na,- Here a,b 



are fundamental representation indices, and = The numbers Nf, Na are limited by the 



condition to keep the theory asymptotically free or conformal: 

bi = 6- Nf ^2Na>0, (4.130) 
where —bi is the one- loop coefficient of the gauge coupling beta- function. 

The effective superpotential of an = 1 supersymmetric SU{2) gauge theory, with 2Nf doublets 
and Na triplets isQ 

2^ l/(4-hi) 



WNf,NAM,X,Z) = -(4- 6i){A-^ipf27v,^[detjv^(r„^)] } 
-t-trAT^/^Af + ^tr2NfmX + --^tr2Ar^F"Za, 



(4.131) 



where 



ro.p{M,X,Z)=Mo.p+tT2NAZaX-^ZpX-^). (4.132) 



The first term in (4.131) is the exact non-perturbative superpotential, while the other terms give the 
tree- level superpotential. A here is a dynamically generated scale, while 'mij and F^" are the bare 
masses and Yukawa couplings, respectively (^q/3 = /^/3a, rriij — —rriji, = Y^). The gauge singlets, 
X, M, Z, are given in terms of the N = 1 superfield doublets, Q", the triplets, as follows: 

X^■j ^Q^aQj,a=l,2,;i,j = l,...,2Nf, 
M^p = <i>S6*^a; a, /3 = 1, ^a; a, 6 = 1, 2, (4.133) 

The gauge-invariant superfields Xij may be considered as a mixture of SU (2) "mesons" and "baryons" , 
while the gauge-invariant superfields Zf^ may be considered as a mixture of SU (2) "meson-like" and 
"baryon-like" operators. 

As an example, consider the theory with Nf — 1, Na = 1. In this case, the superpotential reads 



''Throughout this section, for the sake of brevity we omit the index QCD from the notation of h-QCD- 



pf 1 1 

Wi 1 = -^r^ + fiM + -trmX + -^tiYZ (4.134) 

Here m, X are antisymmetric 2x2 matrices, Y , Z are symmetric 2x2 matrices and 

T = M + tr{ZX-^y. (4.135) 

To find the quantum vacua, we solve the equations: OmW = dxW = dzW ~ 0. The equations 
dW = can be re-organized into the singularity conditions of an elliptic curve: 

y2 ^ x3 + ax^ +bx + c (4.136) 

(and some other equations), where the coefficients a,5, c are functions of only the field Af, the scale 
A, the bare quark masses, m, and Yukawa couplings, Y . Explicitly, 



where 



A3 a 
a = -M, b^—Pfm, (4.137) 



A6 

a = — detF. (4.138) 
The parameter x in the elliptic curve is given in terms of the composite field: 

xeeV (4.139) 

yVi,i has 2 + Nf = 3 vacua, namely, the three singularities of the elliptic curve. These are the three 
solutions, M{x), of the equations: — dy^/dx = 0; the solutions for X, Z are given by the other 
equations of motion. 

Similarly, if one considers a generic SU{Nc) group, the spectral curve for the generic Nf, Nc can 



be put into the form ( 3.41 ) 



y' ^P^^iX)+QNf{X) (4.140) 



with A = ZX-^ [137|. 



From this consideration, we can get an important lesson for the integrability approach; that is, 
one can make an identification of one of the coordinates defining the spectral curve. Indeed, the 
Lax operator of the spin chain (when acting on the eigenfunction, Baker-Akhiezer function) can be 
identified with some known composite operator. To get the pure gauge theory one has to decouple the 
fundamental matter sending its masses to infinity. Therefore, the Lax operator in the Toda system 
can also be identified as a composite operator which involves very heavy "regulator" fermions in the 
fundamental representation. The identification of the second coordinate involved in the equation of 



the spectral curve (4.140) is more tricky. 



Let us note that the relation between JV = 1 superpotentials and J\f = 2 curves can be established 



moving in the opposite direction as well |139|. To this aim one can consider the Affieck-Dine-Seiberg 
superpotential supplemented by the contribution from the heavy adjoint hypermultiplet and use the 
Konishi anomaly relating the vacuum expectation values of different operators 



CI 



1 1 1 TrW2 

^ ' (4.141) 

/ 1 TrT4^2\ 

( 2tr ^iM + -=trrZ + 2 — ^ ) = 



\ ^ V2 8^2 ^ 

where the vacuum expectation value of Tr is the gluino condensate. The combination of the 
equation dW — and the Konishi relations amounts exactly to the condition of vanishing discriminant 
of the Seiberg-Witten curve. 

Generically one can consider more specific degenerations of the Seiberg-Witten curves resulting 
into specific Af = 1 theories. For instance, one can consider all B-cycles on the Af — 2 curves vanishing 



simultaneously |14C(| . In this case, the very peculiar spectrum of strings emerges in A/" = 1 theory [|141 |. 



One more interesting possibility comes from the phenomenon of collisions of vacua at some points in 



the parameter space [142| which is expected to correspond to nontrivial conformal field theories in four 



dimensions |143|. It was recently shown that, using the correspondence between Af — 2 and Af = 1 



theories, one can identify the phase transition at these Argyres-Douglas points as the deconfinement 



phase transition |139|. The analysis of the Argyres-Douglas points within the integrability approach 



remains an interesting problem. 

4.2 Superpotentials from nonperturbative configurations 

We used the analysis above to recognize some meaning of the Lax operator. Now let us make the next 
step and identify the configuration space of the dynamical system under consideration in terms of 
the field theory. This can be done via calculating superpotentials from nonperturbative configurations 



[ 144 , 145 , 146 1 . In the brief description below we follow the paper 1 144 1 . To have the convenient starting 
point, one compactifies the four-dimensional theory to three dimensions on a circle of radius R. The 
most obvious modification is that the theory acquires new scalar degrees of freedom which come from 
the Wilson lines of the four-dimensional gauge field swapped around the compact dimension. Non-zero 
values for these scalars generically break the gauge group down to its maximal Abelian subalgebra, 
SU{Nc) — > {/{l)^""^. Hence, at the classical level, the compactified theory has a Coulomb branch. 
By a gauge rotation, the Wilson line can be chosen to lie in the Cartan subalgebra, (p — ip'^H"'. Each 
scalar (j)"" is related to the corresponding Cartan component of the four-dimensional gauge field as 
0° = (jigi A" • dx. As a consequence, 0° can be shifted by an integer multiple of 27r by performing a 
topologically non-trivial gauge transformation which is single- valued in SU{Nc) / Zjq^. Thus we learn 
that each is a periodic variable with period 27r. 

In addition to the Nc periodic scalars (j)"" , the low energy theory on the Coulomb branch also 
includes Nc massless photons with — 0, where i = 1, 2, 3 is a three-dimensional Lorenz 

index. In three dimensions, a massless Abelian gauge field can be eliminated in favour of a scalar by a 
duality transformation. Thus we can eliminate the gauge fields , in favour of Nc scalars a"" subject 
to X^iii = 0- particular, the fields appear in the effective action as Lagrange multipliers 
for the Bianchi identity: Ca ^ ae^^^diFjk- Specifically, the normalization of ct" is chosen so that 
the resulting term in the action is precisely —ia°'k°' where fc" is the magnetic charge in the U{\) 
subgroup corresponding to the Cartan generator The Dirac quantization condition implies that 
each magnetic charge is quantized in integer units: k'^ G Z ^ and so the path integral is invariant under 
shifts of the form cr" ct" + 27rn" with n°- E Z. Thus, like 0°, each ct" is a periodic variable with 
period 27r. 



It is natural to combine the real scalars (j) and a to form a complex scalar X = —i{a + rt/i) |147 



where r — Airi/g'^ + 9/2Tr is the usual complexified coupling. Let us now concentrate on the Euclidean 
theory with gauge group SU{2). Gauge field configurations in the compactified theory are labelled by 
two distinct kinds of topological charge. The first is the Pontryagin number carried by instantons in 
four-dimensions, 

P = A / Tr[FAF] (4.142) 

The 4d instanton number appears in the microscopic action as the term —ip9. An important feature 
of the compactified theory which differs from the theory on i?* is that p is not quantized in integer 
units. However, for each integer value of p, the theory has solutions with action Sd = 87r^|p|/f/^ — ip9. 
Thus we have Sd — ^nipr for p > and Sd = — 27ripT for p < 0. When the instanton scale size is 
much less than the compactification radius, these solutions are close to their counterparts on i?^. 

In addition to the usual four-dimensional instantons, the compactified theory also has finite action 
configurations which carry three-dimensional magnetic charge, 

I W -B e Z (4.143) 

where Bi — SijkF^^ /2 is the Abelian magnetic field of the low-energy theory. As these configurations 
do not depend either on the time or on the coordinate of the compactified spatial dimension, we will 
refer to them as three-dimensional instantons. In addition to the integer-valued magnetic charge fc, 
these Sd instantons also carry fractional Pontryagin number p = cf) /2tt . The magnetic charge appears 
in the action as the term —ika while, as above, -id instanton number appears as the term ~ip9. 
Including these term the action of a 3d instanton can be written as 5*3^ — —iT(f)\k\ + ika. Thus we 
have Ssd — kX for fc > 0, and S^d = kX for fc < 0. 

Putting everything together we can write out the most general possible superpotential which can 
generated at leading semi-classical order, 

W = 171117121713 ^ UpQ^ + TO1TO2TO3 

(4.144) 

00 00 

+mim2m3 ^ ^ Cfc^nq''" exp(-fcX) + (ife,„g''" exj>{+kX) 

k=l n=l 

The dependence of the superpotential on the masses of the adjoint scalars rrii is fixed by the zero 
modes counting. 

The superpotential of the SU (2) theory is a holomorphic function on the chiral superfield X = 
—i{a + Tcf)). The weak coupling arguments indicate that this function should respect the classical 
periodicity of the variables and a. Thus, W(X) is a holomorphic function defined on the complex 
torus E. The superpotential should also be invariant under the Weyl group of SU{2) which acts as 
X —X. As >V(X) is a non-constant holomorphic function on a compact domain it must have a 
singularity somewhere. As usual such a singularity should correspond to a point at which our low- 
energy description of the theory breaks down due to the presence of extra light degrees of freedom. 
There is only one such point on the classical moduli space: the origin X — Q where non- Abelian gauge 
symmetry is restored. Thus we will assume that W has an isolated singularity at this point and no 
other singularities. Finally, as there is only one elliptic function of order two with a double pole at 



the origin we have the result 1 144 , 



W{X) = TOimama {V{X) + C(t)) (4.145) 
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Importantly, we have obtained a superpotential which is independent of R. 

Starting from an Af — 2 theory in four dimensions, the superpotential of the M — 1 theory on 
X obtained by introducing a mass for the adjoint scalar in the Af = 2 vector multiplet should co- 
incide with the complexified potential of the corresponding integrable system. This connection should 
be a general one which applies to all A/" = 1 theories obtained by soft breaking oi J\f — 2 supersym- 
metry. The key point is that the Coulomb branch of an A/" = 2 theory on i?'^ x coincides with 
the complexified phase space of the integrable Hamiltonian system associated with the corresponding 
theory on R^. The dynamical variables of the integrable system, {X^jP"} yield holomorphic coor- 
dinates on the Coulomb branch of the compactified theory. The dimensional transmutation provides 
then the corresponding superpotential for the Toda case too. The superpotential for the Toda case 
has been also derived from the brane picture for instantons as wrapped branes [146|. 



4.3 J\f = 1 theories and integrability 

We complete this section with general comments concerning the relation between integrable dynamics 
and Af — 1 SYM theories. First of all, one should expect that superpotentials in A/" = 1 theories 
are associated with Hamiltonians of integrable systems^. Therefore, the vacuum in Af = 1 theory 
correspond to a stationary state of the dynamical system 



dW/d4)i ^ dH/dxi = (4.146) 

and the number of vacua in Af = 1 theory coincides with the number of stationary states. Note that 
the locus corresponding to the stationary states of the elliptic Calogero system plays the essential role 
in the relation between the Calogero system and the KdV hierarchy. Namely, the initial data in the 
KdV hierarchy upon this correspondence should follow from the equilibrium locus in the Calogero 
system. 

Let us now say how to come from the Af = 2 theory to the Af = 1 one in integrable terms. First, 
the Af = 2 system has to be placed at the singularities of the Coulomb branch. This means that some 
integrals of motion in the integrable system have to be fixed. If we are sitting at the point of maximal 
degeneracy, all the integrals of motion have to be fixed and the finite-gap solution degenerates to the 
solitonic one. This has been carefully done in the paper [148|. At the next step, one has to give a mass 
to the adjoint scalar that breaks Af = 2 SUSY down to Af = 1. The spectral curve is then "rotated" 
within (C-P'^ and is given by two polynomial equations 



P(w, u, x) = 0, Q{v, u,x) = Q (4.147) 

Peculiar expressions for these polynomials for different J\f = 1 theories are presented in fs^ . These 
equations definitely mean that the integrable system is to be modified. However, the question of 
possible interpretation of the spectral curve of in A/" = 1 theory as a spectral curve of some integrable 
system still remains open. 

We have discussed above the (very restricted) role played by the Af = 2 integrable system in the 
context of A/" = 1 theory. However, it seems that the Whitham integrability (see section 6 below) 
is much more prominent candidate for describing the Af = 1 theories, since it actually governs the 

^*The idea of this identification comes from the fact that supersymmetry is basically broken by the tr<3?^-tcrm whose 
vacuum value is the Hamiltonian of the corresponding Af = 2 integrable system (we discuss the SU{2) case), see also 
below. 

t; A 



dependence of the effective action on Af — 2 SUSY breaking parameters. In a sense, the Whitham 
theory which is nothing but RG flows predicts a very pecuhar "time" -dependence consistent with RG 
flows |149| . All the "times" but the discrete one are SUSY breaking parameters and some predictions 
concerning even Af — theories can be madePj. 



5 Perturbative prepotentials 

After a short discussion oi Af ~ I SUSY theories, we return to the main subject of this survey, Af — 2 
theories. 

In previous sections we established connections between Seiberg-Witten theories and gauge theo- 
ries. However, so far we studied only a kind of general landscape in the theory, since the real, numerical 
results enter the game only with the prepotential calculated. Note that even the correspondence be- 
tween gauge and integrable theories is supported by comparing the perturbative prepotentials. We 
discuss this issue in this section. 

In fact, as we noted already, calculating the prepotential is quite a challenging task. It can be 
done, mostly in the framework of integrable systems, order by order. However, technically one needs 
to use different specific properties of prepotentials to simplify these calculations. 

The other essential point is that the prepotential celebrates many properties that allow one to 
push forward the theory in many directions. This makes the prepotential an important object from 
both mathematical and physical standing points. 



5.1 Perturbative prepotentials in gauge theories 

We start with the perturbative prepotentials, since, on one hand, one is able to compare the pre- 
potentials calculated in physical theories and in integrable systems only at perturbative level and, 
on the other hand, many important properties of the prepotential are observed already at this level. 
Therefore, now we are going to give an explicit description of the perturbative prepotentials for the 
Seiberg-Witten theories discussed in section 3. 

The technical tool that allows one to proceed with effective perturbative expansion of the prepo- 
tential is "the residue formula" , the variation of the period matrix, i.e. the third derivatives of J^{a) 
(see, e.g. [!§): 



datdajduk 



dak 



resd5=o- 



dujidujjdujk 
~SdS ' 



(5.1) 



where SdS = d (^^^ d^i or, explicitly 



SdS = dXd^ for Ad models 

SdS = R^^dlogXdS, for 5d models 



(5.2) 



SdS = dCd^ 



for 6d models 



^^Within the Whitham treatment, one can re-sort the argument identifying superpotentials and Ha milto ni ans of 
integrable systems, at least, in the SU{2) case. Indeed, one just needs to note that g/^i^^yi ~ 



^ (see (|6^)) 



is proportional to the gluino condensate (as a consequence of low-energy theorems). On the other hand, this condensate 
is obtained by taking the derivative of the superpotential w.r.t. A as it follows from Ward identities. 



Here the differential given on the bare torus degenerates into dw/w when torus degenerates into 
the punctured sphere (Toda and spin chains). We remark that although does not have zeroes on 
the bare spectral curve when it is a torus or doubly punctured sphere, it does in general however 
possess them on the covering C. Hereafter, we denote via ^ the variable living on the bare spectral 
curve (torus), while C is the variable living on the compactification torus (in 6 dimensions) or on the 
compactification circle (in 5 dimensions). 

In order to construct the perturbative prepotentials, one merely can note that, at the leading 
order, the Riemann surface (spectral curve of the integrable system) becomes rational. Therefore, the 
residue formula allows one to obtain immediately the third derivatives of the prepotential as simple 
residues on the sphere. In this way, one can check that the prepotential has actually the form ( ^j| ) 
We discuss these calculations later, but now see what the perturbative prepotentials are expected 
from the gauge theory point of view. 

As a concrete example, let us consider the SU{Nc) gauge group. Then, say, the perturbative 
prepotential for the pure gauge theory acquires the form 

^pert = ^ X! ~ "j)^ (oi - aj) (5.3) 

ij 

This formula establishes that when v.e.v.'s of the scalar fields in the gauge supermultiplet are non- 
vanishing (perturbatively are eigenvalues of the vacuum expectation matrix (</))), the fields in the 
gauge multiplet acquire masses m^r' = — a^' (the pair of indices (r, r') label a field in the adjoint 
representation of SU (Nc)). The eigenvalues are subject to the condition J2i = 0. Analogous formula 
for the adjoint matter contribution to the prepotential is 

while the contribution of one fundamental matter hypermultiplet reads as 

^pert ^ log + (5.5) 

i 

Similar formulas can be obtained for the other groups. The eigenvalues of {<j>) in the first funda- 
mental representation of the classical series of the Lie groups are 

Bn {SO{2n+ 1)) : {oi, a„, 0, -oi, -a„}; 

Cn iSp{n)) : {ai,...,a„,-ai,...,-a„}; (5.6) 

Dn (SO{2n)) : {ai, a„, -ai, -a„} 
while the eigenvalues in the adjoint representation have the form 



(5.4) 



Bn ■ {io^jS ± flfe}; j < k < n 

Cn ■ {±2aj;±aj ± Ofc}; j < k <n (5-7) 
Dn : {±aj ± flfc}, j < k <n 
Analogous formulas can be written for the exceptional groups too. The prepotential in the pure gauge 



theory can be read off from the formula (5.7) and has the form 



Bn ■ Tpert = ^ {{ai - ajf log (a^ - aj) + (a^ + ajflog {at + a^)^ ^ X] 

C„ : Tpert = ^ {{ai - ajf log (a.^ - + (a^ + a-,)^ log (a^ + a-,)^ + 2 ^ a^loga^; (58 

i,'L i 
Dn ■ ^pert = ^ ^ ((Oi - a^f log (ci - Gj) + (Oj + flj)^ log (fl^ + flj)^ 

The perturbative prepotentials in Ad theories with fundamental matter are discussed in detail in l3£ 



(see also |150[ ). 

All these formulas can be almost immediately extended to the "relativistic" 5d JV = 2 SUSY gauge 
models with one compactificd dimension. One can understand the reason for this "relativization" in 
the following way. Considering four plus one compact dimensional theory one should take into account 
the contribution of all Kaluza-Klein modes to each 4-dimensional field. Roughly speaking it leads to 
the 1-loop contributions to the effective charge of the form (oy = — aj) 



^ X! masses ^ log ( atj + ^ ) ~ log J| (Rsaij + m) ^ logsinhi?5aij (5.9) 



i?5 

iiidcioi^o fit fit 

i.e. coming from 4c? to 5c? one should make a substitution atj — > sinhi?50ij, at least, in the formula 
for perturbative prepotentials. 

The same general argument can be equally applied to the 6d case, or to the theory with two extra 
compactificd dimensions, of radii i?5 and i?g. Indeed, the account of the Kaluza-Klcin modes allows 
one to predict the perturbative form of charges in the 6c? case as well. Namely, one should expect 
them to have the form 

m ^ n 
i?5 i?6 



Tij ~ ^ log masses ~ ^ lo, 

masses m,n 

~ log Jl (^Rsa.j +m + logS* (^R^a.j 



(5.10) 



i.e. coming from 4c? (5c?) to 6c? one should replace the rational (trigonometric) expressions by the 
elliptic ones, at least, in the formulas for perturbative prepotentials. 

5.2 Perturbative prepotentials in integrable systems 

Now we are going to compare these results with calculations done for integrable models. We start 
with the integrable chains 

As we saw in the previous section (ss. 3. 2, 3. 4, 3. 5, 3. 6, 3. 8), the spectral curve in all cases of spin 
chains (fundamental matter) can be put in the form 

1^ + ^^=2P('^)(C) (5.11) 
w 

or 

W+^^^^^ " (512) 



ran 



In the perturbative limit, only the first term in the r.h.s. of these formulas survives so that the spectral 
curve becomes rational. 

The generating differential dS is always of the form (to unify the notations, we use the same letter 
( for A in Ad and ^logA in 5d) 

dS = CdlogW (5.13) 
The concrete forms of the functions introduced here are: 



a Q a ''^^ 

P(^)-n(C-«0, P(^)^n-^MC-a.), P(6)^n||_^ (5.15) 

i i i ^ ^ 

(in p(^^(C), there is also some exponential of C unless Nf — see ( 3.69| )). 

In order to come to the pure gauge theory in M and 5d, one should just put (^'•''^(C) equal to unity 
(it corresponds to sending hypermultiplet masses to infinity), while in 6d it is still unclear how to 
remove the restriction Nf — 2Nc which prevents constructing the pure gauge theory. 

The perturbative part of the prepotential can be calculated using these manifest expressions and 
the residue formula. 



5.3 Ad theory: a simple example of calculation 



Let us start from the simplest example of the Toda chain. In the perturbative limit (- 



second term in (5.11) vanishes, the curve acquires the form 



(0) 



0) the 



W^2PnA^), PAr^(A) =[](A-A,) 



(a rational curve with punctures PI) and the generating differential (3.19) turns into 



(5.16) 



d4tlt = AdlogPA,^(A) 



(5.17) 



Now the set of the Nc — 1 independent canonical holomorphic differentials is 



and the ^-periods = A; (the independent ones are, say, with i = 1, . . . , Nc — 1) coincide with the 
roots of polynomial P7v^(A), A^. The prepotential can be now computed via the residue formula 

dw.dujjdujk 

J^pert.ijk = res — — - — .. . ,. (5.19) 

■ •' diogPjv,=o dlogPjv (A)dA 



'^'^ These punctures emerge as a degeneration of the handles of the hypereUiptic surface so that the A-cycles encircle 
the punctures. 
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which can be exphcitly done. The only technical trick is that it is easier to compute the residues at 
the poles of dw's instead of the zeroes of dlogP. This can be done immediately since there are no 
contributions from the infinity A = oo. The final results have the following form 



1 6 1 



'ZitOF pertAii / , \ ' \ ' / ^ \ 



n T- 3 2 X - 1 AjTV 

2mTpert,Uj = T + T + 2^ T T : « ^ J> (5.20) 



giving rise to the prepotential formula 



J' pert — n . / J H-^ij) 

ir (5.21) 

/(4)(a;) =a;2loga;2 

Now let us turn to the case oi Ad N = 2 QCD with massive hypermuhiplets. Then, Q^'^\C) + 1 
and there arise additional differentials corresponding to the derivatives of dS with respect to masses. 
They are of the form 

= (5.22) 
It is again straightforward to use the residue formula^ 



dLJidujjduJK fT T ^ r- • ; \ a i 

diog-^=o dlog Y > dX 



and obtain for the prepotentiaQ 

i.j i,a a,j3 

upon the identification Ai = 



(5.24) 



5.4 Prom 4(i to higher dimensions 

Now let us turn to the bd case, to the XXZ chain. In the perturbative limit, the hyperelliptic curve 
( ^.1^ ) turns into the rational curve 



W = 2 = 2 



^^^^ \/n^'(A-A„) 



(5.25) 



^^Let us note that, despite the differentials du>a have the pole at infinity, this does not contribute into the residue 
formula because of the quadratic pole of d\ in the denominator. 

^^The term in the prepotential that depends only on masses is not essential for the standard Seiberg-Witten theory but 
is crucial for the prepotential to enjoy its main properties, similar to the WDVV equations. This term is unambiguously 
restored from the residue formula. 



with Xi — e^°*, J^i'^ '^i = 0, Act = e^™°. First, let us construct the set of holomorphic differentials 
dujj. 

The perturbative differential dS is 

dS = hogX dlog ( = hogXdX I y —V I (5.26) 



or, equivalently. 



dS ^CdCiy coth(C -a^)--y coth(C - m„) + TV, - -Nf (5.27) 

\ 1=1 a=l / 



The variation of ( |5.26| ) gives 



2^5* = -logAd ^ 



SX, 



. , X-X^ 2 

2=1 



1 SXg \ 

2 ^ A - A„ 1 

a— 1 / 



(5.28) 



^ cZA v""^ SXi f Xi \ 1 dA SXo 

^ X U- A» ~ A- AivJ ~ 2T^X V^-Ac 

i—l ^ ^ ' a=l ^ 

where it is implied that only A*", — 1 parameters A^ are independent so that — — X^il^V^ T^- Now 
the set of the N^ — I independent canonical holomorphic differentials is 



dujj 



dX 

T 



A - A, A - A 



AiAT^dA 



(A- A,)(A-AArJ 



i^l,...,Nc-l, A„ =A, -Aj (5.29) 
2ai, while the condition of vanishing the 



and one can check that the ^-periods of dS are logAi 
sum of all residues gives the relation ( 3.66| ) J^i — i S "^a- s.3.5, we discussed how this strange 
constraint should be removed off by using the twisted XXZ spin chain. The analytical explanation of 
this is that one can easily reproduce any lack of residues by putting some additional poles at infinities 
in C or at A = 0, oo. It is analogous to the procedure in Ad. We saw, however, that the same procedure 
does not work in d = 6. 



The additional differentials corresponding to the derivatives of dS (5.26) with respect to masses 
have the form R 



duJa 



dX An 



(5.32) 



^^It can be instructive to describe in detail why ( |5.1 
in the bd case, one should consider differentials on annulus 
. , one rather needs to take 



2A A — Xa 

arises not only in 4d, but also in 5d models. Strictly speaking, 
on sphere with 2 marked points. Therefore, instead of 



E 

m— — oo 



d( 



coth — ai)d( ■■ 



Taking now, instead of (5.18), the differentials 



\ + \id\ 
A - Ai 2A 



X + Xi dX A + Ajv^ dX _ 
A- Ai 2A " A - Ajv, 2A ~ (A - Ai){A - AjvJ 



(5.30) 



(5.31) 



we obtain again the formula (5.18). On the other hand, the mass differentials (5.32) look different from (5.3C). It turns 
out, however, that the difference does not contribute into the results and, therefore, one is free to choose any of these 
two possible mass differentials. 



The better choice for these differentials (see footnote |3^) is, however, 



duja 



dXX + Xo 



coth(^ - ma) 



(5.33) 



4A A - Aa 

Now one can insert these manifest expressions for differentials into the residue formula in order to 
calculate prepotential. 



Similar formulas are also correct for the XYZ chain. A generic XYZ curve (3.127) in the pertur- 
bative limit turns into 



W 



where ho ~ -^fs=. Then, the rescaled generating differential (3.128) is 
V Qo 

dS = (j2 log^* (C - «.) - ^ E log^* - "^-)) + ^°S^odC ' 



(dC I XI 3(C - a^) - 2 XI 3(C - ^a) 

i=l a=l 



(5.34) 



(5.35) 



where the linear piece disappears due to ^ = ^ ^ and 3 denotes the Weierstrass ^-function 



[1091. This expression should be compared with (5.27) in the 5c? case 



Note that the sum of residues vanishes only provided the constraint ( 3.12(: ) is fulfilled. It can not 
be removed by adding an additional pole at infinity, since the variable C, in contrast to the 4d and 5d 
cases, lives on the compact surface with out boundaries (torus) and, therefore, is always finite. 

Now one can take the total variation of dS (we neglect the trivial variation SHq giving rise to the 
holomorphic differential on the bare torus) 



SdS""^^ - X 3(C - ai)5a^dC - 2 £ ^(C " m„)(5m„dC (5.36) 

i—l a—1 

where 3(C|''') = 3(CI''') ~ v{''')C- Therefore, the differentials related to gauge moduli are 



d^t = idiC - cii) - 3(C - oatJ) dC 



while those related to masses are 



(5.37) 



duja = "2 (3(C - rua) - 3(C ~ awj) 



(5.38) 



The second term in this expression is due to the condition ( 3.126| ). Thus, one can see that, in the 6d 
case, mass moduli are practically identical to the gauge ones. 



5.5 Perturbative prepotentials: fundamental matter 

The result of calculation for the prepotentials using all these manifest expressions for the differentials 
and the residue formula is always of the form 



CI 



87ri 



a,/3 



247ri 



(27V, -TV,) (E«' + iE-') 



The explicit form of these functions is 



so that 



= »'iog»-, = (»■ + :^) = 5 
/™(.) = E/™(. + ^ + ^)=i:/'=' 

^Lis^^ ^e^^'"^'^ + quadratic terms j 



,a;^|--Li3(e-2|-l 



1 , 



^ ^5 

-R5 -Re / V -Re 



(5.39) 



(5.40) 



= fog:r, /(^)"(a;) - logsinhx, /(^^"(a:) = log0*(x) 



(5.41) 



Note that, in the 6d case, Nf is always equal to 2N , and, in d = 5, 6, there is a restriction ^ = 
S Ci = 5 "^a which implies that the gauge moduli would be rather associated with shifted by 
the constant j^X)™"- these formulas, Li3(a;) is the tri-logarithm, while Li3^q(a;) is the elliptic 



tri-logarithm [151|. In order to get the prepotential for the pure gauge theories, one just need to omit 



all the terms containing mass moduli ttiq. 

These formulas perfectly match the perturbative prepotentials obtained in s.5.1 in supersymmetric 
gauge theories. This gives another support for their identification with integrable systems. 

Let us discuss now the cubic (in moduli) terms appearing in the perturbative prepotential only 
in higher-dimensional theories. First, take the case, where these terms correspond to the Chern- 
Simons term in the field theory Lagrangian Tr(A A F A F). In our computation they appear due to 



the constant piece {Nc — ^Nf) in the brackets in (5.27). 

In fact, let us choose the special values of the second Casimir functions at m sites so that the 



polynomial Q(A) acquires the form A'"(3^^(A), i.e. m of 2Nc factors in (3.63) turn into A™. In this 
case one gets the constant piece Nc-m- Nf/2 in dS^^\ with an integer m < 2Nc ~ Nf. Therefore, 



the coefficient in front of cubic term for d = 5 in ( 5.39 ) can be made equal to 



6 



{Nc - m - Nf/2) 



(5.42) 



Restoring the dependence on i?5 in d = 5 in ( ^.39| ) (ai aiR^ and ma vriaRb) one can study 
the different limits of the system. The simplest limit corresponds to the 4(i case and is given by 
i?5 0. In this limit, f^^\x) -~> f^^^{x), the cubic terms vanish and we reproduce the perturbative 



4c? prepotential ( 5.24| ). At the level of integrable system it is enough to replace Si R^Si 



in the Lax operators ( 3.57| ), ( |3.67 ) in order to reproduce the Lax operators of the XXX chain. 



Another interesting limit is the limit of flat 5d space-time, i.e. -R5 00. In this limit, only cubic 



terms survive in the prepotential (5.39) for d = 5 (one should carefully fix the branch of /'^■'(x) which 



leading to appearing of the absolute value in ( p.43| )) 




(5.43) 



The first one comes from tlie function 



In fact, there are two different sources of cubic terms [152 

Since this function can be obtained as the sum over the 4(i Kaluza-Klein perturbative contri- 
butions, these cubic terms have perturbative origin and come from the 1-loop (due to the well-known 
effect of generation of the CS terms in odd-dimensional gauge theories) . 



The second source of the cubic terms is due to the bare CS Lagrangian. As it was shown in [152|, 
one can consider these terms with some coefficient CcV- 



Ccl 

6 



(5.44) 



restricted only to satisfy the quantization condition Cc/ + ^ G Z and the inequality \cci\ > Nc — 
^f. In ( ^.42| ), we easily reproduce these conditions so that (5.43) coincides with |152|, provided the 
prepotential is defined in a fixed Weyl chamber. We leave, however, the wall crossing jumps out of 
the discussion. 

In the 6(i case the calculations are quite similar, and one would expect the similar cubic terms 
correspond to Tr(i^ f\ F f\ F). However, although all the difference in calculations is due to slightly 
different mass differentials ( [5.33 ) and (5.38), the bare cubic terms in 6c? are absent. This occurs due 
to the Nf ~ 2Nc condition (since the coefficient in front of cubic terms is proportional to Nc — ■^Nf) 
and, therefore, due to the cancellation of the constant term in the brackets of ( ^.35 ) as compared with 
( p^ ). This is one of the problems with 6d theories: one can not deform (twist) the theory so that 
the constraints Nf = 2Nc and ( |3.126| ) look unavoidable. As for the "quantum generated" terms, i.e. 
those coming from the function /^^^ they are certainly presented, with the coefficient equal to the 5d 
coefficient. 



5.6 SU{2) adjoint matter in 5d 

Now let us consider the adjoint matter case. We discuss here only the 5d theory, i.e. the Ruijsenaars- 
Schneider system, since the limit to the M case is done straightforwardly, while many concrete calcu- 
lations in the Calogero-Moser system look even more tedious. 

Suppose we parameterize a bare curve that is elliptic by the algebraic equation, 

- ax^ - f3x -J. (5.45) 

Then we can use the standard parameterization via Weierstrass p-function and its first derivative: 
X = p(n + a/3 and y = "hp'iO so that = -24?. 

^ - y 

The perturbative (weak-coupling) limit is then given by the imaginary period 2uj' becoming infinite, 
T —> ioo. In this weak coupling limit (without any double-scaling which we did considering the 
degeneration to the periodic Toda chain) the bare curve ( 5.45| ) becomes 



y^=x^{x-l). (5.46) 



CO 



Then J ^ X - 1/3 with x = V V ^ a = 1, /3 = 7 = and -> i.e. 

dS = R^hogXdS, R^hog\^^j£=. At the same Umit, 

, 1 sinfsine , sin£sine T^pf/Q mf*/'? 

Vsin"^ e - sin"^ e ^fiin{€ - sm(e + sin(e ± 

(5.47) 

The corresponding integrable system is then nothing but the trigonometric Ruijsenaars-Schneider 
model we aheady met in s.2.2. Finally the special point e = w = § sees 

piu) = ei ^ ^ (5.48) 



and the particular combination (3.90) becomes 



co(^|e = 7r/2) ^ tanC. (5.49) 
Now let us look at the period matrix in the perturbative limit in the simplest SU{2) case, the 



corresponding spectral curve being ( 3.95 ). As usual, the derivatives of the 1-form dS with respect to 
the moduli are holomorphic differentials. A particular choice of coordinates for the moduli will lead 
to the canonical holomorphic 1-differentials. For the example at hand there is only one modulus and 

If we order the roots 63 < 62 < ei, we may take for the ^-integral 



d. = 1 r dv = ' ^ / + - - 



*7ri?5V^(f + P(e)- 62) (61 - 63 ) VV (x + - e2)(ei - 63) 



(5.51) 



Here K{q) is the complete elliptic integral of the first kind. Dividing dv by the right hand side of this 
expression would then give us the canonical holomorphic 1-differential dQ. Similarly for the B-integral 
we have 



" " iR5^(f + p(e)-e,)(e, -63) VV (t + " e^X'i " es) 

and so the period matrix is 



(5.52) 



One can take the limit e ^ in this expression to obtain T = — t. This agrees perfectly with formula 
^ciass ^ TO? for the classical part of the prepotential and justifies our identification of the gauge 
theory coupling constant and the modulus of the bare spectral curve. 

Let us now calculate the same quantities in the perturbative limit: 



with [/^ — ^.^j ^ + We have 62 = 63 = — 1/3 and ei =2/3 in this hmit and so the A-period in this 
case shrinks to a contour around x — 0. Now dv — \ / {2iR^U) and we may identify the canonical 
differential dil ~ 2iR5Udv. This result also follows from ( |3.5l| ) upon using 

2 

K(q) = -(1 + ^ + ...) 

The B-period of dfl again gives the period matrix. The corresponding integral now goes (twice) 
between x = and x = 1, and the integral dv diverges logarithmically in the vicinity of a: = 0. 
This divergence was to be expected because the period matrix contains a term r on the right hand 
side (coming from the classical part of the prepotential J-'dass = —ra'^) and the perturbative limit is 
given by r — !■ ioo. Upon making the rational substitution x — j^— we obtain 



U dx ,. U dv 



Tpert — <t> dn — — / , . = lim 



TT Jo X^{x~l)iU^-x) e^^lirVlP^Je 



(5.55) 



-— Imi (log- + —log- — 

ITT e^O V 4/ ITT 1 — U'^ 



where e is a small- a; cut-off. Thus, the U dependent part of this integral is finite and can be considered 
as the "true" perturbative correction, while the divergent part just renormalizes the bare "classical" 
coupling constant r. Again the same result follows from ( [3.5S ) upon using 



hm (K{q)-\H-^\ 
The final ingredient we wish are the a-variables, i.e. the A-period of dS itself. This will correspond 



to the integral of ( 5.51 ) with respect to u, which is a rather complicated integral. For our purposes 
the perturbative limit will suffice when there are several simplifications. From the definition of ( |1.3| ) 
and that of dS we find 

_ 1 / logAdx _ logA|j;=o 



i?5 / x^Jx - 1 ii?5 

Now at X = we have that cq = i sine while at the same time A + A^^ = cqu. Together these yield 

2 cos ai?5 — iu sin e 

and 



f/^-^ + -V = -V(l- cos^ «^5) - '^^^ (5-56) 
4 sin e sin e sin e 



Substituting this expression for U into (5.55) we obtain for the finite part of the period matrix 



= -log = llog— '^^^ -. (5.57) 

^ ZTT sin e — sin aR^ in sm (e + an^ ) sm (e — aU^) 



5.7 Perturbative prepotentials: adjoint matter 



It is quite clear that to repeat this calculation for more general curve is a hard job. Instead, let us see 
what is happening with the full spectral curves under the perturbative degeneration^. The results 
turns out to be so strikingly simple, despite quite involved form of the full spectral curve that this 
allows us to immediately calculate the perturbative prepotential for arbitrary SU{Nc). 



We again start with a warm-up example of SU{2), with the curve (3.95) at the special point 
e = 7r/2. Then, upon using ( 5.49 ) we see that ( 3.95 ) reduces in the perturbative limit to 



cot^ 



uX 



A2 + 1 

The rational spectral curve for this situation can be put in the form 



(5.58) 



(A + Ai)(A + Ar^)' 



(5.59) 



upon setting w ~ — e^'^ and iu — (Ai + X-^ ^). 



Similarly in the SU{3) case, in perturbative limit and at the special point e — 7r/2 (3.102) turns 



into 



X^ -vX = i{l-uX^)tan^ (5.60) 

for the choice c — —icQ. Equivalently, 



X{X^+v) 

I uau(; = 

This curve can be put in the form 



^tan^=2y^l^. (5.61) 



(A-AiKA-A^KA-Aa) 

(A + Ai)(A + A2)(A + A3)- ^^-^'^ 
upon the identification w = — e^'^ and with u = Ai + A2 + A3 and v — A1A2 + A2A3 + A3A1 = 

A1A2A3 (Ar' + A^i + A3 1), A1A2A3 = 1. 

Now let us see what is the curve at the special point where e = 7r/2 for the general SU{Nc) case 
ff22l. 



We saw already that substantial simplifications occur at the special point e = uj before even 

considering the perturbative limit. Now in the perturbative limit cr(z) sin(z)e^ Using this we 
easily obtain that 

P(A)-H^^^P(-^ ^^_^2i,_(,N.Pi^ .5 63^ 

*^°t^-p(A) + (-)^cP(_A)' w- e -{ ) (5.63) 

with dS*- logA^. 

We finally turn to the case of generic e in the perturbative limit. The system in this limit is 
described by the trigonometric Ruijsenaars- Schneider model. We have already seen that in this limit 
fT-functions are proportional to sines with exponential factors. However, the periodicity of our spectral 



Similar degeneration for the Calogero curves was studied in |153 



curve means that these exponential factors must cancel amongst themselves. For example we find that 
in this limit 



1 sin(e+(fc-l)6) (-1) 

Dkim = =F 



sin($-e) ■sin'=('=-2) 
where c_ is given by (^.47]). Now the ratio 



sin(C + (k - f )e) 



cos fee + sin ke cot(^ — e) 



(5.64) 



(5.65) 



sin(C - e) 

is expressible in terms of the single function cot(^ — e). This simple observation enables us to separate 



the variables ^ and A in the equation for the spectral curve upon choosinj 
hk = Hk/siVL^'^^~'^\t) we may simplify the spectral curve to give 



-ie^'-c-. With 



i cot(^ 



Jeven 

where P(A) = J2aU k hki—^)^"^^ = Y{^''{^ — e^'*^*) with some constants ai, 



P(Ae-2-) 



P(Ae- 



(5.66) 



the variable w 



=2i(?-c) 



, one finally arrives the spectral curve in the form 



0. Introducing 



w 



P{X) 



with 



dS = logA 



P(Ae- 
dw 



(5.67) 



(5.68) 



Thus, we have shown that our system leads in the perturbative limit to the rational spectral curve 

One may now calculate the corresponding prepotential 



( p.67| ) and the generating differential (5 



using the residue formula. In fact, this calculation has been already done for the spin chains in the 
previous subsections, where the residue formula is applied to a general rational function of the form 



w 



nf(A 



A.) 



with A, 



' Si ~ 0, Aq 



nr(A-A„) 

By choosing Nj — 2Nc, with the hypermultiplets masses 



pairwise coinciding and equal to + e we obtain our curve (5.67). One then finds the prepotential 



from (5.39) gives the prepotential for the Ruijsenaars-Schneider system 



^p^rt = ^E/^'ni?5a.,) - ^E/^'^(^5a., +e) 



(5.69) 

and we are done. Further, upon setting e — 7r/2 we reproduce the results above. 

Now in order to obtain the prepotential in the Ad theory with adjoint matter described by the 
Calogero-Moser system, one can just take the limit we described above, P5 — > and e = MP5. 
Therefore, there is no special point e = ^ in the Calogero-Moser system! The perturbative prepotential 
for the Calogero-Moser system |39, 154, 39) is obtain by degenerating ( f3.69| ) so that f^^^a) — > f^'^\a): 



Note that this separation is not in general possible non-perturbatively where there is no analogue of eg. ( |5.65| ). 
(This is because the ratios "^^^^^^^^^^^"^ transform differently under § — e— +5 — e + 1 and ^ — e^^ — e + r for different 
k.) The two notable exceptions are: 



(i) The case of SU{2), when the variables separated for the standard choice c = cq, see eq.(3.95) 



(ii) The case of e = it/2 (for any SU{Nc)), when separation of variables survives non-perturbatively, see eq.(3.106) 



Sttz^" Sm^' ' - (5-70) 

which can be also calculated from the spectral curve obtained by the degeneration of ( ^.67 ) 



with the same differential ( 5.68| ) 



In section 7, we shall see that this form of the perturbative prepotential persists in the Dell systems. 

6 Theory of prepotential 

So far we discussed only the properties of Seiberg-Witten theory related to underlying integrable 
systems. However, the whole construction turns out to give far more. In fact, in previous section 
we just discussed how to construct different physical data in terms of integrable systems, but almost 
missed meaning of several key notions. Maybe the main object that is to be discussed in far more 
details is the prepotential, the quantity not properly interpreted yet in previous sections. 

Indeed, looking from the integrable standing point, one should find natural both emerging the 
Riemann surfaces with their moduli space as the corresponding spectral curves, as well as the proper 



coordinates on the moduli space given by the action integrals (1.3). However, the notion of the 



prepotential which allows one to make quantitative predictions and defines the low-energy effective 



action was so far introduced by hands, just via formula (1.3). This formula implies that not only A- 
but also i?-periods of the generating differential dS should get a proper interpretation. This means 
that one should better understand meaning of dS and its properties. One has particularly to recover 
the reason for dS to give rise to holomorphic differentials, or, to put it differently, the reason for the 
second derivative of the prepotential to be a period matrix. 

In this section we are going to recover another kind of integrability, Whitham integrability which 
stands behind the notion of the prepotential and generating differential dS and discuss several very 
essential its features, in particular, the WDVV equations ( |1.8| ). 

6.1 Prom Seiberg-Witten theory to Whitham hierarchies: RG dynamics 

Let us consider a generic recipe how the generating differential is introduced in integrable systems 
[H H, 155 1 . It is essentially based on the analogue of the renormalization group flows in field theory, 



where one performs the integration over heavy modes giving rise to the effective action for light 
modes. Within the framework of integrable many-body systems, an analogous procedure looks as 
follows. First, one should start with a solution to the equations of motion with moduli. At the next 
step, one assumes that the moduli adiabatically depend on some new time variables. In this way, one 
obtains the Whitham dynamics on the moduli space, with the generating differential dS playing the 
role of action differential and the prepotential being logarithm of the r-function. For instance, for the 
finite-gap solutions that are in charge of supersymmetric gauge theories without adjoint matter, i.e. 
those described by hyperelliptic curves, it is natural to choose the corresponding branching points as 
these slowly evolving moduli. 

Therefore, we come to the conclusion |Q that Seiberg-Witten theory has much to do with two 
different structures: with finite-gap solutions of integrable systems of the Hitchin type or of spin chains 



(so that the integrals of motion are fixed and related to the vacuum expectation values (vev's) of scalar 
fields in supersymmetric theories) and with the corresponding Whitham hierarchy constructed in the 
vicinity of these solutions and describing the evolution of integrals of motion. 

Since the vev's of scalar fields physically depend on the scale Aqcd, one associates one of the 
Whitham evolutions with the RG dynamics w.r.t. this scale (in fact, Aqcd is identified with the 
first Whitham time). Other slow (Whitham) times describes vev-dependence on perturbations of 
the original SYM action by different operators (from topological sector). To put it differently, while 
Seiberg-Witten theory describes the low-energy effective action, the Whitham dynamics governs the 
behaviour of correlators at low energies. This correlators are from the topological subsector of the 
theory, since only such correlators survive at low energies. Note that this subsector is manifestly 



described by the Donaldson theory, |f56, f57, 158|. 

In fact, this is the prepotential that is the generating function of these (connected) topological 
correlators, i.e. logarithm of the corresponding partition function. Therefore, now the prepotential 
depends on both the integrals of motion and the Whitham times as on independent variables 
Before defining such an object, let us briefly review the very framework of Whitham hierarchies, from 
the very general point of view. 

Following l.Krichever 



f55|, we consider initially a local system of holomorphic functions Qj on 



complex curve, i.e. functions of some local parameter ^ in a neighborhood of a point P. One can 
introduce a set of parameters tj and define on the space {^,ti) a 1-form 



= ^l]/(P,t)<5i7 



(6.1) 



Considering its total external derivative Sco = J^i ^^i{Pj t) A Stj, where Suj = d^rtid(,/\Sti + djil.i6tj A 
St J, one can define general Whitham equations as 



A (5w = 



(6.2) 



so that it is necessary to check the independent vanishing of two different terms: St"^ and St^6^. The 
second term gives 



where [...] 
ordinates 



Y,di^lidjnc]^0 (6.3) 
means antisymmetrization. Fixing now some I = Iq and introducing the "Darboux" co- 



one can get from (6.3) the Whitham equations in their standard form |155| 



(6.4) 



dinj-djni + {ni,nj} = o 



(6.5) 



dx dp dx dp 

The explicit form of the equations (|6.5D strongly depends on choice of the local co-ordinate p. Equa- 



tions (6.2), (py^) and (5.5) are defined only locally and have a huge amomit of solutions. 



Now let us note that the Whitham equations (6^) can be considered as compatibility conditions 
of the system 



(6.6) 



The function A can be used itself in order to define a new local parameter. In this co-ordinate, the 



system (3.5) turns into [155 



dni{X,t) 



dtj 



dnj{\,t) 



X— const 



dti 



X— const 



or for the differentials dfi/ = dy^^li ~ ^^^d\ 

ddQi{X,t) 



ddnj{X,t) 



X— const 



dti 



X— const 



The Whitham equations in the form (E^) imply the existence of a potential 5* such that 



dS{\t) 



dti 



X— const 



With this potential, the 1-form uj can be rewritten as 



5S 



dS{X,t) 
dX 



dX 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



To fill all these formulas with a real content, one can consider interesting and important examples of 
solutions to the Whitham hierarchy arising when one takes the basis functions il/ to be the "quasi- 
energies" of the finite-gap solution; these are the (globally multi-valued) functions whose periods 
determine the "phase" of the quasi-periodic Baker- Akhiezer function corresponding to this finite-gap 
solution. The potential should be then identified with 



dS 



dS_ 
dTi 



ddS 



(6.11) 



where dS is, as usual, the generating differential for the finite gap system. Usually for KP/Toda 
hierarchy these functions are taken to be "half" -multi- valued, i.e. their A-periods are fixed to be zero: 
dVlj — 0. The equations (6.11) are obviously satisfied in a trivial way in the case of the (finite-gap 



solutions to) KP/Toda theory, when {ti} are taken to be external parameters and {51/} do not depend 
themselves on {tj}. The idea, however, is to deform this trivial solution to (|6^), ( |6.3D and ( |6.5| ) into 



a nontrivial one by a flow in the moduli space. Practically it means that the formulas (3.11) should 
be preserved even when {ri/} depend on {t/}, but in a special way - determined by the equations of 
Whitham hierarchy. In other words, the Whitham equations correspond to the dynamics on moduli 
space of complex curves and the moduli become depending on the Whitham times. 

To realize all this technically, one choose d^i to be a set of (normalized) differentials on the 
complex curve holomorphic outside the marked points where they have some fixed behavior. These 
differentials depend on Whitham times only through moduli dependence. 

In fact, for this whole construction to work (i.e. for the Whitham hierarchy to have some non-trivial 
solutions), one needs to choose both differentials (with their normalization) and local parameter in a 
clever way. This is the thing to be added to Seiberg-Witten theory. Indeed, Seiberg-Witten theory is 
given by the set of data that includes Riemann surfaces, their moduli space and the differential dS, 



while the Whitham set-up additionally includes a (non-invariant geometrically) proper choice of the 
local parameter. 

An important example of the A co-ordinate is the hyperelliptic co-ordinate on the hyperelliptic 
curve (corresponding to the KdV equation) 

2g+2 

y'^RiX)^ YliX-rp) (6.12) 

Properly choosing the differentials dflj and expanding the equations (|6.8[ ) near the point A = one 
gets the Whitham equations over the finite-gap solutions of the KdV equation 



dva dQ J 



dtj dVLj 

The set of ramification points {r} — {r^} is usually called Riemann invariants for the Whitham 
equations. 

The particular case of Whitham hierarchy associated with the finite-gap solution of the KdV 
hierarchy described by the hyperelliptic curve 



n 



y^ = R{X)=\[{\-rp) (6.14) 



was considered first in [ 159 



Note also that in the framework of Whitham hierarchy it becomes possible to introduce a coun- 
terpart of the tau-function, again a generating function for the solutions, which can be symbolically 



defined as |155 | 



logrwhitham= / dShdS (6.15) 



s 



where the two-dimensional integral over E is actually "localized" on the non-analyticities of S. In 
fact, this log7\vhitham IS nothing but the prepotential. 

Recently an analogue of Whitham dynamics for deformations of the analytic curves instead of the 



complex surfaces has been found |16C, 161]. At the language of dynamical systems it means that the 



dynamics in the space of real Hamiltonians instead of holomorphic ones is considered. 
6.2 Whitham hierarchy for the Toda chain 

Now let us discuss how one can choose the differentials dfJ/ in order to get non-trivial solutions to 
Whitham hierarchy. They can be defined for the following set of data: 

• complex curve (Riemann surface) of genus g\ 

• a set of punctures Pi (marked points) ; 

• co-ordinates in the vicinities of the punctures Pi] 

• pair of differentials, say, (dA, dz) with fixed periods. 



►71 



We start with considering the simplest case of a single puncture, say at ^ = = 0. This situation is 
typical, e.g., for the Whitham hierarchy in the vicinity of the finite-gap solutions to the KP or KdV 



hierarchies. For instance, in the KdV case the spectral curve has the form (6.14) so that the marked 



point is at A = cx) but, in contrast to the Toda curve (3.8), there is only one infinite point (since it is 
the ramification point). 

Given this set of data, one can introduce meromorphic differentials with the poles only at some 
point Pq such that in some local co-ordinate ^: ^(Pq) = = 

^ = (r""' +Oil))d^, 71 > 1 (6.16) 

This condition defines dfln up to arbitrary linear combination of g holomorphic differentials dui, 
i = 1, . . . ,g and there are two different natural ways to fix this ambiguity. The first way already 
mentioned above is to require that dD,„ have vanishing A-periods, 

dn„ = Vi,n (6.17) 

A, 

Their generating functional {( = ^{P'), £, = ^{P)) 

oo 

W{^, - E nC-'dCdil,,{0 + ... (6.18) 

is well known in the theory of Riemann surfaces. It can be expressed through the Prime form E{P, P') 
(see Appendix A of |8|] for details): 

W{P, P') = dpdp'\ogE{P, P') (6.19) 
Such W{P, P') has a second order pole on diagonal P — > P', 

^) - + = E "^C"-'dC + Oil) (6.20) 



It is the differentials ( 6.16| ) that should be related with the potential differential dS by 



rfa. = ^ (6.21) 

dtn 

The second way to normalize differentials is to impose the condition 

— — = holomorphic (6.22) 
o moduli 

so that these differentials becomes similar to the generating differential dSsw of integrable system. 



For example, for genus 3 = 1 in the co-ordinate 5'~5 + 1~5 + t the formulas acquire the form 

1 ("-"l^+l 

dQi = {p(^) - const)d(, = --p'(Od(, ...,dnri = — d"-'^p{Od^ 

2 n! 

and 

E('_Nn + lc/n — 1 Qn-1 
(n — 1)! og""-'- 

n — l 

where * denotes the (on torus the only one) odd theta-characteristic. For g = 1 i^^iO — ^i{0)du}i is just d^. Let us 
also point out that chosen in this way co-ordinate t; is not convenient from the point of view of Whitham hierarchy since 
its "periods" (t = d^) depend on moduli of the curve. 



TO 



Now let us specialize this description to the simplest case of Seiberg-Witten theory corresponding 
to the 4d pure gauge theory, i.e. to the Toda chain. Therefore, we deal with the family of spectral 



curves (3^). The first problem is that naively there are no solutions to eg. ( 6.22 ) because the curves 
( ^^ ) are spectral curves of the Toda- chain hierarchy (not of KP/KdV type). The difference is that the 
adequate description in the Toda case requires two punctures instead of one. As already mentioned in 
sect. 3. 3, the curves ( |3.6| ) have two marked points and there exists a function w with the A'^-degree pole 
and zero at two (marked) points A = oo-|-, where ± labels two sheets of the hyperelliptic representation 



(3.6), w{X = 00+) = oo, w{X = cx)_) = 0. Accordingly, there are two families of the differentials c?ri„: 
dfl^ with the poles at cx3+ and dQ^ with the poles at cxd_. 

However, there are no differentials dfl^, only dfln — dfl^ + dfl~, i.e. condition ( 6.22| ) requires 



dO„ to have the poles at both punctures. Moreover, the coefhcients in front of w^^^^" at oo+ and 
yj-n/N ( 6.16| ) coincide (in Toda-hierarchy language, this is the Toda-c/iam case with the same 



dependence upon negative and positive times). 

The differentials dO„ for the family ( |3.6|) have the form fs^ : 

df^„ = i?„(A)^=P;/-(A)^ (6.23) 

w ^ w 

The polynomials Rn{X) of degree n in A are defined by the property that P'SKn — R'^SP is a polynomial 
of degree less than - 1. Thus, i?„(A) = P;/^^(A), where {^t^_^Ck\*'^ ^ Yl^oCk\^- For 



exa: 



mple0: 



i?o — 1, 

Ri = A, 

R, = X^-l-u,, 

3 3 (6-24) 



2{N, - 4) 



Ra — X —'U'2X^ ~^3^ ^ I ~ "2 I ) 



These differentials satisfy ( |6.22 ) provided the moduli-derivatives are taken at constant w (not A!). 



Thus, the formalism of the previous section is applicable for the local parameter ^ = w^^/^'^.Q 

Thus, among the data one needs for the definition of solution to the Whitham hierarchy associated 
with the Toda chain, there are the punctures at A = cxdj- and the relevant co-ordinates in the vicinities 
of these punctures ^ = w^^^-^" ^ A^^ at oo^ and f = w^^^^" ^ A^^ at cx)_. The parameterization in 



terms of w, however, does not allow one to use the advantages of hyperelliptic form (3.8) 



^■'In the case of the Toda chain, since there are two punctures, one can consider the meromorphic differential with 
the simple pole at each of them. 

•^^Rule of thumb generally is to choose as the local co-ordinate for the Whitham hierarchy the parameter living on 
the bare spectral curve. For the Toda case, one may choose as a bare curve the lu-cylinder with the corresponding 
generating differential dS = Xdlogw. 
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6.3 Solution to Whitham hierarchy and prepotential 



Now we are going to extend the system of differentials dQi introduced in the previous section to 
include the holomorphic differentials. This will allow us to construct the function that, as a function 
of one set of variables, is the prepotential of the Seiberg-Witten theory (as defined in sect. 2) and in the 
other set of independent variables is the (logarithm of) r-function Twhitham of the associated solution 
to the Whitham hierarchy. On this way, we also get a generic solution to the associated Whitham 
equations. 

Thus, we extend the set of Whitham equations by the equations involving the holomorphic differ- 
entials duJi and moduli a^: 

ddfln ddhji ddiOi ddujj 
dai dtn ' daj dai 

This system is solved by the differential dS that satisfies 

ddS , ddS ,^ 

— — = du)i, — — = dlln (6.26) 



Then, the first equation in (6.26) implies that dS is to be looked for as a linear combination of the 



differentials dQn, satisfying ( 6.22| ). Let us, following ^ Q introduce a generating functional for 



drin with infinitely many auxiliary parameters t„: 

g 

d5 = ^ tndftn = ^ aiduji + ^ tndfln (6.27) 

n>l i—1 n>l 

The periods 

a^= <^ dS (6.28) 
can be considered as particular co-ordinates on the moduli space. Note that these periods do not 



exactly coincide with (1.3), eq.(6.28|) defines ai as functions of hk and tm or, alternatively, hk as 



functions of ai and t„ so that derivatives dhk/dtn are non-trivial. In what follows we shall consider 
ai and t„, 

tn = res4=o CdSiO (6.29) 

as independent variables so that partial derivatives w.r.t. at are taken at constant tn and partial 
derivatives w.r.t. tn are taken at constant ai. 



The differential dS (3.27) determines a generic form of the solution associated to Seiberg-Witten 
type Whitham hierarchy. The Whitham dynamics itself for given dS can be formulated in terms of 
equations ( |6.26D . Note that, if one restricts himself to the Whitham hierarchy with several first times 
(generally, for the genus g complex curve there are g + n — 1 independent times, with n being the 
number of punctures), all higher times in ( |6.27| ) play the role of constants (parameters) of generic 
Whitham solution (see example of SU{2) below). 

Note that the Seiberg-Witten differential dSsw is dSsw = c?^^2, i-S. 



and a-variables are naturally associated with the Seiberg-Witten moduli, while i-variables - with 
the corresponding Whitham times (although aiS can be also considered as variables of Whitham 



dynamics, cf. with (6.25)) 



Now one can introduce the Whitham tau-function (6.15) whose logarithm is a prepotentialJ-{ai, t,; 
logTwhitham by an analog of conditions (1.3): 



dS, 



dt„ 



1 



27rm 



reso C"dS 



Their consistency follows from (6.26) and Riemann identities. In particular. 



(6.31) 



-reso r 



.ddS 



-reso C "dfl^ = 



-reso S.-^'dflr. 



(6.32) 



dtmdtn 2TTin dtm 27rm " ^ "' 27rim 

This calculation implies that the definition ( |6.3l| ) assumes the co-ordinate f is not changed under the 
variation of moduli. It means that this provides a moduli-independent parameterization of the entire 
family - like w in the case of ( |3.6D . Since moduli- independence of ^ should be also consistent with 
( |6.22 ) , the choice of ^ is strongly restricted: to 
previous section). 



w 



±1/Na 



in the case of (3.6) (see the discussion in the 



The last relation in (5.32) (symmetricity) is just an example of the Riemann relations and it is 
proved by the standard argument: 



0=1 dilm. A diln = J dfln 

= + + 



dfln — f dfln f dO,,- 



- — rres i^d^rnd ^<ir2„) 



1 



(6.33) 



-reso ^ "dil,; 



-reso f ™-d^n 



2iTin 2TTim 

where ( |6.17 ) and (3.16) are used at the final stage. The factors like arise since ^""^^c?^ = 
— (i(^^"/n). We shall also use a slightly different normalization dfln ^ ^"n/Nadw _ NjL^yrjn/N^ ^ accord- 



ingly the residues in ( 6.33 ) and ( 6.31 ) and ( 6.32 ) will be multiplied by Nc/n instead of 1/n. 

By definition, the prepotential is a homogeneous function of its arguments ai and f„ of degree 2, 



- 



' doti 



duida. 



+ 2aitr, 



daidtr, 



(6.34) 



Again, this condition can be proved with the help of Riemann identities, starting from (S.31), (6.2S) 
and ( |6.29| ). At the same time, is not just a quadratic function of ai and tn, a non-trivial dependence 
on these variables arise through the dependence of duji and dfln on moduli (like Uk or hk) which in 



their turn depend on Oi and t„. This dependence is obtained, for example, by substitution of (6.27) 
into (HH): 



dui Btn 



dflr, 



(6.35) 



I.e. 




8dil„ 
dui 



dflr, 



(6.36) 



The integral in the l.h.s. is expressed, according to (6.22), through the integrals of holomorphic 
1-differentials, while the integral in the r.h.s. - through the periods of dS. If 



^ = j( dVi = i:^i (6.37) 



then 



trap^^ E^/a. (6.38) 



The relations ( 3.38 ) can be thought of as the other form of the Whitham hierarchy. 

Thus, we have constructed the prepotential for the Whitham hierarchy associated with Seiberg- 
Witten theory. The Whitham equations themselves describe the dependence of effective coupling 
constants given by the period matrix of the Seiberg-Witten curve on bare coupling constants i„. The 
prepotential also depends on characteristics of the vacuum parametrized by vev's Uk- Using this 
interpretation, one expects |Q that the constructed Whitham prepotential has to be associated with 
the generating function of topological correlators, the Whitham times being bare coupling constants 
of different local operators. It issues the problem of calculating derivatives of the prepotential with 
their further comparison with the topological correlators. Indeed, the first two derivatives have been 
calculated in |^ and turned out to be 

1^ = 1^ mi™H™+i,„+i = Tj—hUn+i + Oit^M, ■ ■ ■) (6.39) 



m 

da^dt^ 2TTin da^ 



(6.40) 



- t;- ^m+i,n+i + — ^T- 7r—d,AogSE{0\T) (6.41) 



9t™at" 2TTi V mn da* daJ 

etc. In these formulas parameter (3 — 2Nc, m,n — I, . . . , Nc — 1 and 'Hm,n are certain homogeneous 
combinations of u^, defined in terms of Pn{X): 

n^+l.n+l = -— reSo, (p"/^=(A)dP™/^=(A)) = Hn+l,rn+l (6-42) 

and 

Hn+i = Hn+ia = -— reSooP"/'^iA)dA (6.43) 
n 



The Seiberg-Witten theory itself allows one to evaluate derivatives (6.39), ( |6.41 ) at n = 1, since 
after appropriate rescaling 

hk -> t^ihk, Hk ^ t^iTik (6.44) 

Ti can be identified with Aqcd- Note that, for doing this, one needs to use the Toda curve in the form 
(3.7). Eqs. ( 6.39| ), (6.41) at n = 1 are naturally interpreted in terms of the stress-tensor anomaly, 



+ ^40M ^ ^tr$2 = . . . + -d^pti {g^^G^"" + iG^^G^""^ , (6.45) 



since for any operator O 



^-r{0) = (/3Tr$^0) (6.46) 



aiogA 



Analogous interpretation for n > 2 involves anomalies of 14^„+i-structures. 



For n = 1, eq.(6.39), i.e. eq.(6.46) for O = I, has been derived in |162| in the form 



w 



dlogkqcD 



Eq.( 6.46| ) for O = hm, the analogue of ( 6.41 ) for n = 1 and any m = 1, . . . , A^c — 1, is 



dhr. 



(6.47) 



(6.48) 



Formulas of such a type first appeared for the SU{2) case in [|156[, and literally in this form for 



the SU{2) case in |157[ |158| and for the general SU{Nc) case in ]157[ . It was deduced from its 
modular properties and also from sophisticated reasoning in the framework of the Donaldson theory 
(the topological correlator is an object from the Donaldson theory, i.e. topological subsector of the 
Af — 2 SUSY Yang-Mills theory, which can be investigated with the help of the twisting procedure). 



In (6.48), one can substitute hm by any homogeneous function of /I's, e.g. by 7i„ 



Eq. (S.41) can be also rewritten as: 



T{a,t) - ^T^^^'{a,t) = --I—^^l^:^l^dfAog0E{O\T) (6.49) 



where the prepotential of the Generalized Kontsevich Model |163| 



jrGifMj-^l^^ _ — — t„itn'Hm+l,n+l 



(6.50) 



2N ■ 

m.n 

implicitly depends on the moduli a* through the coefficients of the polynomial P{X). 

Here we considered only the simplest example of the 4c? pure gauge theory. This approach was 



generalized to the theories with fundamental matter and with other gauge groups in [149, 164, 165|). 



Further development on the connections between Whitham hierarchies, r-functions of integrable hier- 



archies of the KP/Toda type and Donaldson topological theories can be found in |166|. These latter, 
very interesting results are still waiting to be re- interpreted in spirit of the approach of |162] to the 
Generalized Kontsevich Model. 

As far as we know, the adjoint matter case has not been done yet. 
6.4 Whitham hierarchy for SU{2): an example 

To conclude our discussion of Whitham hierarchies, let us illustrate the scheme developed above in 
the simplest case of 2-site periodic Toda chairQ One can restrict himself to the first two differentials 
drin and the first two times to, iij^ Then, the Whitham equations 

are reduced to the only equation 

■^^We do not discuss here how to calculate the second derivative of the prepotential in this case, since it is not too 
illuminating. The details can be found in Appendix B of ]82| . 

**^It can be done since the higher flows arc dependent on the two first flows (or, which is the same, u depends only on 
two independent combinations of times). Generically, for a complex curve of genus g, there are g + 1 independent flows. 



ddfln 



ddfli 



dti dto 

To write it down explicitly one should remember that 



(6.52) 



There are two independent differentials (see (6.23) and (6.24)) 



dflo = dz + ^ 



dfli = Adz + 71 — f A + ■^'l 
y V 2A/ 



(6.53) 



and two corresponding Whitham times io and ti. Here we introduced the notation dz = dw/w. 
• The coefficients 7i (i — 0, 1) are fixed, as usual, by vanishing of the correspondent A-periods 



dn, = 



(6.54) 



I.e. 



where 



7i 



a = j> dS 



a 
a 
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1 

a 



\dz, 



dX da 



A V 



du 



(6.55) 



(6.56) 



• The derivatives over moduli in (3.51) are taken at fixed values of z- variable, while the A- variable 
depends on moduli and this dependence is given by ( 3.18| ). This rule leads to the relations 



9A _ J_ ^ _ / - / ^ - cr 



du 2A' du du^~ 2A 



(6.57) 



Using these relations one easily gets from ( |6.5l| ) the following Whitham equation on moduli 



du , s,du 



(6.58) 

The function a{u) is an elliptic integral. Its explicit expression depends on choice of the cycles. 
With the choice used in i.e. when the A-cycles encircle the points A ~ \/u— 2K^ = r~ and 
A = Vu + 2A2 = r+. 



a^-r+E{k), a = -^K{k) (6.59) 

where K{k) and E{k) are complete elliptic integrals of the first and the second kinds respectively and 
the elliptic modulus is fc = 



The Whitham equation (6.58) has the general solution: 



u = F {to + a{u)ti) 

where F(x) is an arbitrary function. This solution can be also rewritten in the form 



(6.60) 
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^u) ^ to + a{u)ti (6.61) 

where $(u) is the function inverse to F{x). 

Thus, we constructed the solution to the Whitham hierarchy just straightforwardly. Now we do 
the same, using the solution for the differential dS ( 6.27 ) found in the previous subsection. We start 
first with only two first non-zero times so that dS has a particular form 



dS = todno + tidni = todz + tiXdz (6.62) 
Then, using formulas of the end of the previous section for the differentials (if2i^2, one can easily obtain 



from (6.35) for n = 0, 1 the following two equations: 



du 
du 



a 
a 
1 
a 



(6.63) 



Any solution to the equations ( 3.6S ) evidently solves simultaneously ( 6.58 ), but not vice versa. 
Let us use the manifest form of functions a{u) and cr(u) from ( 6.59| ): 



du 

du 
'dfo 



8 Eik) 

2tt 1 
'1^ K{k) 



(6.64) 



Then, the first equation of (6.63) has the solution |167| 



with arbitrary function c{to)- The solution to the both equations ( 3.63 ) takes the form 



a{u) 



const — t 







t^ 



(6.65) 



(6.66) 



This particular solution has to be compared with ( 6.61 ) with the constant function <l>(u) (since it has 
no inverse, the second form ( |6.60D does not exist for this concrete solution). 

In order to get the general solution ( |6.61 ) to the Whitham hierarchy, one has to require for all higher 



odd times in (3.27) (i.e. with n > 1) to be non-zero constants. Then, these higher times parameterize 



the general function $(u). Indeed, how it can be easily obtained from (6.35), the solution to the 



Whitham equations with general dS (6.27) is given by formula (6.61) with 



k=l 



(6.67) 



Therefore, the differential dS ( 3.27 ) actually leads to the general solution to the Whitham hierarchy, 
the higher times being just constants which parameterize the solution. 

To complete our consideration of particular system ( |6.63 ), we discuss what happens when one 
considers u as a function of three independent variables: to, ti and a. In this case, one still gets the 
same solution depending on tvjo independent variables. 



Indeed, since 



a = (j) dS — tia + to 



(cf. with eq.( 6.66| )) and using (6.56), (6.63) turns into 



(6.68) 



du du du 

oti a a a a 

du du 

dto da 



with the solution 



(6.69) 



to 



(6.70) 



where \E'(x) is an arbitrary function. Note that (|6.66| ) now can be understood just as the condition of 
constant a. 

Let us now note that a{u) for the given data is some known function of u. Then, the inverse 



(not arbitrary!) function u = ^'(a) can be rewritten in the form (6.70) using (3.6S) with independent 
variables a, to and ti. It means that one can not consider arbitrary function "^(x) with the particular 



anzatz of zero higher times in ( 6.27 ), i.e. solutions to the equations ( 6.69 ) should be additionally 
constrained. 

If now one forgets about the specific anzatz for dS with zero higher times and considers uncon- 



strained dS (|6.27| ), one should also consider the additional Whitham equation 



dd^n ddu! 



da 



dtn 



(6.71) 



that exactly gives rise to the equations (6.69) 



du 

dh ~~ 
du 
dTo 



du 
da 
du 

da 



(6.72) 



These equations are equivalent to (6.5§|) with u being a function of to ^ ol- It means that we again 



obtain the same solution (6.61) that describes u as a function of two independent variables ti and 
to — a. This perfectly matches the number g + n — 1 = 2 of independent variables. 



6.5 WDVV equations 

We started this review with discussing integrability of Seiberg-Witten theory and continued with 
Whitham hierarchies associated both with integrable and topological properties of Seiberg-Witten 
theory. Now we are going to discuss briefly a feature of the prepotential that has much to do more 
with topological theories, while is not that straightforwardly (to our present knowledge) related to 



integrability. This property briefly mentioned in the Introduction is the WDVV equations (1^), which 
are solved by the prepotentials of theories with the spherical bare spectral curve. 

Note that we still do not understand any special reason for the WDVV equations to emerge 
in topological and Seiberg-Witten theories, although, first, the latter two arc connected through the 



Whitham hierarchies (see the previous subsection) and there are also tight connections of the Whitham 



hierarchies with the WDVV equations |40 



To get some insight of general structure of the WDVV equations, we start with examples not 
related to Seiberg-Witten theory but covering the main fields where the WDVV equations have been 
met so far (that is, type A and type B 2d topological theories). 

The simplest non-trivial example of iVc = 3 WDVV (the equations at Nc — 1,2 are trivially 
satisfied) comes from the type B topological theory, namely, Af = 2 SUSY Ginzburg-Landau theory 
with the superpotential W'{X) — X'^ ~ q In this case, the prepotcntial reads as 



1 2 1 2 1 2 
F = -aia2 + -a^as + -0303 

and the matrices Fi (the third derivatives of the prepotcntial) are 



(6.73) 



/ 1 \ 

1 
V 1 / 



/ 1 \ 

1 
\0 q J 



F. 



/ 1 \ 

q 
\0 q J 



(6.74) 



One can easily check that these matrices do really satisfy the WDVV equations (1. 

The second example is the quantum cohomologies of (D-P^ (typc A topological cr-model). In this 
case, the prepotential is given by the formula Q 



1 



1 



2 ^ 2 2 ^ ^ ^ 3fc- 1)! 



(6.75) 



where the coefficients Nk (describing the rational Gromov-Witten classes) counts the number of the 
rational curves in (C-P^ and are to be calculated. Since the matrices F have the form 



/ 1 \ 
1 

V 1 y 



F, 



/ 1 \ 

1 F222 F223 
\ F223 F233 j 



F3 



/ 1 \ 

-^223 -P233 
\ i^233 -P333 / 



the WDVV equations are equivalent to the identity 



^333 — 



223 



^222-^233 



which, in turn, results into the recurrent relation defining the coefhcients N-^ 



(6.76) 



(6.77) 



a^h{ih-\)h{2a-h) 
(3^-4)! (3a -1)1(3^,-1)! 



E 



(6.78) 



The crucial feature of the presented examples is that, in both cases, there exists a constant matrix 
F\. Following |40j, one can consider it as a flat metric on the moduli space. In fact, in its original 
version, the WDVV equations have been written in a slightly different form, that is, as the associativity 
condition of some algebra. We shall discuss this later, and now just remark that, having distinguished 



(constant) metric rj = Fi, one can naturally rewrite (LS) as the equations 



CiCj 



CjCi 



(6.79) 



01 



for the matrices d = rj-^F^, i.e. {Ci)i = Tj^^Fuk- Formula I^J^ is equivalent to (|l|) with j = 1 



Moreover, this particular relation is already sufficient ||39|, 168 1 to reproduce the whole set of the 
WDVV equations (1.8). Indeed, since Fi — FiCi, we obtain 



F,F-'Fk = {QC^'C,) (6.80) 

which is obviously symmetric under the permutation i ^ j. Let us also note that, although the 
WDVV equations can be fulfilled only for some specific choices of the coordinates on the moduli 
space, they still admit any linear transformation. This defines the flat structures on the moduli space, 
and we often call ai flat coordinates. 

to be true, how we explain below. 



In fact, the existence of the flat metric is not necessary for (1 



Moreover, the Seiberg-Witten theories give exactly an example of such a case, where there is no 
distinguished constant matrix. This matrix can be found in topological theories because of existence 
their field theory interpretation where the unity operator is always presented. 

Now note that in the context of the two-dimensional Ginzburg- landau topological theories, the 
WDVV equations arose as associativity condition of some polynomial algebra. We will prove below 
that the equations in Seiberg-Witten theory have the same origin. Now we briefiy remind the main 
ingredients of this approach in the standard case of the Ginzburg-Landau theories. 

In this case, one deals with the chiral ring formed by a set of polynomials {$i(A)} and two co- 
prime (i.e. without common zeroes) fixed polynomials Q(A) and P{X). The polynomials <f> form the 
associative algebra with the structure constants given with respect to the product defined by 
modulo P': 



= C^$fcQ' + (*)P' C%^kQ! 

the associativity condition being 



(6.81) 



i.e. CiCj = Cjd, {Ci)l. — C^j. 



(6.82) 
(6.83) 



Now, in order to get from these conditions the WDVV equations, one needs to choose properly the 
flat moduh lEol: 



(p'/"dQ) , n = ord(P) 



i{n — i) 

Then, there exists the prepotential whose third derivatives are given by the residue formula 



Fijk — res 



P'=0 



(6.84) 



(6.85) 



On the other hand, from the associativity condition (6.83) and the residue formula (6.85), one obtains 
that 



F^,k^{Ci)]FQ,ik, i.e. C,^F,Fq} (6.86) 
Substituting this formula for d into ( |6.83 ), one finally reaches the WDVV equations in the form 



oo 



FiG ^Fj — FjG ^Fi, 

(6.87) 

G^Fq, 

The choice Q' — gives the standard equations ( |1.8| ). As we ah'eady mentioned, in two-dimensional 
topological theories, there is always the unity operator that corresponds to Q' = 1 and leads to the 
constant metric Fq' . 

Thus, from this short study of the WDVV equations in the Ginzburg-Landau theories, we can get 
three main ingredients necessary for these equations to hold. These are: 

• associative algebra 

• flat moduli (coordinates) 

• residue formula 

As we already know, in Seiberg-Witten theory only the first ingredient requires a non-trivial check, 
while the other two are automatically presented due to the proper integrable structures. 



6.6 WDVV equations in Seiberg-Witten theory 

Now we are ready to study the most interesting for us example of Seiberg-Witten theory and determine 
when the prepotential satisfies the WDVV equations. First, let us note that if the prepotential satisfies 
the WDVV equations, its leading perturbative part satisfies them too (since the classical quadratic 
piece does not contribute into the third derivatives). In each case it can be checked by straightforward 
calculations. 

We already discussed the perturbative prepotentials in detail in the previous section, see also 
which also contains a proof of the WDVV equations for the perturbative prepotentials. Here we just 
list some results of adding some newer statements. 

i) The WDVV equations always hold for the pure gauge theories Tpert — ^pert (including the ex- 
ceptional groups) 1^ In fact, in |^ it has been proved that, if one starts with the general prepotential 
of the form 



= ^ X! ~ "j)^ ("^ ~ + + "j)^ {ai + aj)) + 2 X! (6.88) 

the WDVV hold iff a+ = or = 0, 7y being arbitrary. A new result due to A.Veselov ]170( | shows 
that this form can be further generalized by adding boundary terms. For this latter case, however, 
we do not know the non-perturbative extension. 

ii) If one considers the gauge supermultiplets interacting with the nj matter hypermultiplets in 
the first fundamental representation with masses nia 



pert 



pert 



(6.89) 



^^The rank of the group should be bigger than 2 for the WDVV equations not to be empty, thus for example in 
the pure gauge G2-model they are satisfied trivially. It have been checked in 16£ ] (using MAPLE package) that 
the WDVV equations hold for the perturbative prepotentials of the pure gauge F4, and By models. Note that the 
corresponding non-perturbative Seiberg-Witten curves are not obviously hyporcUiptic. 



00 



(where r and K are some undetermined coefficients), the WDVV equations do not hold unless K = 
the masses are regarded as moduli (i.e. the equations (L8) contain the derivatives with respect 

to masses) and 



fF{m) = ^ X! (("^^ ~ ~ '^P^ 

a, (3 

for the SU (Nc) gauge group and 



(6.90) 



(6.91) 

■ ^ TO^logTOa 



Q,/3 

r(r + s) 



for other classical groups, s = 2 for the orthogonal groups and s = —2 for the symplectic ones. 

Note that at value r — ~2 the prepotential ( |6.89 ) can be considered as that in the pure gauge 



theory with the gauge group of the higher rank, rank(G) + nf- At the same time, at value r = 2, like 
Oi's lying in irrcp of G, masses to^'s can be regarded as lying in irrep of some G so that if G = v4„, 
C„, Dn, G = An, Dm C'n accordingly (this is nothing but the notorious (gauge group < — > flavor 



group) duality, see, e.g. [114|). These correspondences "explain" the form of the mass term in the 
prepotential f{m). 

iii) The set of the perturbative prepotentials satisfying the WDVV equations can be further ex- 
tended. Namely, one can consider higher dimensional SUSY gauge theories, in particular, Sd theories. 
Then, the perturbative prepotentials in the SU{Nc) theory (without adjoint matter) is given by for- 



mula ( 5.3£ ) with d — 5. One can check that these perturbative prepotentials satisfy the WDVV 



equations. Note that the presence of the cubic terms in ( ^.39 ) turns to be absolutely crucial for the 



WDVV equations to hold. This means that the coefflcient Cc in ( 5.43| ) is fixed by the equations. This 
is not surprising since that the WDVV equations arc not necessarily satisfied when Cc ^ Nc — Nf/2. 
Since it corresponds to degeneration of an integrable system, one fixes some moduli (masses) which 
are the parameters to vary in the WDVV equations (see |3|] for details). 

One can also consider other classical groups. Then, the perturbative prepotentials acquire the 



form (5.8) with all x^logx substituted by ^x^ + ^Lia (e "^^^^y One can easily check, along the fine of 



[ p9| that these prepotentials satisfy the WDVV equations. 

iv) If in the id theory the adjoint matter hypermultiplets are presented, i.e. J-pert = ^pert~^-^pert~^ 
fA{m), the WDVV equations never hold. At the same time, the WDVV equations are fulfilled for 
the theory with matter hypermultiplets in the symmetric/antisymmetric square of the fundamental 
representation^ iff the masses of these hypermultiplets are equal to zero. 



v) Our last example |122, p9l has the most unclear status, at the moment. It corresponds to the 5d 



theory with adjoint matter. It is described by the perturbative prepotential ( 5.69| ). This prepotential 



^■^ These hypermultiplets contribute to the prepotential 

1 \ - 



S X 2 



pert 



(6.92) 



O A 



satisfies tlie WDVV equations iff e = 7r/2 [122|. Moreover, it gives the most general solution for a 
general class of perturbative prepotentials J-pert assuming the functional form 



= ^ /(a • a), (6.93) 

where the sum is over the root system $ of a Lie algebra. The mystery about this prepotential is 
that, on one hand, it never satisfies the WDVV equations unless e = tt and we do not know if the 
corresponding complete non-perturbative prepotential satisfies the WDVV even for e = 7r/2. On the 



other hand, in the limit i?5 ^ and e ^ MR for finite M, (when the mass spectrum ( 3.84 ) reduces 
to the two points Af„ — and M„ = M), the theory is the four dimensional SYM model with A/" = 4 
SUSY softly broken to A/" = 2, i.e. includes the adjoint matter hypermultiplet. It corresponds to item 
iv) when the WDVV always do not hold. By all these reasons, this exceptional case deserves further 
investigation. 

From the above consideration of the WDVV equations for the perturbative prepotentials, one can 
learn the following lessons: 

• masses are to be regarded as moduli 

• as an empiric rule, one may say that the WDVV equations are satisfied by perturbative pre- 
potentials which depend only on the pairwise sums of the type (a^ ± bj), where moduli Oi and 
bj are either periods or masses^. This is the case for the models that contain either massive 
matter hypermultiplets in the first fundamental representation (or its dual), or massless matter 
in the square product of those. Troubles arise in all other situations because of the terms with 

± ± Cfe ± . . .. (The inverse statement is wrong - there are some exceptions when the WDVV 
equations hold despite the presence of such terms - e.g., for the exceptional groups.) 

Note that, for the non-UV-finite theories with the perturbative prepotential satisfying the WDVV 
equations, one can add one more parameter to the set of moduli - the parameter Aqcd that enters all 
the logarithmic terms as x^logx/AgcD- Then, some properly defined WDVV equations still remain 
correct despite the matrices J^^^^ no longer exist (one just needs to consider instead of them the 



matrices of the proper minors) |171|, see footnote ^ 

Now let us check the WDVV equations for full Seiberg-Witten prepotentials. As we already 
discussed, in order to derive the WDVV equations along the line used in the context of the Ginzburg- 
Landau theories, we need three crucial ingredients: flat moduli, residue formula and associative al- 
gebra. However, the first two of these are always contained in the SW construction provided the 
underlying integrable system is known. Thus, the only point to be checked is the existence of the 
associative algebra. The residue formula ( |3.1[ ) hints that this algebra is to be the algebra oi the 
holomorphic differentials dwi. In the forthcoming discussion we restrict ourselves to the case of pure 
gauge theory, the general case being treated in complete analogy. 

Let us consider the algebra fl-^ and fix three differentials dQ, du, dX E fl^. The product in this 
algebra is given by the expansion 

duJidLjj = C^.duJkdQ + {*)duj + {*)dX (6.94) 



^^This general rule can be easily interpreted in D-brane terms, since the interaction of branes is caused by strings 
between them. The pairwise structure (a^ it bj) exactly reflects this fact, and bj should be idcntiflcd with the ends 
of string. 



that should be factorized over the ideal spanned by the differentials doj and dX. This product belongs 
to the space of quadratic holomorphic differentials: 



-n^ en^ = ■ {dQ®duj®dX) (6.95) 

Since the dimension of the space of quadratic holomorphic differentials is equal to 3g — 3, the l.h.s. 



of ( 6.94 ) with arbitrary duJi's is the vector space of dimension 3g — 3. At the same time, at the 



r.h.s. of ( |6.94 ) there are g arbitrary coefficients C,^ in the first term (since there are exactly so many 



holomorphic f-differentials that span the arbitrary holomorphic 1-differential C^jdcuk), 5 — 1 arbitrary 
holomorphic differentials in the second term (one differential should be subtracted to avoid the double 
counting) and g — 2 holomorphic 1-differentials in the third one. Thus, totally we get that the r.h.s. 



of ( 6.94 ) is spanned also by the basis of dimension g + [g — 1) + (g — 2) 3(7 — 3. 

This means that the algebra exists in the general case of the SW construction. However, generally 
this algebra is not associative. This is because, unlike the LG case, when it was the algebra of 
polynomials and, therefore, the product of the two belonged to the same space (of polynomials), 
product in the algebra of holomorphic 1-differentials no longer belongs to the same space but to the 
space of quadratic holomorphic differentials. Indeed, to check associativity, one needs to consider the 
triple product oiVl^: 

-Vl} - n^ = n^- [dQf ® ■ duj ® Vt^ ■ d\ (6.96) 

Now let us repeat our calculation: the dimension of the l.h.s. of this expression is bg — b that is the 
dimension of the space of holomorphic 3-differentials. The dimension of the first space in expansion of 
the r.h.s. is the second one is — 4 and the third one is 2g — ^. Since g + (3(7 — 4) + {2g — A) = 6(7 — 8 



is greater than 5g — 5 (unless 5 < 3), formula ( 6.96 ) does not define the unique expansion of the 



triple product of fl^ and, therefore, the associativity spoils. 

The situation can be improved if one considers the curves with additional involutions. As an 
example, let us consider the family of hyperelliptic curves: j/^ = -Po^2g+2(A). In this case, there is the 
involution, (t : ?/ — > —y and fi^ is spanned by the cr-odd holomorphic 1-differentials ^ ^ , i = I, ...,g. 
Let us also note that both dQ and dio are cr-odd, while dX is tr-even. This means that dX can be 
only meromorphic on the surface without punctures (which is, indeed, the case in the absence of mass 
hypermultiplets). Thus, 

dX omits from formula ( ]6.94D that acquires the form 



nl = nl-dQ®nl- dw (6.97) 

where we expanded the space of holomorphic 2-differentials into the parts with definite cr-parity: 
VL^ = Vl\ © , which are manifestly given by the differentials - — ^.f^^ , i = l,...,2g— 1 and - — , 
i — 1, (7 — 2 respectively. Now it is easy to understand that the dimensions of the l.h.s. and r.h.s. 



of (6.97) coincide and are equal to 2g — 1. 



Analogously, in this case, one can check the associativity. It is given by the expansion 

n^^ = n\- {dQf ®VL\-dijj (6.98) 

where both the l.h.s. and r.h.s. have the same dimension: "ig — 2 — g + {2g — 2). Thus, the algebra 
of holomorphic 1-differentials on hyperelliptic curve is really associative. This completes the proof of 
the WDVV equations in this case. 



Now let us briefly look at cases when there exist the associative algebras basing on the spectral 
curves discussed in the previous section. First of all, it exists in the theories with the gauge group 
An, both in the pure gauge 4d and 5d theories and in the theories with fundamental matter, since, 
in accordance with the previous section, the corresponding spectral curves are hyperelliptic ones of 
genus n. 

The theories with the gauge groups SO{n) or Sp{n) are also described by the hyperelliptic curves. 
The curves, however, are of higher genus 2n — 1. This would naively destroy all the reasoning of 
this section. The arguments, however, can be restored by noting that the corresponding curves (see 
( 3.55 )) have yet another involution, p : A ^ —A. This allows one to expand further the space of 

holomorphic differentials into the pieces with definite p-parity: il^ = il^ ® ft^^ etc. so that the 

proper algebra is generated by the differentials from il^ . One can easily check that it leads again to 

the associative algebra. 

Consideration is even more tricky for the exceptional groups, when the corresponding curves are 
looking non-hyperelliptic. However, additional symmetries should allow one to get associative algebras 
in these cases too. 

There are more cases when the associative algebra exists. First of all, these are 5d theories, 
with and without fundamental matter One can also consider the SYM theories with gauge 

groups being the product of several factors, with matter in the bi- fundamental representation p| . 
As we discussed earlier, these theories are described by sl{p) spin chains and the existence of the 
associative algebra in this case has been checked in [172|. 

The situation is completely different in the adjoint matter case. In four dimensions, the theory is 
described by the Calogero-Moser integrable system. Since, in this case, the curve is non-hyperelliptic 
and has no enough symmetries, one needs to include into consideration both the differentials du and 
dX for algebra to exist. However, under these circumstances, the algebra is no longer to be associative 
how it was demonstrated above. This can be done also by direct calculation for several first values of 
Nc (see M). This also explains the lack of the perturbative WDVV equations in this case. 



6.7 Covariance of the WDVV equations 

After we have discussed the role of the (generalized) WDVV equations in SYM gauge theories of the 
Seiberg-Witten type, let us briefly describe the general structure of the equations themselves. We look 
at them now just as at some over-defined set of non-linear equations for a function (prepotential) of 
r variables^ (times), F{t^), i — I, . . . ,r, which can be written in the form (6.87) 



FiG Fj — FjG Fi. 



G^Y^T^'^Fk, Vi,j = l,...,r and \/ii''{t) 



(6.99) 



fe=i 



d^F 



Fi being r x r matrices {Fi)jk — F^ijk = gtiQpgtk ^""^"^ "metric" matrix G is an arbitrary linear 
combination of F^'s, with coefficients r]'^{t) that can be time-dependent. ^ 



The WDVV equations imply consistency of the following system of differential equations |173| 



''^We deliberately choose different notations for these variables, t instead of a in the gauge theories, in order to point 
out more general status of the discussion. 
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We already remarked in the previous subsection that one can add to the set of times (moduli) in the WDVV 



equations the parameter h^QCD [L71|. In this case, the prepotential that depends on one extra variable t'-^ = K can be 



Contracting with the vector ry' (t) , one can also rewrite it as 



dtp 
dt 



(6.100) 



(6.101) 



where 



C, 



(6.102) 



(note that the matrices Ck and the differential D depend on choice of {ri\t)}, i.e. on choice of the 



metric G) and (6.99) can be rewritten as 



[C„Q]=0, Vz,j 



(6.103) 



As we already discussed, the set of the WDVV equations ( |6.99 ) is invariant under linear change 
of the time variables with the prepotential unchanged According to the second paper of ]4C|] 



and especially to [174|, there can exist also non-linear transformations which preserve the WDVV 



structure, but they generically change the prepotential. In |175|, it is shown that such transformations 



are naturally induced by solutions of the linear system (6.100) 



Fit) 



so that the period matrix remains intact: 



(6.104) 



d^F 



jp . — jr^^ 



(6.105) 



Now let us make some comments. 



1. As explained in |173|, the linear system (6. IOC) has infinitely many solutions. The "original" 
time- variables are among them: ip^{t) = f . 

naturally considered as a homogeneous function of degree 2: 



,r) = (t°fF{e/t°), 



see |155, 5^1 general theory. As explained in |171|, the WDVV equations (6.99) for F{t^) can be also rewritten 

in terms of jF(t^): 

J'lG'^^Fj = Tjg-^J'l, VI,J = 0,l,...,r;{r,^(t)} 
where this time jFj are (r + 1) X (r + 1) matrices of the third derivatives of and 



k=o 



V^^Fk, Q-^ = (detg)g-l 



Note that the homogeneity of .7-" implies that t^-derivatives are expressed through those w.r.t. t', e.g. 



-7^ i-^ 
Aik'' 



-r^Umt'^th"' etc. 



Thus, all the "metrics" Q are degenerate, but Q^^ are non-degenerate. One can easily reformulate the entire present 
section in terms of J^. Then, e.g., the Baker- Akhiezer vector-function ip{i) should be just substituted by the manifestly 
homogeneous (of degree 0) function ipi^' /i^)- The extra variable should not be mixed with the distinguished "zero- 
time" associated with the constant metric in the 2d topological theories which generically does not exist (when it does, 
see comment 3 below, we identify it with f). 



2. Condition ( |6.105 ) guarantees that the transformation ( |6.104 ) changes the hnear system ( 3.10C| ) only 



by a (matrix) multiplicative factor, i.e. the set of solutions {ip^(t)} is invariant of (6.104). Among other 



things this implies that successively applying ( 6.104 ) one does not produce new sets of time- variables. 

3. We already discussed that, in the case of 2d topological models 174|, there is a distinguished 

time- variable, say, , such that all Fij^. are independent of t^: 



d_ 



ijk 



yi,j,k = i,. 



(6.106) 



(equivalently, -^Frjk — Vi, j, k). Then, one can make the Fourier transform of (6.100) with respect 
to and substitute it by the system 



_d_ 



(6.107) 



where -^^ (i"^, . . . , f ^) = / "0^ (^"^j . . . , T \ r)e^*'^(ir. In this case, the set of transformations ( 3.104 ) 
can be substituted by a family, labeled by a single variable z: 



In the limit z - 
transformation 



f tl=r.{t) (6.108) 
and for the particular choice of the metric, G = F^., one obtains the particular 



dP ~ ' 



const, 



(6.109) 



discovered in |174|. (Since — djC, one can also write = Clh'', CI = [F^ ) Fzfc.) 



4. Parameterization like ( 6.1G9| ) can be used in the generic situation ( |6.104[ ) as well (i.e. without 



distinguished ^'"-variable and for the whole family (3.104)), the only change is that h'' is no longer a 
constant, but a solution to 



id,-Dc,y,h'^^o 



(6.110) 



{h^ — Dip^ is always a solution, provided ip'' satisfies (6.101)). 

Note also that, although we have described a set of non-trivial non-linear transformations which 
preserve the structure of the WDVV equations ( |6.99 ) , the consideration above does not prove that 



all such transformations are of the form (6.104), (3.105). Still, (6.104) is already unexpectedly large, 
because ( |6.99 ) is an over-defined system and it could seem to be very restrictive, if to have any 
solutions at all. 

To conclude this short survey of the prepotentials and WDVV equations, let us emphasize that 
there are a lot of problems to be solved. We already mentioned the problem of lack of the WDVV 
equations for the Calogero-Moser system. The way to resolve this problem might be to construct 



higher associativity conditions like it has been done by E.Getzler in the elliptic case |176], that is to 
say, for the elliptic Gromov-Witten classes. Even if this will succeed, it is not an explanation why at 
all generalized WDVV equations should be satisfied for the Seiberg-Witten prepotentials. 



on 



7 Dualities in the many-body systems and gauge theories 



In this section, we are going to discuss duality, a property of integrable systems which has much to 
do with the possibihty of exact solving SUSY gauge theories at low energies and with the possibility 
of re- interpreting them in terms of strings/branes (see also the Introduction). This property also will 
allow us to fill in the right bottom corner of table 1, i.e. to construct a double elliptic system. 

In fact, as in string theory there are several dualities (S-, T- and U-dualities), as in integrable 
systems there are several kinds of duality [177|. We briefly survey all of them here, and then concen- 
trate on the most interesting, important and still remaining mysterious one which can be described 



via canonical transformations [178, Note that all the dualities at the moment arc attributed only 
with the Hitchin type systems, i.e. the Seiberg-Witten theories with adjoint matter, while spin chains 
are still not embedded into this treatment. 

Let us point out that the dualities under discussion interchanges different theories, while there 
is also a duality that interchanges descriptions of the same integrable system. This is the duality 
between Nc x Nc and 2x2 Lax representations, we discuss it in the next section. 



7.1 D-branes, gauge theory and separation of variables 

Thus, similarly to branes/string theories, the Hitchin type systems are involved into three different 
dualities. Let us start with considering a counterpart of the T-duality which can be identified with 



the separation of variables in dynamical systems |177|. The separation of variables (SoV) in integrable 



system allows one to reduce the system with N degrees of freedom to N identical systems with one 
degree of freedom each. 

Recently, E.Sklyanin formulated a recipe for the SoV in a large class of quantum integrable models 
[ |179| . The method is to study the pole dynamics of the Baker-Akhiezer function ^'(z. A): 

L(e)*(e,A) = A(0*(e,A) (7.111) 

with some clever choice of normalization which is the main, difhcult and artistic part of the method. 
The poles of 5'(^, A) together with the eigenvalues Xi — A(^i) are then identified with the (canonically 
conjugated) separated variables. 



In | ]177[ there was constructed a geometrical picture behind the recipe (see also [180, |181| ). Their 
construction of the separation of variables in integrable systems on moduli spaces of holomorphic 
bundles with some additional structures that can be described as a symplectomorphism between the 



moduli spaces of bundles (more precisely, torsion free sheaves) having different Chern classes. In [177|, 
there were discussed three equivalent realizations of the Hitchin integrable systems. 

The first one uses non-Abelian gauge fields on the base curve E embedded into a symplectic surface 
S. The phase space of the integrable system is the moduli space of the stable pairs {£, (/)), where £ is 
the rank r vector bundle over E of degree while is the holomorphic section of (8) End(£). In 
the second realization, one deals with the moduli space of pairs (f , </>), where C is the curve (divisor) 
in S which realizes the homology class r[E] and C is the line bundle on C. The third realization is 
the Hilbert scheme of points on S with the number of points determined by the rank of the gauge 
grouppl. 
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The r ole o f the Hilbert schemes in the context of Hitchin systems was first pointed out in 1 182 ] (see als o the orig inal 



context of |p_24]) for the surfaces without marked points and then generalized to the Calogero systems in [183, 177]. 



The equivalence of these three reahzations can be interpreted in string theory terms. To this end, 
let us specify the moduli space Mv of stable torsion free coherent sheaves £ on S. Let As = I — [pt] G 

H*{S,X) be the v4-roof genus of 5. The vector w = C/i (5)^1^ = {r;w;d-r) e Z)> G T^.i^ 

corresponds to the sheaves with the Chern numbers: 

c/io(f) = re 

chi{£) ^ w e H\S;X) (7.112) 
ch2i£) = deH^iS;1.) 

Now let us describe geometrically string theory objects. In type IIA string theory compactified 
onto S", there are Dp-branes with even p which wrap different supersymmetric cycles in S labelled by 
V £ H*{S,'^). The BPS states correspond to the cohomology classes of the moduli spaces, of 
brane configurations. These moduli spaces can be identified with the moduli spaces of appropriate 
sheaves, M^j. 

The string theory compactified on S, has the following vacuum moduli space 

Ma = (r-*'20) \0{4, 20; i?)/0(4; R) x 0(20; R) 

where the arithmetic group ©(F^'^") is the group of discrete automorphisms that maps between the 
states with the different v. The if^ is the only invariant of the group action. 

Now, the equivalence of the first and the second realizations of Hitchin systems corresponds to the 
physical statement that the bound states of N Z?2-branes wrapped around S are represented by a 
single Z)2-brane which wraps a holomorphic curve C which is an A'^-sheet covering of the base curve 
E. The equivalence of the second and the third descriptions is naturally attributed to T-duality. 

Let us mention that the separation of variables above provides some insights on the Langlands du- 
ality which involves the spectrum of the Hitchin Hamiltonians. An attempt to reformulate Langlands 



duality as a quantum separation of variables has been presented in [ 185 |, while in |178] the authors 
proposed to identify the proper classical version of the Langlands correspondence with the map of the 
phase space of the Hitchin system to the Hilbert scheme of points on four-dimensional manifold. 



7.2 S-duality 

Now let us see that S-duality well-established in field theory also has a clear counterpart in the 
holomorphic dynamical system. In field theory, S-duality just maps one theory onto another so that 
the weak coupling regime of one theory maps onto the strong coupling regime of the other one and 
vice versa. Therefore, a hint of what is S-duality in integrable systems can come from Seiberg-Witten 
theory where this duality just maps the weak coupling, perturbative limit of the gauge theory (due 
to the asymptotic freedom at large energies, where the gauge bosons are almost massless, with their 
masses proportional to a^) onto the strong coupling limit (where the BPS monopole is almost massless, 
with the mass proportional to af ). This means that S-duality interchanges theories with small a and 
small aj). Note now that, in the integrable set-up, formula ( |l.3| ) it just corresponds to interchanging 
A- and i?-cycles. Therefore, the duality transform is just the SP{2g, 2) group acting on the 2g cycles 
and the whole set of variables {a,j, af} is just local coordinates on the SP{2g, ^) bundle. 

It is an easy exercise on the properties of the period matrix of Abelian varieties to check that the 
two form 

m 



da' Adaf (7.113) 
vanishes. Therefore one can always locahy find a function T ^ prepotential such that 

To illustrate this scheme with some formulas, let us look at the two-body system, relevant for the 
SU (2) Af = 2 supersymmetric gauge theory. It is the periodic Toda chain with the Hamiltonian (see 
s.3.2) 

i?=y + cos(q) (7.115) 



The action variable is given by one of the periods of the differential pdq (3. IE) on the curve ( p. 18 ). 
Then, the two periods a and obey the Picard-Fuchs equation 



Ju^ 4(m2 - 1)7 V a^. 
which can be used to obtain the asymptotic expansions of the action variable and the prepotential 
in the vicinity ofw = oo or u = ±1. The S-duality manifests itself in the fact that near u = oo 
(which corresponds to high energies in the two-body problem and also to the perturbative regime of 
the SU{2) gauge theory) the appropriate action variable is a, while near u = 1 (which corresponds to 
the dynamics of the two-body system near the top of the potential and to the strongly coupled SU (2) 



gauge theory) it is a . The monodromy invariant combination of these periods |162 



aa^ -2T^u (7.116) 

(which originates from the periods of Calabi-Yau manifolds, on one hand, and from the properties of 
anomaly in field theory, on the other) can be chosen as a global coordinate on the space of integrals 
of motion. At u — )■ oo the prepotential has an expansion of the form 



T ^ ^ulogu +...r^ a^loga + ^ /na^-^" 

n 

Let us emphasize that S-duality maps the dynamical system to itself, just at a different energy. 
7.3 "Mirror" symmetry in dynamical systems 

Now we come to the most interesting type of duality, a duality between pairs of dynamical systems 



(which is self-duality for some systems) ||6^, 124, 178, This duality plays a role similar to the 
mirror symmetry. 

To start with, let us see how this symmetry is formulated in field theory. The initial motivation 



comes from the 3d theory example [186|, where the mirror symmetry interchanges the Coulomb and 
Higgs branches of the moduli space. The specifics of three dimensions is that both the Coulomb and 
Higgs branches are hyperkahler manifolds and the mirror symmetry can be formulated as a kind of 
hyperkahler rotation. 



no 



Taking into account the relation between dynamical systems and low-energy effective actions, one 
can formulate a counterpart of this mirror symmetry for gauge theories in different dimensions (see 



further details in |175]). 

Let us consider a system with the phase space being a hyperkahler four dimensional manifold. Let 
this manifold at most involve two tori, i.e. be the elliptically fibered K3 manifold. One torus is where 
momenta live, and the second one is where coordinates do. The duality at hands actually interchanges 
momentum and coordinate tori. However due to the nontrivial fibration this interchange is far from 
being naive. Note that typically one deals with degenerations of these tori and only in the double 
elliptic case two elliptic moduli survive. 

One can also consider a theory which is self-dual. This integrable system is called Dell (double 
elliptic) systemp^ and describes the Seiberg-Witten theory for the six-dimensional theory with adjoint 
matter. 

The other cases correspond to some degenerations of Dell system. Say, the degeneration of the 
momentum torus to C jZ^ leads to the five-dimensional theory, while the degeneration to i?^ leads to 
the four-dimensional theory. Since the modulus of the coordinate torus has meaning of the complexified 
bare coupling constant in field theory, interpretation of the degenerations of the coordinate torus is 
different. In particular, the degeneration to the cylinder corresponds to switching off the instanton 
effects, i.e. to the perturbative regime. 

We shall see below that the "mirror" symmetry maps theories in different dimensions to each 
other. Instanton effects in one theory "get mapped" into an additional compact dimension in the 
dual counterpart. We discuss here the pg-duality mostly for classical integrable systems only. Since 
the wave functions in the Hitchin like systems can be identified with some solutions to the KZ or 
qKZ equations, the quantum duality would imply some relations between solutions to the rational, 
trigonometric and elliptic KZ equations. Recently the proper symmetries for the KZ equations were 



discussed in |75 



Now let us formulate the same mirror symmetry in terms of dynamical system. The idea of duality 
here expresses a relationship between two completely integrable systems S\ , 5*2 on a fixed symplectic 
manifold with given symplectic structure (Af, lS) and goes back to |^ . We say the Hamiltonian 
systems are dual when the conserved quantities of S\ and 5*2 together form a coordinate system for 
M. Consider for example free particles, h'^^^ = 'YliiV\l^- For this system the free particles momenta 
are identical to the conserved quantities or action variables. Now consider the Hamiltonian = 

c^jk with conserved quantities g^. Together {pi, qi\ form a coordinate system for phase space, and 
so the two sets of Hamiltonians are dual. Duality then in this simplest example is a transformation 
which interchanges momenta and coordinates. For more complicated interacting integrable systems 
finding dual Hamiltonians is a nontrivial exercise. Note that this whole construction manifestly 
depends on the particular choice of conserved quantities. A clever choice may result in the dual 
system arising by simply interchanging the momentum and coordinate dependence, as in the free 
system. 

Some years ago Ruijsenaars ]62| observed such dualities between various members of the Calogero- 
Moser and Ruijsenaars families: the rational Calogero and trigonometric Ruijsenaars models were 
dual to themselves while trigonometric Calogero model was dual with the rational Ruijsenaars system. 
These dualities were shown by starting with a Lax pair L = L(p, q) and an auxiliary diagonal matrix 
A — v4(q). When L was diagonalized the matrix A became the Lax matrix for the dual Hamiltonian, 

*^In fact, this is an example of the Dell system with 2 degrees of freedom. One may equally consider the Dell system 
with arbitrary many degrees of freedom, see ss.6.3,6.4. 



while L was a function of the coordinates of the dual system. Dual systems for a model possessing 
a Lax representation are then related to the eigenvalue motion of the Lax matrix. We discuss this 
approach in more details below. 

Another approach to finding a dual system is to make a canonical transformation which sub- 
stitutes the original set of Poisson-commuting coordinates g^, {<Zi,(7j} = 0, by another obvious 
set of the Poisson-commuting variables: the Hamiltonians hi{p,q) or, better, the action variables 
ai{h) = ai{'p,q). It will be clear below that in practice really interesting transformations are a little 
more sophisticated: hi are identified with certain functions of the new coordinates (these functions 
determine the Ruijsenaars matrix A{q)), which - in the most interesting cases - are just the same 
Hamiltonians with the interactions switched-off. Such free Hamiltonians are functions of momenta 
alone, and the dual coordinates substitute these momenta, just as one had for the system of free 
particles. 

The most interesting question for our purposes is: what are the duals of the elliptic Calogero and 
Ruijsenaars systems ? Since the elliptic Calogero (Ruijsenaars) is rational (trigonometric) in momenta 
and elliptic in the coordinates, the dual will be elliptic in momenta and rational (trigonometric) 
in coordinates. Having found such a model the final elliptization of the coordinate dependence is 
straightforward, providing us with the wanted double-elliptic systems. 



7.4 Duality as the canonical transformation: two-body system 

The way of constructing dual systems via a canonical transformation is technically quite tedious for 
many-particle systems, but perfectly fits the case of SU{2) which, in the center-of-mass frame, has 
only one coordinate and one momentum. In this case the duality transformation can be described 
explicitly since the equations of motion can be integrated in a straightforward way. Technically, 
given two Hamiltonian systems, one with the momentum p, coordinate q and Hamiltonian h{p, q) and 
another with the momentum P, coordinate Q and Hamiltonian H{P, Q) we may describe duality by 
the relation 

Kp,q)^nQ), 

(7.117) 

H{P,Q)^F{q). 
Here the functions f{Q) and F{q) are such that 

dP /\ dQ ^ -dp /\ dq, (7.118) 
which expresses the fact we have an (anti-) canonical transformation. This relation entails that 

F'iq)^ = /'(Q)^^g^. (7.119) 

At this stage the functions f{Q) and F{q) are arbitrary. However, when the Hamiltonians depend 
on a coupling constant and are such that their "free" part can be separated and depends only on 
the momenta,^ the free Hamiltonians provide a natural choice for these functions: F(q) = h(){q) and 
f{Q)^HoiQ) where 

Note that this kind of duality relates the weak coupling regime for h{p, q) to the weak coupling regime for H{P, Q). 



With such choice the duaUty equations become 



(7.121) 



hoiQ) = h{p,q), 
Ho{q)^H{P,Q), 

dh{p,q) dH{P,Q) 
-Ho(q) = K{Q)- 



(7.122) 



dp 



dP 



Free rational, trigonometric and elhptic Hamiltonians are /io(p) — ^o(p) = coshp and ho{p) 
cn(p|fc) respectively. 

Note that from the main duality relation, 



Ho{q)^HiP,Q) 



it follows that 



which together with (7.122) implies 



dq 



dq 
dP 



1 dHjP, Q) 
H'M dP 



1 dh{p, q) 



Q KiQ) dp 

When compared with the Hamiltonian equation for the original system, 



dq 
dt 



dhjp, q) 
dp 



(7.123) 



(7.124) 



(7.125) 



(7.126) 



we see that P — hQ{Q)t is proportional to the ordinary time- variable t, while Hq^Q) — h(p,q) = E 
expresses Q as a function of the energy E. This is a usual feature of classical integrable systems, 
exploited in Seiberg-Witten theory p^: in the SU{2) case the spectral curve q(t) can be described by 



E. 



(7.127) 



where p is expressed through dq/dt and q from the Hamiltonian equation dq/dt = dH/dp. In other 
words, the spectral curve is essentially the solution of the equation of motion of integrable system, 
where the time t plays the role of the spectral parameter and the energy E that of the modulus. 

Let us consider several simple examples. The simplest one is the rational Calogero. In this case, 



the duality transformation connects two identical Hamiltonians, (7.12C). Therefore, this system is 



For example, in the rational Calogero case 



h{p, q) = 
H(P, Q) : 



2 g2 



,2 q2 



(7.120) 



We recall that, on the field theory side, the coupling constant v is related to the mass of adjoint hypermultiplet and 
thus remains unchanged under duality transformations. 



self-dual. Somewhat less trivial example is the trigonometric Calogero-Sutherland model. It leads to 
the following dual pair 



Hp, 9) = it + = -IT 

i sm a i 



sm q 

(7.128) 



/ 2iy2 

H{PM) = Wl- -^cosP^cosg 

One may easily recognize in the second Hamiltonian the rational Ruijsenaars Hamiltonian. Thus, the 
Hamiltonian rational in coordinates and momenta maps onto itself under the duality, while that with 
rational momentum dependence and trigonometric coordinate dependence maps onto the Hamiltonian 
with inverse, rational coordinate and trigonometric momentum dependencies. Therefore, we can really 
see that the pg-duality exchanges types of the coordinate and momentum dependence. This implies 
that the Hamiltonian trigonometrically dependent both on coordinates and momenta is to be self-dual. 
Indeed, one can check that the trigonometric Ruijsenaars- Schneider system is self-dual: 



hip, g) = W 1 — ^ cosp = cos Q 



(7.129) 



/ 2zy2 

H{P,Q) ^ Jl- — cosP^cosg 
Now we make the next step and introduce into the game elliptic dependencies. 



7.5 Two-body systems dual to Calogero and Ruijsenaars 

We begin with the elliptic Calogero Hamiltonianpl 



*^We use below the Jacobi functions. Recall the standard relations: 



sn(o)= /_!i3_ ^ ^ J_eiiX9!l) 



cn(g) = J^^®^ = W-^ = W-^H1M (7.130) 



Here the Jacobi moduli k"^ and fc'^ = 1 — are the cross-ratios of the ramification points of the (hyper-) elliptic 
representation of the torus, 

3 3 

= 'W{^-ea{T)), ^ea = 0, x = p{q), y = ip'(g). (7.131) 

a— 1 a=l 

Then 

^2 ^ £23 ^ ^ ^ e!o(0|r) ^,2 ^ 1 _ ^2 ^ fH = = ^01 (0|^) .7 l32^ 

ei3 i?|(0) e4,(0|r)' ei3 i?|(0) ef,^{i)\r) ^ ' 

and 

eij=ei-ej, q = 2K{k)q, q = 2ojq, eis = — (7.133) 

where K{k) is the complete elliptic integral of the first kind. 



and seek a dual Hamiltonian elliptic in the momentum. Thus ho{p) = ^ and we seek H{P, Q) = HqIq) 
such that HQ{q) = cn(g|fc). Eqs.(7.123) become 
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sn2((7|fc) ' 



cn{q\k)^H{P,Q), p-cn'{q\k)^Q 



dHjP, Q) 
dP ' 



(7.135) 



Upon substituting 



cx)'{q\k) = -8YY{q\k)dn{q\k) = - ^ {1 - m){k'^ + k^m), (7.136) 
(this is because sn^q = 1 — cn^q, dn^q = fc'^ + fc^cn^q, fc'^ + k'^ — I and cnq = H) we get for ( [7.135 ): 



dP 



(1 -i/2)(fc'2 + fc2^2^ 



(7.137) 



Now from the first eqn.(7.135) can be expressed through Q and sii^{q\k) = 1 — C'n'^{q\k) = 1 — 
as 



p2 ^ 



Q2(l-i/2)' 



so that 



(7.138) 



Therefore H is an elliptic function of P, namely |178, p2| 



(7.139) 



i/(P, Q) = cn(g|fc) = a(Q) • cn P^k'^ + k^a^Q) 



with 



ka{Q) 



(7.140) 



(7.141) 



In the limit = 0, when the interaction is switched off, aijq) — 1 and i?(P, Q) reduces to 



Hq{P) — cn(P|fc), as assumed in ( 7.135 ). 

We have therefore obtained a dual formulation of the elliptic Calogero model (in the simplest SU{2) 
case). At first glance our dual Hamiltonian looks somewhat unusual. In particular, the relevant elliptic 
curve is "dressed" : it is described by an effective modulus 



ka(Q) 



ka{Q) 



^/k^TWiQ) ^l-k^l-a^Q)y ^^'^^^^ 

which differs from the "bare" one A; in a Q-dependent way. In fact k^ff is nothing but the modulus 
of the "reduced" Calogero spectral curve ^ . 



Now one can rewrite (7.140) in many different forms, in particular, in terms of ^-functions, in 
hyperelliptic parameterization etc. This latter is rather illuminating and allows one to picture mapping 
between the bare and dressed tori. Therefore, we refer the reader for the details to the paper | ]9^ . 

Let us now construct the action-angle variables for the 2-particle Calogero system and see how 
they related with the dual variables P and Q. 



Since this system possesses one degree of freedom, we proceed similarly to the 2-particle Toda 
case of s.3.2 and construct the action variable, a as the integral of differential dS — pdq on the level 
submanifold of the phase space. Then, using formula (7.135), one obtains 




dq 



(7.143) 



which expresses a through Q quite transcendentally. Using this relation, one obtains 

Q 



da = dQ 



zdq 



(7.144) 



sn^ (q) 



We already met this integral in s.5.6, formula ( 5.51 ). It can be manifestly calculated giving 



da = 



KikeffiQ)) 



-.dQ 



(7.145) 



Kik) ^k'-^ + k^a^Q) 

where we used the formulas of footnote ^ for the both bare and dressed tori as well as the definitions 
( [7.14l[ ) and ([7.142|) . Note that da = dQ at large Q ^ V2E. 

Since the coordinate on the Jacobian differs from the argument of the Jacobi function by a factor 
2K{k) (see footnote |4^), one expects for the angle variable p-^'^'^ (looking at the Hamiltonian ( 7.140[ )) 



Jac 



Py^k'^ + k^a^{Q) 



(7.146) 

2K{k,ff{Q)) 

This is, indeed, the case, up to inessential constant normalization factor 2K{k) (emerging here since 
the original Calogero symplectic structure is defined here as dp A dq with q being the argument of the 
Jacobi function, but not function) , since it follows from ( 7.145 ) and ( 7.146 ) that 

1 



dP AdQ 



dp-^"" A da 



2K{k) 



(7.147) 



To conclude our discussion of the Calogero model, we make some comments on its solution. Ac- 
cording to the argument of the previous subsection, our Hamiltonian (7.140) should be simply related 
to the solution q{t) of the equations of motion of the Calogero Hamiltonian, which in the case of SU (2) 
are immediately integrated to give 



H 



Cal 



^9 1 p 



(7.148) 



More explicitly, the equation 



has a solution |187] 



dq 
'dt 



^ = J2E- 



2v'^ 
sv?[q\k) 



cn(q|fc) = y^l - • cn [uj 2E - 2v'^k'^ 



(7.149) 



E - z/2fc2 



(7.150) 



This may be derived straightforwardly by differentiating both sides and applying ( 7.136| ). Note that 
the Calogero equation ( 7.149| ) and the family of Calogero spectral curves are essentially independent 



of the value of coupling constant v^: it can be absorbed into rescaling of moduli (like E) and the 
time- variables (like t). 

In order to see that ( 7.150 ) is identical to ( [7.140 ) one needs to put E = 1^2. and make the 
rescaling P = h'Q{Q)t = Qt. With these substitutions we find that 



and 




2 = arat(Q), 



(7.151) 



tV'2E~2v^k^ =P\ 1- 



E-iy^ 



= k\ 



l-2i/2/Q^ 



2 =P.^fc'2 + fc2a2^,(Q), 

karatiQ) 



We then see that ( |7.15G| ) is identical to ( |7.140| ). 

We remark that the relevant symplectic structure here ia^ 



(7.152) 



dE hdt^ QdQ hdt = -dP A dQ. 



(7.153) 



All of the above formulae are straightforwardly generalized from the Calogero (rational-elliptic) 
system to the Ruijsenaars (trigonometric-elliptic) system. The only difference ensuing is that the q- 
dependence of the dual (elliptic-trigonometric) Hamiltonian is now trigonometric rather than rational 



[ p2| (see also |178| for a more involved formula in the action-angle variables; this, however, makes 
formula not convenient for the double elliptic generalization): 



trig 



(9) 



1 



2i.2 



sinh^ q 



(7.154) 



Rather than giving further details we will proceed directly to a consideration of the double-elliptic 
model. 



7.6 Double-elliptic two-body system 



In order to get a double-elliptic system one needs to exchange the rational Q-dependence in (7.135) for 
an elliptic one, and so we substitute a^^f{Q) by the obvious elliptic analogue a1ii{Q) = 1 



sn2(Q|fc) 



Moreover, now the elliptic curves for q and Q need not in general be the same, i.e. k ^ k. 



Instead of (7.135) the duality equations now become 



^" In what follows we prove that not only the elliptic-rational (the dual of the elliptic Calogero model) but also 
the elliptic-trigonometric (the dual of the elliptic Ruijsenaars model) and the elliptic-elliptic (our new double-elliptic) 
Hamiltonians have the same form (7.150), but the latter with the identifications E = ^ sinh^ Q and E = isn^(Q|fc). 
Thus they are also related to Calogero equation {7.1A9). However, the relevant symplectic structures - which are always 
given by dP A dQ = hQ{Q)dt A dQ = dho{Q) A dt - are no longer equivalent to dE A dt (since E 7^ ho{Q), i.e. E is no 
longer associated with the proper Hamiltonian). 



nn 



cniq\k) = HiP,Q\k,k), 
cn{Q\k)^H{p,q\k,k), 

(Q\k) gp - cn iq\k) — ^ 

and the natural anzatz for the Hamiltonian (suggested by ( 7.140D ) is 



(7.155) 



H 



(p, q\k, k) = a{q\k, k) ■ cn (^p (3{q\k, k) \ 7(g|fc, k)^ = acn(p/3|7), 



H{P, Q\k, k) = a(g|fc, k) ■ cn [p f3{Q\k, k) \ j{Q\k, k)) = Scn{PP\^) 



(7.156) 



Substituting these anzatz into ( 7.155| ) and making use of (7.136), one can arrive, after some 
calculations il, alF] 



a^{q\k,k) = a^iq\k) = 1 - 



sn^{q\k) ' 



[3\q\k, k) = k'^ + ea^{q\k), 7^(g|A:, k) = 



Pa'{q\k) 

(7.157) 



with some constant v and finally the double-elliptic duality becomes 



H{P, Q\k, k) = cn{q\k) = Q!(Q|fc)cn P + k'^a^{Q\k) 



ka{Q\k) 



k'^ + k^a^{Q\k) ^ 



(7.158) 



H{p, q\k, k) = cn((5|fc) = a(g|fc)cn /c'^ + k'^a'^[q\k) 



ka{q\k) 



k'^ + k^a^{q\k)^ 



(7.159) 



We shall now consider various limiting cases arising from these and show that the double-elliptic 
Hamiltonian ( 7.159| ) contains the entire Ruijsenaars-Calogero and Toda family as its limiting cases, 
as desired. (Of course we have restricted ourselves to the SU{2) members of this family.) 

In order to convert the elliptic dependence of the momentum p into the trigonometric one, the 
corresponding "bare" modulus k should vanish: fc ^ 0, fc'^ = 1 — fc^ ^ 1 (while k can be kept finite). 
Then, since cn(a;|fc = 0) = cosh a;, 



fjdeiii^p^ g) — > a{q) coshp = H^{p, q) (7.160) 

with the same 

Thus we obtain the SU (2) elliptic Ruijsenaars Hamiltonian.^ The trigonometric and rational Ruijse- 
naars as well as all of the Calogero and Toda systems are obtained through further limiting procedures 

^^For ease of expression hereafter we suppress the dependence of a,/3, 7 on k and k in what follows using C((q) for 
a{q\k,k) and a{Q) for a{Q\k,k) etc. 
^2 Indeed, 

F-iq) = cH^\k) (p(e-) - = ( .,.,/^.„. - ..,./^_,„J = ( 1 - ^'^'^^"^ 

where q = 2ujq^/ei3 and 2u^ = sn^(e|fc). 



sn2(^ei3e|fc) sn^ ei3q\k) I sn^{q\k) \ sn'^{q\k) 
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in the standard way. 
The other hmit k 



(with k finite) gives a{q\k) atrig{q) = 1 — 



cosh q 



and 



at„giq) ■ cn ( p^/ fc'2 + k^a^^ig{q) 



katrig(q) 



H^{p,q). 



(7.162) 



This is the elhptic-trigononietric model, dual to the conventional elliptic Ruijsenaars (i.e. the tri- 
gonometric-elliptic) system. In the further limit of small q this degenerates into the elliptic-rational 
model with atrigiq) — > otrat{q) = 1 — ^^-j which is dual to the conventional elliptic Calogero (i.e. the 
rational-elliptic) system, analyzed in some detail above. 

Our approach has been based on choosing appropriate functions f{q) and F{Q) and implementing 
duality. Other choices of functions associated with alternative free Hamiltonians may be possible. 
Instead of the duality relations ( [7.155 ) one could consider those based on ho{p) = sn(p|fc) instead of 
cn(p|fc). With this choice one gets somewhat simpler expressions for /3s and 7^: 



/3s = 1, 7s((j|fc,A;) = fcas(g|fc), as{q\k) = I - 



2v' 



cv?{q\k) 



(7.163) 



and the final Hamiltonian is now 



Hs{p, q\k, k) = aMk) ■ sn{p\kas{q\k)). (7.164) 

Although this Hamiltonian is somewhat simpler than our earlier choice, the limits involved in obtaining 
the Ruijsenaars-Calogero-Toda reductions are somewhat more involved, and that is why we chose to 



present the Hamiltonian (7.159) first. 

One might further try other elliptic functions for /lo (p) • Every solution we have obtained by making 
a different anzatz has been related to our solution ( 7.159| ) via modular transformations of the four 
moduli k, k, k^ff = 7 and k^ff = 7. 

Another approach to the 2-particle double elliptic systems was proposed in [[178| . It is suggested 
there to start with the four dimensional manifold, namely, the elliptically fibered K3 manifold. This 
manifold is the space where the separated variables live. Since it is a phase space manifold, it is 
provided with nontrivial Poisson brackets. They follow from the natural Poisson structure on the 
KB manifold. The Hamiltonian can be chosen as a linear function on the base manifold while the 
dual Hamiltonian as a linear function on the fiber. Using explicit formulas for the Poisson brackets of 
coordinates on the elliptically fibered K3 manifold, in principle, one may obtain explicit Hamiltonians. 
We are planning to report on these calculations as well as on equivalence of the two approaches in a 



separate pubhcation |188| 



7.7 Duality, r- functions and Hamiltonian reduction 

Now one should extend the results of the two-particle case to the generic number of particles. Un- 
fortunately, it is still unknown how to get very manifest formulas. Here we describe a method that 
ultimately allows one to construct commuting Hamiltonians of the iV-particle Dell system, but very 
explicit expressions will be obtained only for first several terms in a perturbation theory. 

However, in order to explain the idea, we need some preliminary work. In fact, we start from 



the old observation 189, [72, 190 that, if one study the elliptic (trigonometric, rational) solution 



to the KP hierarchy, the dynamics (dependence on higher times of the hierarchy) of zeroes of the 
corresponding r-function w.r.t. the first time is governed by the eUiptic (trigonometric, rational) 
Calogero system. Moreover, if one starts with the two-dimensional Toda lattice hierarchy and studies 
the zeroes of the r-function w.r.t. the zero (discrete) time, the corresponding system governing zeroes 



is the Ruijsenaars system [190 



Let us illustrate this with the simplest example of the trigonometric solution of the KP hierarchy. 
This is, in fact, the A^-solitonic solution, the corresponding r-function being 

T{{tk}\{l^^,fI^,X,})^detNxN - 6^-= ^ ) , L,, EE ^' (7.165) 

Here Vi, fii and Xi are the soliton parameters. Let us impose the constraint Vi = Hi + e. Then, 
the standard Hamiltonian structure of the hierarchy on the soliton solution gives rise to the Poisson 
brackets (in |ll3|, Q similar calculation was done for the sine-Gordon case) 



{X^„X,} = — {M.,/ij} = 0, {X„^,,} = XA, (7.166) 



Note that, upon identification /ij ~ qi, Xi = e^* Jl/^^i ( {q -qiY^ ~ «^)' matrix Lij in ( 7.165 ) 



becomes the rational Ruijsenaars Lax operator, with the proper symplectic structure {pi,qj} = Sij. 



At the same time, (7.165) can be rewritten in the form 



T{{tk}\ {i^^,ti.,,X,}) ^ det^xAT (^e^"% - i.^e^Mi "^^^ (7.167) 

where x = ti. This determinant is the generating function of the rational Ruijsenaars Hamiltonians. 
On the other hand, the N zeroes of the determinant as the function of ex are just logarithms of 



the eigenvalues of the Lax operator. These zeroes are governed, as we know from |73, 189| by the 



trigonometric Calogero system, and are nothing but the N particle coordinates in the Calogero system: 

N 

r{{tk}\ W^,^i^,X,}) ^Y.^'^'^^k - X{s\n{tx-q,) (7.168) 

k i 

Therefore, the eigenvalues of the (rational Ruijsenaars) Lax operator are the exponentials e''' of 
coordinates in the (trigonometric Calogero) dual system, while the r-function simultaneously is the 
generating function of the Hamiltonians in one integrable systems and the function of coordinates in 
the dual system. This is exactly the form of duality, as it was first realized by S. Ruijsenaars |p2|. Its 



relation with r-functions was first observed by D.Bernard and O.Babelon |Q (see also [191|). 

Similarly, one can consider the A^-solitonic solution in the two-dimensional Toda lattice hierarchy. 
Then, 



rn{{tk}\ {l^^.^i^,X,}) = dct^x^ (^<5y - L,, gS. (^.-^^ )+*^ (^T'' ) ^ ^ ^ _A_ 

(7.169) 

This time one should make a reduction i/j = e'^^i. Then, Lij becomes the Lax operator of the 
trigonometric Ruijsenaars, and the zeroes w.r.t. to tn are governed by the same, trigonometric 
Ruijsenaars system. This is another check that this system is really self-dual. 

1 no 



Thus constructed duahty can be also interpreted in terms of the Hamiltonian reduction [124|. 
In doing so, one starts with the moment map which typically is a constraint for two matrices or 
two matrix-valued functions. In order to solve it, one diagonalizes one of the matrices, its diagonal 
elements ultimately being functions (exponentials) of the coordinates, while the other matrix gives 
the Lax operator (its traces are the Hamiltonians) . Now what one needs in order to construct the 
dual system is just to diagonalize the second matrix in order to get coordinates in the dual system, 
while the first matrix will provide one with the new Lax operator (its traces are the Hamiltonians in 
the dual system). 

As an illustrative example, let us consider a Hamiltonian reduction of the trigonometric Calogero- 
Sutherland system which differs from that discussed in s.2.1 Now one starts from the free 

Hamiltonian system H = tiA^ given on the phase space {A, B) of two N x N matrices, with the 
symplectic form ti{SA A B~^5B) (i.e. A lies in the gl{N) algebra, while B lies in the GL{N) group). 
Then, one can impose the constraint 

A- BAB-^ = P) (7.170) 

where v is & constant, I is the unit matrix and P is a matrix of rank 1. Now one has to make the 
Hamiltonian reduction with this constraint. First of all, one solves the constraint diagonalizing matrix 
B and using the gauge freedom of the system to transform P into the matrix with all entries unit. It 
gives^ Aij — piSij + (1 — where bi are the diagonal elements of B. This is nothing but the 

Lax operator of the trigonometric Calogero-Sutherland model. Traces of powers of this Lax operators 
give the Calogero Hamiltonians. 

One can also diagonalize the other matrix, A. Then, ( 7.170| ) can be rewritten in the form 



AB - BA = i^{B - P) (7.171) 

where P is another rank 1 matrix which is can be gauged out to the matrix with all entries depending 
only on the number of row (let denote them Xi/v). Solving then this constraint, one obtains B = 
^ j^' _^ , where are the diagonal elements of the matrix A. This is nothing but the Lax operator 
of the rational Ruijsenaars-Schneider system, traces of its powers give a set of the dual Hamiltonians. 



Certainly, this picture |124| nicely explains the Ruijsenaars observation of duality ||6^ and is just 
the T-function approach told in different words. 

The same scheme, certainly, works in for the trigonometric Ruijsenaars too, where one can use a 
similar procedure in order to show that diagonalizing either of the matrices A and B, one arrives at 
the same system, trigonometric Ruijsenaars model, checking again that it is self-dual. 

In fact, similarly to the trigonometric Calogero system, there are several equivalent ways to obtain 
the trigonometric Ruijsenaars model. Say, instead of the Hamiltonian reduction described in section 
2, one can make the Hamiltonian reduction of the gauge theory given on a two-torus T"^ = x 
[ |178| . The other way of doing is to consider the Poisson reduction |^ . Namely, one introduces a set of 
commuting functions on the space of graph connection on a graph corresponding to a moduli space of 
flat connections on a torus with one hole. Being reduced to a particular symplectic leaf of the moduli 



space of flat connections on the torus in the symplectic structure |192[ , this set of functions turns out 
to be a full set of commuting Hamiltonians. This picture is far more geometrical and the self-duality of 
trigonometric Ruijsenaars system gets within this framework rather transparent geometrical meaning. 



Since (7.170) does not fix the diagonal elements of A, they are parametrized by arbitrary numbers Pi 



That is, this Poisson reduction physicaUy is nothing but the three dimensional Chern-Simons theory 
on X (the hue gives the temporal variable) with the insertion of an appropriate temporal Wilson 
line and spatial Wilson loop. It is the freedom to choose the place of the latter one that is in charge of 
the duality. The group of (self-) dualities of this model is therefore generated by the (modular) group 
acting on the non-trivial torus cycles. 



7.8 Dual Hamiltonians for many-body systems 

Unfortunately, the scheme described in the previous subsection does not that immediately applied to 
elliptic models with arbitrary number of particles. The reason is that even having the dual Lax oper- 
ator calculated, one typically meets in elliptic cases a problem of constructing the dual Hamiltonians, 
i.e. invariant combinations involving the Lax operator. 

Instead, in order to obtain dual Hamiltonians, here we use the other approach discussed above, that 
dealing with zeroes of the r- functions | ]9^ , . This approach gives no very explicit form of the dual 
Hamiltonians. Therefore, we construct later a kind of perturbative procedure for the Hamiltonians 
which allows one to get them absolutely explicitly term by term ||9^ . 

Consider SU [N) {N — g + 1) system. The whole construction is based on the fact that the spectral 
curve of the original integrable system (Toda chain, Calogero, Ruijsenaars or the most interesting 
double elliptic system) has a period matrix Tij(a), i, j = 1, . . . , N with the special property: 

N 

ET„(a)=T, Vi (7.172) 

where r does not depend on a. As a corollary, the genus- theta-function is naturally decomposed 
into a linear combination of gcnus-(; theta-functions: 




N \ N-l 



(7.173) 



riieZ \ i,j = l i=l I k=0 



where 



&k^&i'\p,\f,,)^ J2 exp 




(7.174) 



and Pi = C+p^, J2^^iPt = 0; T,j = + ^^^i = 0, V«. 

Now we again use the argument that zeroes of the KP (Toda) r-function (i.e. essentially the 
Riemannian theta-function) , associated with the spectral curve ( [7.191 ) are nothing but the coordinates 
Qi of the original (Calogero, Ruijsenaars) integrable system. In more detail, due to the property ( 7.172 ), 
Q(^^p\T) as a function of — jj X^i^i Pi elliptic function on the torus (1, t) and, therefore, can 
be decomposed into an iV-fold product of the gcnus-onc thcta-functions. Remarkably, their arguments 
are just C ^ Qf- 



N 



e(^) ip\T) = cip, T,T)Y[0iC- ip, T)\t) (7. 175) 

(In the case of the Toda chain when r — )■ ioo this "sum rule" is implied by the standard expression for 
the individual e'^' through the KP r-function.) Since one can prove that qi form a Poisson-commuting 



set of variables with respect to the symplectic structure (7.176) [Q, this observation indirectly justifies 
the claim |9^, ^ that all the ratios 6fc/9; (in order to cancel the non-elliptic factor c(p, T, r)) are 
Poisson-commuting with respect to the Seiberg-Witten symplectic structure 



JV 



^dpf'^'Adai (7.176) 

1=1 

(where oa? = -ai - ... - a^^i, i.e. Ylf=i = 0)- 
It can be equivalently written as 

e,{ej, Qk} + OjiOk, e,} + efe{e„ e,} = o wij, k 

or 

{ioge„ioge,} = |iog|j,iogefc| Vz, 

The Hamiltonians of the dual integrable system can be chosen in the form Hk = O^/Oo, k = 1, . . . , g. 
However, these dual Hamiltonians are not quite manifest, since they depend on the period matrix of 
the original system expressed in terms of its action variables. Nevertheless, what is important, one can 
work with the ^-functions as with series and, using the instantonic expansion for the period matrix 
Tij, construct Hamiltonians term by term in the series. 

7.9 Dual Hamiltonians perturbatively 

Let us see how it really works. We start with the simplest case of the Seiberg-Witten family, the 
periodic Toda chain. 

Perturbative approximation. First we consider the perturbative approximation which is nothing 



but the open Toda chain, see ( 5.16 ) and ( 5.17 ). This is the case of the perturbative Ad pure N = 2 



SYM theory with the prepotential 



1 ^ 

^pertia) = ^Y1 4log"y (7-178) 

Then, the period matrix is singular and only finite number of terms survives in the series for the 
theta-function: 

eW(p|T) = J2 (7.179) 

(7.180) 



fe=0 



Here 



i,[i]=kiei jg7 
zif(a) = e-^^'=^ (7.181) 



and / are all possible partitions of N indices into the sets of fc = [/] and N — k = [I] elements. 
Parameter Aqcd becomes significant only when the system is deformed: either non-perturbatively 
or to more complex systems of the Calogero-Ruijsenaars-double-elliptic family. The corresponding 



r-function Q'-^\p\T) , eq.(7.179), describes an iV-soliton solution to the KP hierarchy and is equal to 



the determinant ( [7.165 ). The Hamiltonians \ eq.( 7.18C ) dual to the open Toda chain, are those 



of the degenerated rational Ruijsenaars system, and they are well-known to Poisson-commute with 



respect to the relevant Seiberg-Witten symplectic structure (7.176). This is not surprising since the 
open Toda chain is obtained by a degeneration from the trigonometric Calogero system (s.3.3) which 
is dual to the rational Ruijsenaars. 

The same construction for other perturbative Seiberg-Witten systems ends up with the Hamilto- 
nians of more sophisticated systems. 

For the spectral curves (5.71) of the trigonometric Calogero system (perturbative Ad J\f — 4 SYM 



with SUSY softly broken down to A/" = 2 by the adjoint mass M) the Poisson-commuting (with respect 
to the same ( 7.176| )), whose perturbative prepotential is given by (|7^), Hamiltonians are given 



by (|7.180| ) with 



M2 



(7.182) 



i.e. are the Hamiltonians of the rational Ruijsenaars system which are really dual to the trigonometric 
Calogero model. 



For the spectral curve ( 5.67 ) of the trigonometric Ruijsenaars system (perturbative 5dAf = 2 SYM 
compactified on a circle with an e twist as the boundary conditions) whose perturbative prepotential 
is given by ( f3.7C| ), the Hamiltonians are given by (7.180) with 



^ sinh(ay -f e) sinh(ay — e) 



sinh a, 



(7.183) 



i.e. are the Hamiltonians of the trigonometric Ruijsenaars system, which is, indeed, self-dual. 

Finally, for the perturbative limit of the most interesting self-dual double-elliptic system (the 
explicit form of its spectral curves is yet unknown except the case of two particles, see the next 
subsection) the relevant Hamiltonians are those of the elliptic Ruijsenaars system, given by the same 



(7.180) with 



252 



-i] v-y sn£(ay) 



(7.184) 



where r is the modulus of the second torus associated with the double elliptic system and are aij 



rescaled, [109]. 

In all these examples = 1, the theta-functions 0*^^^ are singular and given by determinant 
(solitonic) formulas with finite number of items (only terms with Ui = 0, 1 survive in the series 
expansion of the theta- function), and Poisson-commutativity of arising Hamiltonians is analytically 
checked within the theory of Ruijsenaars integrable systems. 



Beyond the perturbative limit. Beyond the perturbative limit, the analytical evaluation of O'^^) 
becomes less straightforward. 



1 na 



The non-perturbative deformation (3.6) of the curve (5.16), is associated with somewhat sophisti- 
cated prepotential of the periodic Toda chain, 



The period matrix is 



k=l 



(7.185) 



-2Wy +log^ + ^A' 

k=l 

and r in this case can be removed by the rescahng of A. Then, 



,9^.FW(a) 
daiduj 



(7.186) 



N-l 



JV-1 



A:=0 



N 



N 



_i_ \^ „27ri(pi+pj-pfc-p,) J74 jr4 ;r4 jr4 TT jz. 7-. 



N 



(7.187) 



with 



1-3 



daidaj 



(7.1E 



The first few corrections Z''^^ to the prepotential are exphcitly known in the Toda-chain case [193 , 
for example, 



N 



2iTr ^-"-f 

i—1 k^i 



(7.189) 



The coefficients Qk in ( 7.187j ) are expanded into powers of (A/a)^^ and the leading (zeroth-order) 



terms are exactly the pertm'bative expressions (7.18C ). Thus, the degenerated Ruijscnaars Hamilto- 



nians (7.18C) are just the perturbative approximations to the Hk — Ofc/0o - the full Hamiltonians of 
the intcgrable system dual to the Toda chain. 



1 Ci^-r 



The same expansion can be constructed for the other, eUiptic systems. Note that the period matrix 
in these cases is expanded into powers of q 



— Anil Uij 

OQiOaj 



^ daidaj 

k—l 



(7.190) 



and the dimensional constant /i depends on the system. E.g., in the Calogero model /i = iM etc. The 
limit to the Toda system corresponds to g — + 0, g/i^^^fixed. 

Note that the first (instanton-gas) correction to the prepotential is indeed known to be of the same 



universal form (7.189) not only for the Toda chain, but also for the Calogero system [154 . For the 



Ruijsenaars and double-elliptic systems eq.(7.189) is not yet available in the literature. 

In ]9^ the perturbative expansions for the Hamiltonians were tested for Poisson commutativity, 
up to few first terms and small number of particles N . In particular, it was checked that ( 7.189| ) is 
true for both the Ruijsenaars and Dell systems. 

The main lesson we could get from this consideration is that if the universal expressions like ( [7.1891 ) 



in terms of perturbative Z^i^\ the same for all the systems, will be found for higher corrections to the 
prepotential^^, this will immediately give an explicit (although not the most appealing) construction 
of the Hamiltonians dual to the Calogero and Ruijsenaars models and, especially important, the 
self-dual double-elliptic system which was explicitly constructed above only for N = 2. 

Let us note now that the SU{2) Hamiltonian ( 7.159| ) has a rather nice form in terms of the p,q 
variables or their duals P, Q. However, this our general SU{N) construction of the Hamiltonians 
as ratios of genus g 6'-functions uses another kind of canonical variables - angle-action variables 
Qi. The flat moduli play the central role in Seiberg-Witten theory theory, while Qi are rather 
generalizations of the algebraic moduli. 

It is an interesting open problem to express the Hamiltonians in terms of Pi and Qi, perhaps, they 
can acquire a more transparent form, like it happens for SU{2). Another option is to switch to the 



"separated" variables Vi, Qi |179[ such that £^ VidQi = aiSij (while generically Oi = 
Generically, they are different from Pi, Qi but coincide with P, Q in the case of SU{2). 



7.10 Dell systems: 6d adjoint matter 

Thus, we have constructed the Dell system and a set of dualities that acts on the Toda-Calogero- 
Ruijsenaars-Dell family. They are all put together in Fig. 2. 

To avoid possible confusion, the recurrent relation |157, 149, 194 for the Toda-chain prepotential, 



-log©o 



p=0 



does not immediately provide such universal expressions. Already for Z^^^ this relation gives 



2 

i<j k^i.j 



(0) 
jk 



which coincides with (r.l89) for the Toda-chain \ eq.(7.181), but is not true (in variance with (7.189)) for Calogero 
\ eg. ( 7.18^ ). Meanwhile, the recurrent relations of ref!" ^95| for the Calogero system provide more promising expan- 



sion 



All this emphasizes the need to study extension of the Whitham theory and the WDVV equations from the Toda 
chains to the Calogero and Ruijsenaars systems. 
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Figure 2: Action of the coordinate-momentum duality on the Calogero-Ruijsenaars-Dell family. 
Hooked arrows mark self-dual systems. The duality leaves the coupling constant v intact. 



Since the theories of this family describes low energy limits of different SUSY gauge theories, this 
table describes the dualities between different gauge theories too. Say, perturbative limit of the 4(i 



theory (trigonometric Calogero) is dual to a special degeneration 1 148 of the perturbative limit of the 
5cZ theory (rational Ruijsenaars), while full, non-perturbative 5d theory (elliptic Ruijsenaars) is dual 
to the perturbative limit of the 6d theory (Dell system) . 

To conclude our discussion of dualities and Dell systems, we calculate the perturbative prepotential 



of the 2-particle Dell system (7.159) and show that it coincides with what is to be expected for the 6d 



theory with adjoint matter, supporting the identification of the Dell system with the 6d gauge theory. 

Let us start from the full Dell model for 2-particles. Then, it is given by the Hamiltonian ( 7.159| ) 
parametrized by two independent (momentum and coordinate) elliptic curves with elliptic moduli k 
and k. As usual, for two particles we can construct the full spectral curve from the Hamiltonian ([p7|. 



see also ss.3.2, s.3.7 and (|7.127| )) 



HiC,C\k,k)^u (=cn(Q|fc)) 
It is characterized by the effective elliptic moduli 



(7.191) 



ka{q\k) 



^eff 



ka{q\k) 



(7.192) 



where the functions a and (3 are manifestly given in (7.157). Coordinate- momentum duality inter- 
changes k ^ k, keff ^ keff (and q,p ^ Q,P), while in general SU{N) case, the model describes an 
interplay between the four tori: the two bare elliptic curves and two effective Jacobians of complex 
dimension g = N — 1. 

From our previous discussions it should be clear that the generating differential in 6d theories is 
of the form dS — C^Cj with ^ living on the coordinate torus and C ~ on the momentum torus. 

Now we calculate the leading order of the prepotential expansion in powers of k when the bare 
spectral torus degenerates into sphere. In the forthcoming calculation we closely follow the line of 
s.5.6. 

When k —^ 0, sn(q\k) degenerates into the ordinary sine. For further convenience, we shall param- 
eterize the coupling constant 21/"^ = sn^(e|A;). Now the spectral curve ( [7.191 ) acquires the form 



«(0 



sn^(e|fc) u 



cn C/3 



Here the variable ^ lives in the cylinder produced after degenerating the bare coordinate torus. So 
does the variable x ^ 1/ sin^ ^. Note that the A-period of the dressed torus shrinks on the sphere to 
a contour around a; = 0. Similarly, i3-period can be taken as a contour passing from x = to a; = 1 
and back. 

The next step is to calculate variation of the generating differential dS = r/d^ w.r.t. the modulus 
u in order to obtain a holomorphic differential: 

dv = (-isn(e|A:) + k'^uA — (7.194) 
^ ^ x^y (x ~ 1){U^ ~ x) 

where t/^ = V.",!^ . Since 

da r ddS__ f 1 

we deduce that u = cn{a\k) and U = . The ratio of the B- and A-periods of dv gives the period 

matrix 

T = - / = = -— lim log- + —log- —r 7.196 

Jo Xy^{x-l)iU^-x) ITTK-^OV '^4/ ITT '^l-[/2 

where k is a small- a; cut-off. The U dependent part of this integral is finite and can be considered 
as the "true" perturbative correction, while the divergent part just renormalizes the bare "classical" 
coupling constant r, i.e. classical part of the prepotential (see s.5.6 for further details). Therefore, 
the perturbative period matrix is finally 

„ i, sn"^ (elk) — sn"^ (a\k) i, Oiia + e)6iia — e) 

Tf^mte = -log ^ ^ -I- const -\og^ -2^-^^^ ^ 7.197 

TT sn^(a|fc) TT "iW 



and the perturbative prepotential is the elliptic tri-logarithm. Note that it is again of the form (5.7C) 



and ( 5.69| ) but with the function /'^^ instead of Z*^^^ and /^^^ Remarkably, it lives on the bare 



momentum torus, while the modulus of the perturbative curve ( [7. 193 ) is the dressed one 



7.11 Quantum duality 

Our discussion so far concerned the classical systems. However duality is expected to be a powerful 
tool to deal with quantum problems too. This point has been first recognized by Ruijsenaars in his 
early papers. The very idea is that Hamiltonians and their duals have common wave-functions and 
this fact could be exploited to consider one or another set of eigenvalue problems, with respect to 
different variables. Since the wave functions in the systems under consideration can be realized in 
group terms, one could expect that quantum duality can be formulated purely in group terms too. 



Recent activity in this direction supports such expectations \19t\. 

Here we work out a couple of simple quantum 2-particle examples of dual systems in order to 
illustrate the general picture. Let us start with the most standard oscillator example with the Hamil- 
tonian 



9^ uP-q^ 

Its normalized eigenfunctions are 



H = -l/2— + ^ (7.198) 



Hipn = w(n + 1/2)V„ 

(7.199) 

where Hn{^) are the Hermite polynomials -ffn(^) = e^^(— 9j)"'e~^^. Using this representation for the 
wave-function one can easily obtain a recurrence relation 

y/n+ltpn+i{q) + Vntpn-i{x) = \/2uqtpn{x) (7.200) 
It means that tpnil) is the eigenfunction of the difference operator 

Hd = T+sfR+sfRT-, T± = e^^ (7.201) 
acting on the subscript n, which is nothing but the quantized dual Hamiltonian. 

Less trivial example whose classical version we discussed throughout this whole section is the 
quantum (2-particle) Hamiltonian of the trigonometric Calogero model 



2d^^ 2sm^{q) 

Its normalized eigenfunctions are 



n2 



V„(a) = s^n^{q)^nj^^^^Tf-_\ (cos(,)) (7.203) 

For the sake of simplicity, wc take half-integer v and n. One can make the change v —v — 1 and 
get another eigenfunction with the same eigenvalue. Using the fact that the generating function for 

nO's 

fr'o V'l - 2xy + t/2 

satisfies the two obvious equations 

(x - y)dxZ = ydyZ 

(7.205) 

{l-2xy + y'')dyZ={x-y)Z 

which implies 

{ydy - m)d^Z ={x- y)d^+^Z (7.206) 
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and, hence, yields 



((1 - 2xy + y^)dy + y- x)dZ'Z = m(l + 2ydy)d'^'^Z (7.207) 
one derives the recurrence relations 



_[^l_m l + m 

"^^^^ - 21+1 "'+i + 2rri"'-i 



cos(g V« = W ? — --r^V'n+i + \ t r^V-n-i 

n(n + 1) y n(n — 1) 



(7.208) 



i.e. Tpn is an eigenfunction of the finite-difference operator acting on the subscript n 

Ho^iq) - cosiq)^p{q) (7.209) 




Hn-T,,ll- 4^ + .1 - (7.2f0) 
" n[n — f ) y ?i(n — 1) 

which is a quantum version of the rational Ruijsenaars model. 

Summarizing, when the system with trigonometric dependence on momentum is quantized, its 
Hamiltonian becomes a finite-difference operator. The wave-functions become the functions of dis- 
crete variables. The origin of this is in the Bohr-Sommerfeld quantization condition. Indeed, since the 
trigonometric dependence on momenta implies that the leaves of polarization arc compact and, more- 
over, non-simply connected, covariantly constant sections of the prequantization connection along the 
polarization fiber generically ceases to exist. It is only for special "quantized" values of the action 
variables that the section exists. In the elliptic case, the quantum dual Hamiltonian is going to be a 
difference operator of infinite order. 



8 Branes as low energy degrees of freedom in supersymmetric 
gauge theories 

8.1 Field theory on D-brane worldvolume 

We begin this section with the description of nonperturbative states in string theory which play the 
key role for constructing the low-energy degrees of freedom in SUSY gauge theories. In particular, 
the phase space of the integrable systems, considered in the previous sections is associated with the 
collective coordinates of branes of different dimensions. Besides it, nonperturbative dualities in string 
theory suggest several different formulations for the same single field theory. 

The central objects introduced into string theory by exact solving the Af = 2 SYM theory are 
the duality (between the strong and week coupling regimes) and the exact spectrum of BPS solitons. 
It appeared that in order to consistently describe duality in string theory, one has to introduce 
solitonic objects of different dimensions, "p- branes" with p space coordinates [Q. In quasi-classical 
approximation, their masses behave as m — g^^, where g is the string coupling constant, and the 
mass formulae for the BPS states, being related to the central extension of the SUSY algebra remains 
exact at quantum level. At the strong coupling regime, the solitons become the lightest objects in the 
theory and define the dynamics in that region. Technically, it is possible to introduce dual variables 
which yield a week coupling theory in terms of the solitonic fields. 

1 1 o 



All this picture is literally inherited from Seiberg-Witten theory. On the other hand, Seiberg- 
Witten theory itself gets a new description in string/brane terms. Indeed, the appearance of Riemann 
surfaces in Seiberg-Witten theory can be explained in terms of branes. It is this surface which is the 
part of the worldvolume of the larger brane. The structure of the brane intersections is fixed by the 
SUSY and the set of configurations yielding the field theory become finite or even single. 

As we already explained in the Introduction, the key point which admits the use of branes for 
obtaining field theories is the U{1) gauge theory on the brane worldvolume. It is known that branes 
are charged with respect to the RR sector of the theory and can serve as a space of endpoints of 
strings. Using this property Witten formulated the new view on the gauge theory as the theory on the 
worldvolume of coinciding branes which carry U{1) factor each. The open strings connecting branes 
give rise to non-Abelian theory, rank of the gauge group being defined by the number of branes. In 
what follows we shall follow this line: SUSY field theories are defined on the brane worldvolume and 
the parameters are introduced via the geometry of the brane configuration. 

Some parameters of TV = 2 SUSY theories (condensates of scalars) coordinatize the Coulomb 
branch of the corresponding vacuum moduli space, while other parameters (bare masses and coupling 
constants) remain external parameters. All these parameters are defined in the brane terms (see the 
Introduction): scalar condensates yield the distance between branes (or determine embedding of the 
single M-theory brane into the target space), coupling constant is just the distance between additional, 
background branes (or also determines geometry of the M-theory brane), while the hypermultiplet 
mass is the coordinate of another brane (often of different dimension) in the proper direction. In what 
follows we shall present examples of particular brane configurations. Coordinates of branes appear to 
be associated with the degrees of freedom of the integrable system governing the low energy action in 
the SUSY gauge theory. 

Let us emphasize that the gauge theories of the Hitchin type responsible for the integrable many- 
body systems arise in the "momentum" space or "the space of fields" and should not be thought of as 
those obtained by the dimensional reduction from the Ad SUSY gauge theory. In fact, these theories 
emerge from different degenerations of the toric geometry. 

It is useful to start with the softly broken JV = 4 theory. In A/" = 1 terms, there are 3 adjoint 
hypermultiplets whose complex scalar components parameterize the embedding of branes into six 
dimensions. The potential of scalars is 

^^^Tr ^ (8.1) 

^ 1=1,2,3 

and and is induced by the kinetic term (as a result of integration over auxiliary fields) . There is also a 
non-trivial superpotential, W = Tr$i [$2, ^3] + -M^Tr($| -|- <!>§), where are A/" = 1 hypermultiplets 
with lowest components 4>i respectively and M is the mass of the hypermultiplet. The minimum of 
the potential energy amounts to the vacuum configuration with non-trivial commutation relations 
between the 3 scalars. In the large mass limit, the fluctuations orthogonal to the algebraic surface in 
the two dimensional complex space 

Tr($^ -I- <Pl) = const, (8.2) 

are frozen. Note that, since one of the scalars remains massless, it jointly with the gauge superfield 
forms one J\f — 2 supcrmultiplet, while the two massive scalars are the components of the M = 2 
adjoint massive hypermultiplet, i.e. the A/" = 4 SUSY is broken down to the N = 2 SUSY. In this 



J\f ~ 2 theory, there are vaUeys for the massless scalar (pi due to the potential term ( ^.iD , i.e. this 
scalar may have a condensate which parameterizes the Coloumb branch. 

There are two equivalent approaches to M = 2 theories in the string theory framework which start 
from the theory in ten dimensions. In the first, geometrical engineering approach, one considers the 
compact Calabi-Yau (CY) threefold near the singularity providing the desired gauge group |5^. In 
the second approach proposed in , the theory on the worldvolume of D4-branes stretched between 
two NS 5-branes is considered. The location of the branes is summarized in the following table 



NS5 
DA 
D6 



1 

+ 
+ 



2 

+ 
+ 



6 7 



+ 



(8.3) 



where plus stands for the extended direction of each brane and minus for the transverse direction. 
The fourbranes are finite in the xq direction stretching between the two NS branes, therefore, the 
long range theory is a four-dimensional theory as we need. While in this type IIA picture the original 
SUSY (32 real supercharges) is broken down to N = 2 in Ad (8 real supercharges) due to the concrete 
brane configuration^, in the IIB picture (related to CY), the original SUSY is broken via geometry 
of the CY manifold. D-branes wrapped around shrinking cycles provide the field theory and break 
half of the supersymmetries indeed. Note that the IIA picture can be lifted to M-theory in lid which 
has the lid SUGRA as its low-energy limit. 



8.2 Branes and degrees of freedom in many-body systems 

Let us start with the description of degrees of freedom for the 4rf SU{Nc) theory without matter 



in IIB/F picture [197 . It emerges as the worldvolume theory on D7-branes wrapped around an 
elliptically fibered K3 . The bare coupling of the gauge theory yields the complex structure of the 
elliptic fiber, therefore, the spectral curve of the Calogero model enters the picture automatically. 

To identify the degrees of freedom using the approach of section 2.3, one has to define the holo- 
morphic one-form (j) as well as the anti-holomorphic vector potential A on the spectral curve. The 
holomorphic form can be identified if one takes into account that the torus is embedded into the K3 
space as a supersymmetric cycle. Therefore the embedding of D7 branes is described in terms of 
the one-form which after all turns out to be the Lax operator for the many-body system. The 
anti-holomorphic vector field A whose eigenvalues are the coordinates of particles in the many-body 
system comes from the Nc branes wrapped around the spectral curve with the worldvolume gauge 
fields. Finally the marked point on the spectral curve with nontrivial monodromy is introduced via 
additional "background" branes localized on the spectral curve. Such "background" branes play the 
role of regulators in field theory. Let us emphasize that in F-theory the 11-th and 12-th dimensions 
just provide the bare coupling constant. 

Since the background brane does not change the interpretation of degrees of freedom let us con- 
sider the simplest case of the free many-body system. The trigonometric Ruijsenaars system, i.e. 
the 5d SUSY theory is the most illustrative, since in this case the brane geometry is in one-to-one 
correspondence with the Hamiltonian reduction procedure of section 2. Let us remind that the system 
originates from the Chcrn-Simons theory given on the product of torus and real (time) line. There is 

^^The presence of D4-branes breaks symmetry down to 16 supercharges and two NS branes breaks it down further to 
8 supercharges. 



the relation for the monodromies of the gauge field around the torus cycles gA,gB 



9AgBgA^9B^ = 1 

In the fi'ee case the matrices can be diagonalized simultaneously and their diagonal elements represent 
the coordinates and momenta of the dynamical system — diag{e^'^^), gs — diag{e^P'). Note that 
from the F-theory point of view we have degenerated the fiber torus to a cylinder. 

Let us compare the geometrical descriptions following from the Chern-Simons Lagrangian and the 
brane configuration. Two coordinates x^^xq are essential for the description of the five dimensional 
theory. The typical scale along 0:5 is R^^, while that along xq is g~^ in the string units [198|. We 
shall assume that both coordinates are compactified onto circles of the corresponding radii. 

The (coordinate and momentum) circles in the Ruijsenaars model are related to the circles in the 
brane configuration via the T-duality transformation along the both directions. Therefore, the radii of 
circles in the Ruijsenaars model are i?5, g^, and after the T-duality transformation the eigenvalues of 
the monodromy matrix transform into the brane positions along the dual circle. There are precisely 
Nc 5-branes in the IIB theory extended along xq and localized at pi along xq. The coordinates 
in the dynamical system are provided by Nc 1-branes localized at qi along xq. The SUSY gauge 
theory is defined on {xq, Xi, X2,X3,X4, Xq), and the coordinates of 5-branes along x^ define the vacuum 
expectation values of the real scalar field. 

To switch on the interaction in the dynamical system, one has to add a massive hypermultiplet in 
adjoint in the field theory language, to introduce a marked point with non-trivial monodromy around 
in the Chern-Simons description or to add the background branes in the string theory approach. 

The new relation for the monodromy matrix is 



gAgsgA^gs^ = 5c, 

where gc corresponds to inserting the orbit at the marked point. After the T-duality transfor- 

mation, a brane configuration with non-trivial flux arises. The color orientation of the fluxes follows 
from the structure of the monodromy matrix gc- 

Another approach to description of fleld theories, for deflniteness let them be 4d theories, is based 
on the consideration of field theory on the worldvolume of M5-brane, the solitonic state in M-theory 
[ pT| . We look for an M5-brane configuration that provides a theory with 8 real supercharges, i.e. in 
Af — 2 theory. This can be achieved by considering the product of and a fixed genus compact 
surface E as the worldvolume of M5-brane. After compactiflcation onto the circle (11-th dimension), 
the M5-brane in the IIA theory looks as a collection of D4-branes located between two NS-branes. 

To get the SU{Nc) theory, one has to introduce Nc D4-branes of flnite size along one of the 
coordinates, say, xq. Then, the plane x^ + ix^ in the perturbative regime is where the momentum 
variables in the many-body system live. Moreover, there are DO-branes on the D4-worldvolume, one 
per each D4-brane, which are ordered along the xg coordinate. In the low-energy limit, the DO-brane 
coordinates are defined 



q,^jA + (b,, (8.4) 

where A fixes their equilibrium locations, while (jfi correspond to small fluctuations and coincide with 
the variables in the Toda system. Let us note that Qi can be also considered as the coordinates in the 
Calogero system, without any ordering. Comparing with the field theory variables one immediately 
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arrives at NcA = ^2 which perfectly fits the suggested ordering of the DO-branes along the xq 
direction. 

A qualitative interpretation of the 2x2 Toda Lax operator is then as follows. The element Ln 
corresponds to the DO-brane confining with the D4-brane so that they have the same {x^ + ix^)- 
coordinate. Since Lax operators are related with fermionic zero modes in a nontrivial background 
field, the size of the Lax matrix indicates the presence of the two dimensional space of zero modes. 
Elements L12, L21 fix the matrix elements of the spectral parameter between the fermionic zero modes 
localized on the neighbour DO-branes. 

Let us note once again that qualitatively in any brane configuration, there are three essential 
elements presented. First, there are the branes that provide the worldvolume for the field theory, 
Nc D7-branes in IIB/F or one M5-branc in IIA/M picture. Second, it is necessary to introduce the 
branes responsible for the variables in the dynamical system, which are zero size instantons in the 
field theory language [ |199[ . Finally, the account of the regulators in field theory corresponds to an 
additional Wilson line in the IIB/F picture or fluxes in the IIA/M description. 

Recently the correspondence between integrable systems and branes has been pushed forward even 
further. Namely, one can start with the Ad theory with one compact dimension [200| and map the 
system to the Hitchin phase space by a chain of dualities. On one hand, such a moduli space coincides 
with the moduli space of the non-commutative instantons on T'* or K3 manifolds |201|. On the other 
hand non-commutative instantons introduced in [ 202 1 allow a deformed ADHM description, moreover, 
the explicit map between the moduli space of non-commutative instantons and the trigonometric 
Calogero model has been elaborated |203|. 



8.3 Branes in five dimensional theory 

While considering 5d gauge theories the role of the brane configuration becomes the most transparent. 
It appears that in five dimensions the brane configuration fixes the structure of the spectral curve of 
the integrable system. We shall consider below the configuration of branes in the IIB theory which is 
related to the toric manifold data providing the compactification of M-theory ^ ^ . 

We start from the geometrical engineering picture ||9^. As a compact manifold we take the 
CY manifold in the vicinity of singularities. Being interested only in this region of parameters, we 
discuss exclusively a local model implying that its embedding into a compact Calabi-Yau is always 
possible. To this end, we consider A„_i ALE spaces fibered over a set of the base spheres, each sphere 
corresponding to an SU (n) gauge theory^ The size of the base sphere V oc is governed by the 
coupling constant of the corresponding SU{n) factor. Kahler moduli of P^'s associated with blowing 
up the ALE space define the Cartan scalars a.; on the Coulomb branch of the moduli space. This 
non-compact Calabi-Yau space is a hypersurface in the holomorphic quotient and can be represented 
as a set of Kahler cones, or, in other words, by the toric polyhedron A. 

For example, one can start from the Af = (2, 2) linear sigma model construction of Calabi-Yau 
[ ^04| and define U{1) charge vectors of matter fields 

Q'^ = (g?, •■■,'??) (8.5) 

From the a-model point of view these vectors define gauge-invariant variables, while from the Calabi- 

^^Throughout this subsection we consider the gauge theory whose gauge group is a product of fc — 1 simple factors, 
typically SU{ni). 
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Yau standing point they act on the homogeneous coordinates in projective space and fulfill r linear 
relations 



EQ>J=0 (8.6) 

J 

between vectors j/j G N = Z'' which define coordinate patches in each Kahler cone. Vertices vj in 
lattice N form a toric polyhedron A ~ {i^j}- For the toric variety to be smooth, A must be convex. 
Thus, having defined the original Calabi-Yau in terms of holomorphic quotients, we now turn to the 
solution of the model via mirror symmetry. 

As it was explained in Q , Kahler moduli of the original Calabi-Yau specify Coulomb moduli 
of the low-energy effective theory. Under mirror transform they get mapped to the complex structure 
of the mirror manifold, which is just Riemann surface E having the same moduli space. There is a 
profound reason in this symmetry, because the Riemann surface is the spectral curve of the associated 
integrable system, XXZ spin chain, section 3.5. 

In terms of charge vectors this mirror manifold is 

Ilyf-^ (8.7) 

J 

For the models we consider here it is just a spectral curve of underlying integrable system. 

Let us now turn to the description of 5d theories in terms of brane configurations. One can derive 
the brane configuration using T-duality along X4 from the IIA picture [pl|. Thus, one obtains the IIB 



picture |198] with the brane configuration |198 



NS5 
Db 



123456789 

+ + + + + -- -- (8.8) 

+ + + + - + -- - 



where the field theory is realized on the D5-brane worldvolume. The holomorphic solution to the 2d 
Laplace equation in the space {xa, x^), As = (s = xg + ixio): s = log(x4 -I- ix^) transforms into the 
solution to the Id Laplace equation s(= a;^) — \\x^\ — ^sign(a;5). If one demands for SUSY to be 
broken up to 1/4, a restriction on the brane intersection arises. In particular, when charged n D5- 
brane meet charged k NS5-branes they must merge into (/c, n) bound state with the slope tanc/) = — 



on (xsjXg) plane [198| 



As usual we divide moduli of brane configurations into two groups. The first group of parameters, 
external parameters (e.g. coupling constants, bare quark masses) define the five-dimensional theory, 
while the rest are dynamical moduli. In brane terms, changing the parameters of the first group one 
changes the asymptotic (i.e. infinite in all their non-zero dimensions) branesQ, while changing the 
parameters of the second group does not. 

This is a good point to make a contact with the geometrical engineering approach In 
that language each link in the toric polyhedron A is associated with a 1-cycle on the corresponding 
Riemann surface E. Periods along 1-cycles in the compact part of E (i.e. internal links in A) define 
the dynamical moduli whereas bare parameters arise from 1-cycles that wrap the non-compact part 
of E (i.e. related to the links on the boundary of A). One can go further and note that the number 
of such external links in A corresponding to some low-energy effective field theory coincide with the 
number of infinite fivebranes in the brane construction of the same theory. 

^^In our case, NS5 branes. 



The reason for this is that the brane picture and the toric polyhedron corresponding to a five- 
dimensional field theory are intimately related to each other. To establish other connections let us visu- 
alize brane diagrams in other terms. Any such configuration depicts the kinematics of two-dimensional 
(zs, xq) free particles. The charge conservation in the vertices where branes merge together is nothing 
but the momentum conservation law for free particles. The integer- valued momentum p^^'^') = {k,n) 



of particle is just the SL(2,'^) charge of the brane. The vertical branes in terms of [198| carry the 
(1,0) charge (momentum) vector, and horizontal ones - (0,1). Notice that the number of vertical 
(NS5) branes specifies the number of simple gauge factors whereas the number of horizontal (D5) 
branes defines the rank of each SU{ni) factor. 

For instance, for SU{n)'^ gauge theory there are fc + 1 NS-branes and n D-branes in each stack. 
And, similarly, the number of vertical links in the toric polyhedron, i.e. the height of A equals fc -I- 1 
and the number of horizontal links, i.e. its horizontal size is n. At generic point on the moduli space, 
the vertices of A lie on a two-dimensional integer lattice. In these conventions we can endow each link 
in A with an integer-valued two-dimensional vector {ks,ns). 

One can exploit the fact that type IIB brane configurations for five-dimensional field theories 
are in one-to-one correspondence with toric polyhedrons of Calabi-Yau M-theory compactifications 
describing these field theories. The exact relation between the two is provided by identification of the 
{ks,ns) charge (momentum) vector with a link vector in A. This powerful tool allows one to solve 
any five-dimensional Y[ SU (n) gauge theory with arbitrary matter. 

The toric polyhedron A for the case of pure gauge SU (n) theory can be found in . Adding a 
massive hypermultiplet, we break one link on the boundary of A into two pieces and add one more in- 
ternal link. Therefore, the number of internal and external links representing massive hypermultiplets 
increases by one. Within the brane approach, the picture is exactly the same: adding a fundamental 
matter hypermultiplet to the theory one attaches a new semi-infinite brane. It can be done in four 
different ways: between the pairs of the four semi-infinite branes already presented in the configura- 
tion.. The same possibilities can be found in breaking one of the four external links of the polyhedron 
into two parts. For the Yl'^^i SU{ni) case, one has to take a fibration of iii spheres over each of k 
base spheres. The total number of external links is 2{k -\- 1) as the number of semi-infinite branes in 
the brane configuration. The number of finite branes is X^iLil^"-* — 1) — 1 which matches exactly the 
number of internal links in A. The relation between toric diagrams and brane configurations has been 



also considered in | 205 | 



8.4 Branes and equations of motion in many-body systems 

Our next step in clarification of the role of brane moduli as degrees of freedom in integrable systems 



is to interpret the equations of motion in the brane terms [197|. To this end, let us relate the brane 
configuration with the monopolc moduli space and identify the equations of motion with the Nahm 
equations defining the structure of the moduli space. 

DO-branes localized on D4-branes can be considered as the degrees of freedom corresponding to 



the BPS monopoles. They can be related to other D4-branes by open strings. Following |206|, the 
theory on the string worldvolume can be obtained by the dimensional reduction from the lOd SYM 
theory. Moreover, there is one-to-one correspondence between the equation defining the supersym- 
metric vacuum state in cr-model and the Nahm equation for monopoles. Now, in order to get the 
equations of motion in the Toda system, one should take into account the connection between the 
spectral curve in the Toda system and the spectral curve amounted from the Nahm equations for the 



charge Nc cychc monopole [207]. The spectral curve for the cyclic configuration can be derived from 
the generic monopole spectral curve 



A" + A"-iaiH + .... + a„(u;) = (8.9) 

where ai{w) is a polynomial of degree 2i. Then, assuming the center of mass is fixed at the origin, 
the overall U{1) phase is unit and imposing the (order Nc) cyclic group condition, one arrives at the 
Toda curve. 

One can also connect the Lax operator for the Toda system with the Nahm matrices Ti known in 



the monopole theory |20S 



T{w) = T+ + 2Taw + T^w^ (8.10) 

where T± are respectively the (strictly) upper and low triangular parts of the Toda Lax operator ( p.4[ ), 
while To is its diagonal partQ. Adopting the notation T± = Ti ± 1X2, one can identify the Toda Lax 
equation (i.e. the equations of motion) 

^ = .,,,[T(u;),M] (8.11) 
where M = T+w^^ — T-w, with the Nahm equations 

^^2te.,k[T„Tk] (8.12) 

To develop a more geometrical picture, it is useful to consider fermions in the background monopole 
field 



{aD + (j)-t)^ = 0. (8.13) 

where t is an external parameter. The equation has k linear independent solutions for the magnetic 
charge k. It is the variable t which becomes the time variable in the Toda system and simultaneously 
plays the role of the space variable in the brane interpretation. The Nahm matrices in terms of the 
fermionic zero modes are 

J xj^+^d^x; j = l,2,3 (8.14) 

These expressions have the structure of the Berry connection if the momentum space is considered as 
the parameter space. There exists also an inverse construction which allows one to derive the multi- 
monopole fields starting from the auxiliary one-dimensional problem for the fermions in the external 
Nahm connection 

d ^ ^ ^ ^+ N 
[i— — h i(7x + aT )v = 

/ . (8.15) 

Aj = / v+djvdt = / tv+vdt 

where the superscript plus of T means Hermitian conjugation. The limits of the integration are fixed 
by the asymptotics of the scalar field 0. Note that the fermions at a; = are identified with the 
spectral fermions of the Toda (or, more generally Nahm) system. 

^*In monopole terms, the variable w is the coordinate on the twister manifold ([^P^. 



One may also consider the generalization of the Nahm construction for the SU{n) case pOq |. Then, 
one needs to fix asymptotics of the n components of the scalar field to be some (ordered) numbers 
/ii, i = I, ...,n, X^jMi = 0- It naturally divides the domain where the Nahm system is given onto the 
n — 1 intervals (/ii,/i2), On each a-th interval the Nahm system for the three pk x 

matrices Ti^a{t) is 

^=2e,,feK-„r,J (8.16) 

Thus, the SU{n) monopole is characterized by the set of topological charges pi, ...,p„_i. 

The correspondence between solutions to the Nahm equations and the spectral curves of integrable 



systems is generalized to the SU{Nc) case [206 , where the connection with the higher spin magnets 
emerges . The solution to the Nahm equations is a product of n — 1 components, each of them 
being a spectral curve for the SU{2) monopole of charge pa, while the whole monopole spectral curve 
is a special case of the curve for the sl{p) spin magnet with p =^ X^aPa- "^^^ configuration of the 
magnet is specified by a peculiar pair clasterization, each of them describing the sl{pa) magnet. There 
are also additional matching conditions for the correlations between the two-site groups. 

The peculiarity of the spin chain configuration has a clear explanation in the brane language. 



The monopole moduli space is described via the configuration of Dl- and D3-branes [206|. The 
clasterization in the brane language means that the Dl-branes have their endpoints only on the 
neighbor D3-branes. Thus, appealing to the spin chain interpretation, one can identify the D3-branes 
with the sites of the spin chain, while the Dl-branes give rise to the SU{p) group, p being the full 
number of the Dl-branes. Let us note that the twisting procedure used in section 3 is equivalent to 
adding semi-infinite Dl-branes to the very left and the very right D3-branes. 



8.5 Branes and two Lax representations 

Now let us consider two distinct Lax representations of the Toda chain that we discussed in section 3.2 
and explain their interpretation in terms of brane diagrams [ p8[ . Remind that in the 2x2 represen- 
tation the system is formulated in terms of the transfer matrix, while for the x Nc representation 
the notion of the transfer matrix is missed. 

From the brane pictures point of view the two Lax representations are nothing but the choice of 
parametrization of the fivebrane worldvolume E. Of two holomorphic coordinates (which are Cartan 
elements of scalars) A = X4 -I- 1x5 and s = xq + ixio, we can regard any we wish as a spectral 
parameter. In the "standard" 2x2 representation, w — exp ^ is a spectral parameter. Another 
possibility is to choose A as a spectral parameter. It leads to the Nc x Nc representation. 

Geometrical image of this phenomenon is the base-fiber symmetry of the Calabi-Yau threefold 
[ p4| . In the "standard" 2x2 language w parameterizes the base Am singularity and counts the 
number m of gauge factors. A is a coordinate in the fiber. If each coefficient at is a polynomial 
of degree n, we deal with the SU{n)"^ gauge theory with n extra fundamentals at each end of the 
chain of the gauge factors. The duality SU(n)'^ ^ SU{m + 1)"^^ proposed in is nothing but 
the TO+lxTO+l^nxn symmetry of the Lax representations. Under this symmetry the role of 
coordinates on the Newton polygon describing the manifold is changed. 

Note that the generating differential 

^^Equivalcntly, s = + ixg in the Type IIB language. 
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dS = X- 



dw 



Xds 



-sdX 
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also suggests this symmetry. 

Let us extend this argumentation onto description of the fundamental matter in 4(i gauge theory. 
There are two different realizations of the matter in the brane language: via semi-infinite D4-branes 
and via D6-branes. On the integrable side, the first possibility corresponds to the additionally twisted 



Toda chain |20£], while the second one to the XXX magnet. 

Indeed, the curve for the theory with the fundamental matter can be also obtained from the Nc x Nc 
representation by deforming the Lax operator. To this end, one has to add non-trivial entries to the 
first and the last columns which corresponds to attaching the semi-infinite branes to the end branes of 
the configuration. On the contrary, in 2 x 2 representation D6-branes correspond to the fundamentals 
which are localized at the points corresponding to the masses. It is in agreement with the relation 
between the masses and inhomogeneities in the spin chain. The both representations can be related 
if one takes into account the phenomenon of the brane creation when branes pass through each other 

[ill- 

In the authors considered the three dimensional theory in the common longitudinal directions 
of branes (in the type IIB theory): 
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(8.18) 



+ 



By T-duality in the third direction, this setup gets mapped onto the Type llA brane picture pl| ] 
so that the statement of creating a new D3-brane when the D5-brane passes through the NS5-brane 
transforms into the claim that, moving the D6-brane through the NS5-brane creates a new D4-brane 
stretched between them. 



8.6 Whitham dynamics via branes 

So far we discussed the brane interpretation of the Toda/spin chains. Now we are going to say some 
words about the other, Whitham integrability of Seiberg-Witten theory and start with the brane 
configuration in F-theory, i.e. on the elliptically fibered K3 manifold or, equivalently, on oricntifold 
in IIB theory. 

According to | plO| , Af = 2 theory in Ad can be considered as a theory on the worldvolume of 3- 
branes in the background of four 7-branes located at the points ±A on the complex plane of it = Tri^^ 
for the SU (2) gauge group. We assume that masses of all fundamentals are large and, therefore, we 
can work with the pure gauge theory. 

Consider now the dynamics of the 3-branes along the directions transverse to the background 
7-branes. We shall restrict ourselves to the SU{2) gauge group. However, since the whole picture is 
very transparent, the generalization to higher rank groups is straightforward. 

We are going to study the Whitham dynamics of moduli parametrized by branching points. In 
the SU{2) case, one can write the Seiberg-Witten curve in the hyperelliptic form ( ^.8[ ) 



-I o-i 



Y^ = iX^-ur-At,cD, dS = 2X'^ (8.19) 



i.e. there are 4 branching points zt^uzt Aq^j-,. Since the Whitham dynamics governs the u- 
behaviour, all the four points move under the Whitham flows. Therefore, it is more convenient 
to use another parameterization that goes back to the original paper M 



dx 

Y 

This parameterization is obtained from ( 8.19| ) by change of variables x = \^ — u. This is a different 



^{x + u){x-klcD)i^ + KlcD), dS^{x + u)— (8.20) 



curve but the new differential dS has on this curve the same periods as the differential (S.19) 



The parameterization ( B.2C ) has the evident advantage that only one of its branching points is 



moved by the Whitham flows. This Whitham dynamics is nothing but the Gurevich-Pitaevskii solution 



first emerged in the problem of modulation of the one-zone solution of the KdV equation |211]. In 
the brane language this solution looks as follows. 

Initially, the 3-branes coincide with a pair of 7-branes and then they scatter during the evolution 
process. After all, they are absorbed by the remaining pair of 7-branes. Therefore, the exact metric 
on the moduli space is provided by the 3-brane exchange process which generalizes the closed string 



exchange picture suggested in 212 



Returning back to the Whitham dynamics, the coordinate u defines positions of the two 3-branes 
located at the points ±^/u on the a-plancQ while the other branching points fix locations of the back- 
ground branes. The branching point corresponding to the dynamical 3-branes moves in accordance 
with the Whitham flows for the one-zone solution from u = Ag,^^ to m = — Aq(^^, while 7 branes 
remain at rest. 

Let us emphasize the close analogy between the dynamics of branes and the low-energy monopole 
scattering. Indeed, there is a geodesies which corresponds to the transformation of the monopoles into 



dyons |213 . The brane evolution above is equivalent to this process in the following sense. At the 
initial moment, when branes are located at the points ±Aq(^^, there are massless monopoles in the 



theory in accordance with the exact formula (O) for the BPS spectrum. After the scattering process 
at the right angle at m = 0, the branes acquire an electric charge which can be explained as follows. At 
the endpoint of the evolution when branes are at u = — Aq^^, there are massless dyons in the theory. 
However, the dyons can be treated as (1,1) strings stretched between the dynamical and background 



branes |214|. Therefore, the dynamical branes have to acquire the electric charge necessarily. 



9 Conclusion 

We have tried to clarify the role and interpretation of classical integrability in the vacuum sector of 
supersymmetric gauge theories which we basically understand as obtained via string compactifica- 
tions. Briefly this can be formulated as follows. The vacuum sector is described in terms of moduli 
spaces which describe the phase spaces of two different integrable systems. An integrable many-body 
system defines a dynamics on a moduli space which can be formulated in terms of the geometry of 
the compactification manifold. On the other hand, the second integrable dynamics, the Whitham 
dynamics can be interpreted as governing the renormalization group flows in the theory and develops 

Since the whole our consideration here is done at the quasi-classical (perturbative) level, u = o? . 
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on the Coulomb branch of the moduh space. Both integrable systems can be interpreted in terms of 
the dynamics of collective coordinates of branes in different dimensions. 

The effectiveness and the validity of the approach is out of doubts at the moment, despite the lack 
of its derivation from the first principles. The most important problem to be solved seems to be to 
sum up the infinite number of instantons into the finite number of degrees of freedom that provide the 

phase space for the integrable system. One more question essential from the conceptual point of view 
is to clarify the meaning of quantization of the integrable system. Role of the corresponding "Planck 
constant" is uncertain at the moment. 

It is absolutely clear that the branes are extremely important for dealing with nonperturbative 
contributions and play the role analogous to the Feynman diagrams in the perturbation theory. We 

have shown that, in simplest situations, integrability can be described in brane terms, however, the 
extensive formulation of these relations is still missed and, therefore, further work in this direction is 
desired. 
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